Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,113 T→C 5.3% D25G (GAT→GGT)  lptE ← LPS assembly OM complex LptDE, lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913672,1130TC 5.3% 662.3 / 17.0 203D25G (GAT→GGT) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base T (88/104);  new base C (6/5);  total (94/109)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 206 total aligned reads displayed.

GGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTA  >  NC_000913/671965‑672258
                                                                                                                                                    |                                                                                                                                                  
ggATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTAt                                                                                                                                                    >  1:105516/1‑149 (MQ=255)
 gATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATc                                                                                                                                                   <  3:416564/149‑1 (MQ=255)
    ggAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACg                                                                                                                                                >  1:430334/1‑149 (MQ=255)
        cGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAgg                                                                                                                                            <  2:355327/149‑2 (MQ=255)
         gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                           <  2:1147796/149‑1 (MQ=255)
          tCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATg                                                                                                                                          >  2:833911/1‑149 (MQ=255)
             ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       >  1:837111/1‑149 (MQ=255)
             ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       >  4:729095/1‑149 (MQ=255)
                 gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                   >  3:1283270/1‑149 (MQ=255)
                       tttCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcg                                                                                                                             >  1:1544217/1‑149 (MQ=255)
                          cTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcg                                                                                                                          >  4:88527/1‑149 (MQ=255)
                             tATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                       <  1:109459/149‑1 (MQ=255)
                              aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      <  3:1558590/149‑1 (MQ=255)
                              aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      <  4:1445495/149‑1 (MQ=255)
                              aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      <  2:1571634/149‑1 (MQ=255)
                                cAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                    <  2:92256/149‑1 (MQ=255)
                                 aaGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                             <  4:1124443/123‑1 (MQ=255)
                                   gCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑c                                                                                                                >  1:1240247/1‑149 (MQ=255)
                                   gCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑c                                                                                                                >  1:463811/1‑149 (MQ=255)
                                     aaCTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑cac                                                                                                              >  1:735449/1‑149 (MQ=255)
                                       cTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCg                                                                                                            >  2:1391725/1‑149 (MQ=255)
                                           acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTG‑TTAATCACCGccgg                                                                                                        <  3:402109/149‑3 (MQ=255)
                                             acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCagag                                                                                                      >  3:509343/1‑149 (MQ=255)
                                               cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGAt                                                                                                    >  2:1424457/1‑149 (MQ=255)
                                                 ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                  >  4:1334147/1‑149 (MQ=255)
                                                 ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                  >  1:557989/1‑149 (MQ=255)
                                                    aGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaacaa                                                                                               >  1:1478536/1‑149 (MQ=255)
                                                       cGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAAcaacaa                                                                                            >  2:520444/1‑149 (MQ=255)
                                                       cGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAAcaacaa                                                                                            <  3:1065261/149‑1 (MQ=255)
                                                        gTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAAt                                                                                           <  4:1081497/149‑1 (MQ=255)
                                                           aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  1:89860/149‑1 (MQ=255)
                                                           aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  4:1045379/149‑1 (MQ=255)
                                                           aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  3:196162/149‑1 (MQ=255)
                                                             tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGc                                                                                      <  1:19296/149‑1 (MQ=255)
                                                              ggTTACGCTCCGCAGGGCTTAATGGCCCGTTAGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGATCCTGCGTGGTATCACGCCGATGCCAGCCACACCCAGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGcc                                                                                     <  3:971050/149‑1 (MQ=255)
                                                               gTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                        <  2:1208857/114‑1 (MQ=255)
                                                                ttACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAg                                                                                   >  3:1378826/1‑149 (MQ=255)
                                                                   cgcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGata                                                                                <  1:1207961/149‑1 (MQ=255)
                                                                    gcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGatat                                                                               >  2:1107089/1‑149 (MQ=255)
                                                                     caccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                              <  2:894210/149‑1 (MQ=255)
                                                                     caccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                              <  2:1478536/149‑1 (MQ=255)
                                                                      accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCg                                                                             >  2:129598/1‑149 (MQ=255)
                                                                        cgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGAAAGAAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGGTATCGCa                                                                           <  2:434675/149‑1 (MQ=255)
                                                                              gCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                     <  1:828891/149‑1 (MQ=255)
                                                                                ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcct                                                                   <  1:1150636/149‑1 (MQ=255)
                                                                                 tAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcctc                                                                  <  2:424153/149‑1 (MQ=255)
                                                                                     ggCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgc                                                              >  3:1221448/1‑149 (MQ=255)
                                                                                      gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGat                                                                                 >  3:153303/1‑129 (MQ=255)
                                                                                      gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGat                                                                                 <  4:153303/129‑1 (MQ=255)
                                                                                      gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                             >  3:596347/1‑149 (MQ=255)
                                                                                      gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                             >  1:804756/1‑149 (MQ=255)
                                                                                       cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGCAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACACTGTTGCCAGATATCGCACGCTTCCTCCCgcgc                                                            <  3:1324789/149‑1 (MQ=255)
                                                                                        ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       >  1:117086/1‑74 (MQ=255)
                                                                                        ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                         <  2:1047826/88‑1 (MQ=255)
                                                                                        ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCt                                                           >  3:788679/1‑149 (MQ=255)
                                                                                         cGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCtt                                                          <  3:1142922/149‑1 (MQ=255)
                                                                                          gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                         >  2:1371261/1‑149 (MQ=255)
                                                                                          gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                         <  3:133710/149‑1 (MQ=255)
                                                                                           ttCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAg                                                        <  1:649183/149‑1 (MQ=255)
                                                                                            tCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGc                                                       <  4:1270101/149‑1 (MQ=255)
                                                                                               gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                    <  1:520444/149‑1 (MQ=255)
                                                                                               gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                    <  1:1391725/149‑1 (MQ=255)
                                                                                                aTCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                        <  2:117905/81‑1 (MQ=255)
                                                                                                 tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      <  2:124619/82‑1 (MQ=255)
                                                                                                 tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGc                                                       <  2:816964/144‑1 (MQ=255)
                                                                                                 tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACg                                                  >  4:18809/1‑149 (MQ=255)
                                                                                                  cGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGa                                                 <  4:1561042/149‑1 (MQ=255)
                                                                                                     cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                              <  1:1028185/149‑1 (MQ=255)
                                                                                                       gAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGa                                            <  2:292500/149‑1 (MQ=255)
                                                                                                         gTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAtt                                          <  2:1540366/149‑1 (MQ=255)
                                                                                                           ccAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTga                                        <  3:1094594/149‑1 (MQ=255)
                                                                                                            cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                       >  2:367678/1‑149 (MQ=255)
                                                                                                            cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACAACCAGATTgag                                       <  3:891975/149‑1 (MQ=255)
                                                                                                            cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATCAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCACCGACCAGATTgaa                                       <  4:1016049/149‑2 (MQ=255)
                                                                                                             aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                         >  2:1055941/1‑130 (MQ=255)
                                                                                                             aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                      >  3:1199516/1‑149 (MQ=255)
                                                                                                             aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                      >  1:1377135/1‑149 (MQ=255)
                                                                                                              ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                     <  2:235462/149‑1 (MQ=255)
                                                                                                              ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                     <  3:729095/149‑1 (MQ=255)
                                                                                                                 tCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                  >  3:1444371/1‑149 (MQ=255)
                                                                                                                  cATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  2:973843/94‑1 (MQ=255)
                                                                                                                      aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                             <  1:304160/149‑1 (MQ=255)
                                                                                                                      aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                             >  1:1321693/1‑149 (MQ=255)
                                                                                                                        cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                           >  3:632656/1‑149 (MQ=255)
                                                                                                                           cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTAc                        <  4:622232/149‑1 (MQ=255)
                                                                                                                           cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTAc                        <  4:1511786/149‑1 (MQ=255)
                                                                                                                                 ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                  >  4:90049/1‑149 (MQ=255)
                                                                                                                                  gaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATc                 >  1:1256749/1‑149 (MQ=255)
                                                                                                                                   aaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCa                >  4:652376/1‑149 (MQ=255)
                                                                                                                                    aggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCAc               >  4:1242320/1‑149 (MQ=255)
                                                                                                                                     ggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTa                                <  1:121471/131‑1 (MQ=255)
                                                                                                                                       aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACttt            <  2:981315/149‑1 (MQ=255)
                                                                                                                                        atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                     >  2:1384843/3‑123 (MQ=255)
                                                                                                                                        atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                     <  1:1384843/121‑1 (MQ=255)
                                                                                                                                           tGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                  <  4:516997/59‑1 (MQ=255)
                                                                                                                                            gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTa       <  4:1378826/149‑1 (MQ=255)
                                                                                                                                               tGGTATCACGCAGATGTCAGCCACACCCTGCGGTGCTTAA‑CACCGCCAGAGATAACAACAACGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGCTTGAGGAGTCTACCTGGTACGTAAATCACTTTACGTACAg    <  2:735449/149‑1 (MQ=255)
                                                                                                                                                 gTATCACTCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa           >  2:869511/1‑140 (MQ=255)
                                                                                                                                                 gTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTATGTAAATCACTTTACGTACAGTa  >  4:382083/1‑149 (MQ=255)
                                                                                                                                                 gTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa           <  1:869511/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTA  >  NC_000913/671965‑672258

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