Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,342,315 G→A 6.1% S132N (AGC→AAC)  pyrF → orotidine‑5'‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,342,3150GA 6.1% 212.5 / 3.1 65S132N (AGC→AAC) pyrForotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (27/34);  new base A (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATT  >  NC_000913/1342168‑1342459
                                                                                                                                                   |                                                                                                                                                
aCGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAg                                                                                                                                                  >  4:157729/1‑148 (MQ=255)
      tCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGa                                                                                                                                            <  6:246340/148‑1 (MQ=255)
        gctgctGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGaa                                                                                                                                           <  5:98796/147‑1 (MQ=255)
                 gCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGAcc                                                                                                                                 >  7:242373/1‑148 (MQ=255)
                     aCTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTgg                                                                                                                              <  8:160139/147‑1 (MQ=255)
                     aCTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTg                                                                                                                               <  8:94860/146‑1 (MQ=255)
                     aCTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGt                                                                                                                             <  4:272418/148‑1 (MQ=255)
                               gtgGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTgg                                                                                                                    >  7:305261/1‑147 (MQ=255)
                                tgGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                   <  1:388752/147‑1 (MQ=255)
                                tgGATGGTGAAGGTTCTTGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                   <  7:207288/147‑1 (MQ=255)
                                     ggTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGaca                                                                                                             >  6:293232/1‑148 (MQ=255)
                                      gTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGTCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGacac                                                                                                            <  2:281053/148‑1 (MQ=255)
                                      gTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGaca                                                                                                             <  8:242373/147‑1 (MQ=255)
                                            gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAACATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                   >  4:163212/1‑135 (MQ=255)
                                            gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAACATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                   <  3:157579/135‑1 (MQ=255)
                                            gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAACATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                   <  3:163212/135‑1 (MQ=255)
                                            gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAACATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                   >  4:157579/1‑135 (MQ=255)
                                               cATGCCTCTGGTGGGGCGCGTATGATGACCGCAGTGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGAGACCTGGTCGATCTTGGCATGACACTGTCACCTg                                                                                                   >  6:209837/1‑148 (MQ=255)
                                               cATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCa                                                                                                       <  6:329059/144‑1 (MQ=255)
                                               cATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCa                                                                                                       >  5:329059/1‑144 (MQ=255)
                                                  gCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTg                                                                                                   >  4:325026/1‑145 (MQ=255)
                                                   ccTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAg                                                                                                >  5:88501/1‑147 (MQ=255)
                                                         gtggGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTAt                                                                                           <  4:349497/146‑1 (MQ=255)
                                                         gtggGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTAt                                                                                           <  3:325026/146‑1 (MQ=255)
                                                         gtggGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATg                                                                                          <  4:152521/147‑1 (MQ=255)
                                                            gggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGTAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGa                                                                                      <  3:336526/148‑1 (MQ=255)
                                                            gggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGa                                                                                      <  7:153732/148‑1 (MQ=255)
                                                            gggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGa                                                                                      <  8:305261/148‑1 (MQ=255)
                                                                 cgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTc                                                                                 <  6:260522/148‑1 (MQ=255)
                                                                    atgatgACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGacac                                                                                                            <  2:23334/118‑1 (MQ=255)
                                                                    atgatgACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGacac                                                                                                            >  1:23334/1‑118 (MQ=255)
                                                                       atgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcgg                                                                           >  3:234540/1‑148 (MQ=255)
                                                                       atgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcgg                                                                           >  2:392631/1‑148 (MQ=255)
                                                                                    gTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCaaaa                                                              <  5:55781/148‑1 (MQ=255)
                                                                                             tGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCCGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGCACGTCTGGAGGCACTGACGCAAAAATGTGGcgt                                                     >  6:187581/1‑146 (MQ=255)
                                                                                             tGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCt                                                     >  8:383764/1‑148 (MQ=255)
                                                                                              ggTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCtt                                                    >  3:280482/1‑148 (MQ=255)
                                                                                              ggTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCtt                                                    >  4:136291/1‑148 (MQ=255)
                                                                                                ttCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGa                                                  >  2:78283/1‑148 (MQ=255)
                                                                                                ttCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGa                                                  >  1:38926/1‑148 (MQ=255)
                                                                                                  ccGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCTTGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAt                                                 <  4:163447/147‑1 (MQ=255)
                                                                                                  ccGTTTGGCAAAGATGCACCGCGTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtg                                                >  7:56889/1‑148 (MQ=255)
                                                                                                       tGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgt                                            <  8:232846/147‑1 (MQ=255)
                                                                                                        ggCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtg                                           <  2:38926/147‑1 (MQ=255)
                                                                                                         gCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtgg                                          <  5:25455/147‑1 (MQ=255)
                                                                                                          cAAAGATGCACCGCTTGTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGTCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTTGCCTTGAtggtgtggtg                                        <  6:83427/148‑1 (MQ=255)
                                                                                                            aaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGCCCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGgtgtggtgtg                                      >  7:308077/1‑148 (MQ=255)
                                                                                                             aGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtg                                                           >  1:297834/1‑126 (MQ=255)
                                                                                                             aGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAAtg                                                           <  2:297834/126‑1 (MQ=255)
                                                                                                                      gCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAg                            >  7:10187/1‑148 (MQ=255)
                                                                                                                              ttGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgt                    >  8:152654/1‑148 (MQ=255)
                                                                                                                                   gtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt               >  1:290870/1‑148 (MQ=255)
                                                                                                                                   gtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt               >  6:36408/1‑148 (MQ=255)
                                                                                                                                    tgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt               <  1:392631/147‑1 (MQ=255)
                                                                                                                                     gACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCtt              >  4:95431/1‑147 (MQ=255)
                                                                                                                                      aCAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCttt             <  3:157729/147‑1 (MQ=255)
                                                                                                                                      aCAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTa            <  3:136291/148‑1 (MQ=255)
                                                                                                                                       cAGTGTGGACCAGCATGGAAGCCAGGGTCCTGGGCGATCTGGGCATGTGACTGTCACGTGCAGATTATGCAGAACGGCTGGCGGCACTGACGCAAAATTGTGGCCTTGATGGTGTGGGGTGTTCTGCTGAGGAAGCTGGGCGCTTTa            <  6:117201/147‑1 (MQ=255)
                                                                                                                                         gtgtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTaaa          >  3:154697/1‑147 (MQ=255)
                                                                                                                                           gtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtg                                                <  7:210568/107‑1 (MQ=255)
                                                                                                                                           gtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtg                                                >  8:210568/1‑107 (MQ=255)
                                                                                                                                           gtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACa        <  1:78283/147‑1 (MQ=255)
                                                                                                                                               aCCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAt   <  7:383764/148‑1 (MQ=255)
                                                                                                                                                ccAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAtt  <  4:234540/148‑1 (MQ=255)
                                                                                                                                                  aGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcg                                                                            <  5:286787/72‑1 (MQ=255)
                                                                                                                                                  aGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcg                                                                            >  6:286787/1‑72 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATT  >  NC_000913/1342168‑1342459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: