Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,685,561 T→C intergenic (‑54/‑274) glyA ← / → hmp serine hydroxymethyltransferase/fused nitric oxide dioxygenase/dihydropteridine reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,685,5610TC100.0% 286.1 / NA 78intergenic (‑54/‑274)glyA/hmpserine hydroxymethyltransferase/fused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/41);  total (37/41)

GACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTC  >  NC_000913/2685432‑2685695
                                                                                                                                 |                                                                                                                                      
gACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACaa                                                                                                                               >  1:571197/1‑139 (MQ=255)
  cGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCCTTATAGGCGGTCCCGTTGGACAACg                                                                                                                             <  2:152376/139‑1 (MQ=255)
     aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  2:1667621/139‑1 (MQ=255)
     aCTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTGAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  2:1071506/139‑1 (MQ=255)
       tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        <  2:520846/139‑1 (MQ=255)
       tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        >  2:1704515/1‑139 (MQ=255)
       tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        >  1:1297950/1‑139 (MQ=255)
       tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        <  1:1704515/139‑1 (MQ=255)
       tttttCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        <  1:1616820/139‑1 (MQ=255)
               gCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATaa                                                                                                                >  1:68719/1‑139 (MQ=255)
                        ccTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCAt                                                                                                       <  2:1050508/139‑1 (MQ=255)
                          tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        >  2:1410663/1‑120 (MQ=255)
                          tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        <  1:1410663/120‑1 (MQ=255)
                          tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                     >  1:1638212/1‑139 (MQ=255)
                          tGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtg                                                                                                     >  2:1604867/1‑139 (MQ=255)
                            ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                  <  2:1120463/108‑1 (MQ=255)
                            ccACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                   <  2:571197/139‑1 (MQ=255)
                            acacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                  >  1:1120463/2‑108 (MQ=255)
                             cacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                  <  1:1064864/107‑1 (MQ=255)
                             cacaGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGa                                                                                                                                  >  2:1064864/1‑107 (MQ=255)
                                  ttCGGCTTCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          >  1:1576067/1‑110 (MQ=255)
                                  ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCg                                                                                                                          <  2:1576067/110‑1 (MQ=255)
                                  ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             <  2:1638212/139‑1 (MQ=255)
                                  ttCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             <  2:68719/139‑1 (MQ=255)
                                                    aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        <  2:323310/94‑1 (MQ=255)
                                                    aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGaa                                                                                                                        >  1:323310/1‑94 (MQ=255)
                                                    aaTGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGaa                                                                           >  2:1694826/1‑139 (MQ=255)
                                                      tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGc                                                                                                  <  1:39156/114‑1 (MQ=255)
                                                      tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGc                                                                                                  >  2:39156/1‑114 (MQ=255)
                                                      tGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAAtt                                                                         >  2:782241/1‑139 (MQ=255)
                                                             ttCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaa                                                                  >  1:857824/1‑139 (MQ=255)
                                                              tCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCaaa                                                                 <  2:857824/139‑1 (MQ=255)
                                                                  gCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                             <  1:782241/139‑1 (MQ=255)
                                                                   cTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAtt                                                            >  2:1783954/1‑139 (MQ=255)
                                                                      tAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTAcc                                                         >  1:456386/1‑139 (MQ=255)
                                                                        aCATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCtt                                                       <  1:1694826/139‑1 (MQ=255)
                                                                            ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   >  1:776291/1‑139 (MQ=255)
                                                                            ccGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   <  2:937182/139‑1 (MQ=255)
                                                                             cGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAAtt                                                                         <  2:591874/116‑1 (MQ=255)
                                                                             cGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAAtt                                                                         >  1:591874/1‑116 (MQ=255)
                                                                                   tcCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAAGAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTAAGCa                                            <  1:1157073/139‑1 (MQ=255)
                                                                                   tcCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCa                                            >  1:454598/1‑139 (MQ=255)
                                                                                     cTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAg                                          >  1:794128/1‑139 (MQ=255)
                                                                                      tGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   <  2:1296809/129‑1 (MQ=255)
                                                                                      tGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   >  1:1296809/1‑129 (MQ=255)
                                                                                          tCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                   >  2:1522985/1‑77 (MQ=255)
                                                                                          tCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgtg                                                                                                   <  1:1522985/77‑1 (MQ=255)
                                                                                          tCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCttt                                     <  1:1783954/139‑1 (MQ=255)
                                                                                           cAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTg                                    <  2:1297950/139‑1 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             >  2:425191/1‑81 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATaa                                                                                             <  1:425191/81‑1 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                             <  2:396121/113‑1 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                             >  1:396121/1‑113 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   <  2:1569578/123‑1 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   >  1:1569578/1‑123 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                   <  2:794128/139‑1 (MQ=255)
                                                                                            aGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                   <  2:454598/139‑1 (MQ=255)
                                                                                              cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGga                                  <  2:1723954/138‑1 (MQ=255)
                                                                                              cTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGga                                  >  1:1723954/1‑138 (MQ=255)
                                                                                               tAACAATAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGgaga                                >  2:648610/1‑139 (MQ=255)
                                                                                                   aaTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgt                                                                                                    <  2:620945/67‑1 (MQ=255)
                                                                                                   aaTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATtgt                                                                                                    >  1:620945/1‑67 (MQ=255)
                                                                                                    aTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCAt                                                                                                       <  1:600371/63‑1 (MQ=255)
                                                                                                    aTAAAATTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCAt                                                                                                       >  2:600371/1‑63 (MQ=255)
                                                                                                       aaaTTTTTGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAagag                        <  1:648610/139‑1 (MQ=255)
                                                                                                             ttGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAATCGAATCATTGTGCGATAACAGGTCTTGACaaa                                                                               <  2:103473/78‑1 (MQ=255)
                                                                                                             ttGGCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAATCGAATCATTGTGCGATAACAGGTCTTGACaaa                                                                               >  1:103473/1‑78 (MQ=255)
                                                                                                               ggCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACa                                                                                 <  1:644422/74‑1 (MQ=255)
                                                                                                               ggCCTTTATAGGCGGTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACa                                                                                 >  2:644422/1‑74 (MQ=255)
                                                                                                                            ggTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGt                                                                                        >  1:351235/1‑54 (MQ=255)
                                                                                                                            ggTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGt                                                                                        <  2:351235/54‑1 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                             >  2:389176/1‑80 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGAt                                                             <  1:389176/80‑1 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   >  2:533570/1‑90 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAgg                                                   <  1:533570/90‑1 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                   <  1:1369980/106‑1 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTgg                                   >  2:1369980/1‑106 (MQ=255)
                                                                                                                             gTCCCGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTc  >  1:599402/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                      
GACGTACTTTTTCCTGCTCCATAGCCTGCCACAGTTCGGCATCATAATCGGCAATGTTCATTTCACGCTTTAACATCCGCATCTCCTGACTCAGCTAACAATAAAATTTTTGGCCTTTATAGGCGGTCCTGTTGGACAACGGCGAACAGTATAACCGAATCATTGTGCGATAACAGGTCTTGACAAAGGAATTTACGCAAACGATTACCTTCAGGCTACGCAAGGCTTTGGAGAATAAAGAGCTTGCAACCGGAAACGGATTTC  >  NC_000913/2685432‑2685695

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