Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,110 G→A 5.4% T26I (ACC→ATC)  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,1100GA 5.4% 681.3 / 18.4 203T26I (ACC→ATC) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base G (87/105);  new base A (6/5);  total (93/110)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

Only 100 of 206 total aligned reads displayed.

AGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTA  >  NC_000913/671964‑672258
                                                                                                                                                  |                                                                                                                                                     
aGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTa                                                                                                                                                     <  1:1066323/149‑1 (MQ=255)
 ggATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTAt                                                                                                                                                    >  1:105516/1‑149 (MQ=255)
         cGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGa                                                                                                                                            <  1:1550533/149‑1 (MQ=255)
          gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                           >  3:1333047/1‑149 (MQ=255)
          gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                           >  1:1219220/1‑149 (MQ=255)
              ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       >  4:729095/1‑149 (MQ=255)
                aCGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagc                                                                                                                                     <  3:78400/149‑1 (MQ=255)
                aCGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagc                                                                                                                                     <  1:1389758/149‑1 (MQ=255)
                  gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       <  2:1264183/145‑1 (MQ=255)
                  gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       >  1:1264183/1‑145 (MQ=255)
                       gTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAccc                                                                                                                              >  4:1558590/1‑149 (MQ=255)
                          tCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggc                                                                                                                           >  1:1382803/1‑149 (MQ=255)
                          tCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggc                                                                                                                           >  3:80177/1‑149 (MQ=255)
                           cTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcg                                                                                                                          >  4:88527/1‑149 (MQ=255)
                            tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                         >  1:1123642/1‑149 (MQ=255)
                             ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                        <  1:723997/149‑1 (MQ=255)
                             ttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGGGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                        <  4:1029299/149‑1 (MQ=255)
                               aTCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      <  3:1558590/149‑1 (MQ=255)
                                tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAt                                                                                                                     <  2:1123642/149‑1 (MQ=255)
                                 cAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                    <  2:498672/149‑1 (MQ=255)
                                  aaGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                             >  3:1124443/1‑123 (MQ=255)
                                      aaCTCGACACCATTCAGACGTAACTTGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAa‑cac                                                                                                              >  3:1170287/1‑149 (MQ=255)
                                            acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTG‑TTAATCACCGccgg                                                                                                        <  3:402109/149‑3 (MQ=255)
                                              acCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTGTTAAT‑CACCGccgggt                                                                                                      >  1:1343248/1‑145 (MQ=255)
                                                cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTG‑TTAATCACCGccgggtgt                                                                                                    <  2:1343248/149‑7 (MQ=255)
                                                cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGAt                                                                                                    <  3:662738/149‑1 (MQ=255)
                                                 aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCAGAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATa                                                                                                   >  3:1066069/1‑149 (MQ=255)
                                                   tCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaac                                                                                                 <  2:975181/149‑1 (MQ=255)
                                                            aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  1:49975/149‑1 (MQ=255)
                                                            aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  4:1333047/149‑1 (MQ=255)
                                                            aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        <  3:196162/149‑1 (MQ=255)
                                                             cTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTg                                                                                       <  1:435439/149‑1 (MQ=255)
                                                              tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGc                                                                                      <  2:88019/149‑1 (MQ=255)
                                                               ggTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGcc                                                                                     <  4:225954/149‑1 (MQ=255)
                                                                gTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                        >  1:1208857/1‑114 (MQ=255)
                                                                   aCGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGat                                                                                 <  1:765957/149‑1 (MQ=255)
                                                                    cgcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      >  3:600185/1‑112 (MQ=255)
                                                                    cgcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGata                                                                                <  2:1188293/149‑1 (MQ=255)
                                                                     gcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGatat                                                                               <  2:1219220/149‑1 (MQ=255)
                                                                      caccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                              <  3:1399597/149‑1 (MQ=255)
                                                                       accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCg                                                                             >  1:1157610/1‑149 (MQ=255)
                                                                       accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCg                                                                             >  1:235462/1‑149 (MQ=255)
                                                                       accgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCg                                                                             >  1:895847/1‑149 (MQ=255)
                                                                         cgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCa                                                                           <  3:498151/149‑1 (MQ=255)
                                                                          gcacGGCTTAATGGCCCGTTCGGGTCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCAc                                                                          <  4:164852/149‑1 (MQ=255)
                                                                              ggCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCtt                                                                      >  4:1439245/1‑149 (MQ=255)
                                                                                 ttAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcct                                                                   >  1:1540366/1‑149 (MQ=255)
                                                                                       gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGat                                                                                 >  3:153303/1‑129 (MQ=255)
                                                                                       gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                             >  2:649183/1‑149 (MQ=255)
                                                                                       gCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcgcg                                                             >  2:959154/1‑149 (MQ=255)
                                                                                        cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CAcc                                                                                                             >  1:375765/1‑100 (MQ=255)
                                                                                        cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CAcc                                                                                                             <  2:375765/100‑1 (MQ=255)
                                                                                        cccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCgcgc                                                            <  2:837111/149‑1 (MQ=255)
                                                                                         ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgcca                                                                                                                                       <  2:117086/74‑1 (MQ=255)
                                                                                         ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCt                                                           >  4:133710/1‑149 (MQ=255)
                                                                                         ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCt                                                           >  4:1062152/1‑149 (MQ=255)
                                                                                         ccGGTCGAATCGCCTGAGTCCATGATCACGCCCTTCATAGTGGAACGAACCTGCGTGGTATCAGGCAGATGCCAGTCACACCCGGCGGTCATTAA‑CACCGCCAGAGATAACGACAATGTTGCCAGATAGCGCACGCTCGCTCCCACGCt                                                           <  2:830701/149‑1 (MQ=255)
                                                                                           gTTCGGATCGCCTGAGTCCAGGATCATGATCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTCCCTCCCGCGCTTa                                                         <  4:1102001/149‑1 (MQ=255)
                                                                                           gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                         <  3:1062152/149‑1 (MQ=255)
                                                                                            ttCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAg                                                        >  2:1390906/1‑149 (MQ=255)
                                                                                               ggAGCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCa                                                     >  2:274411/1‑149 (MQ=255)
                                                                                                gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                    <  3:913478/149‑1 (MQ=255)
                                                                                                 aTCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGt                                                                                                                        >  1:117905/1‑81 (MQ=255)
                                                                                                  tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                      >  1:124619/1‑82 (MQ=255)
                                                                                                  tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGc                                                       >  1:816964/1‑144 (MQ=255)
                                                                                                    gCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGAc                                                <  2:151123/149‑1 (MQ=255)
                                                                                                      cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                              <  1:1028185/149‑1 (MQ=255)
                                                                                                      cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                              >  4:1527099/1‑149 (MQ=255)
                                                                                                      cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCa                                              >  3:1312251/1‑149 (MQ=255)
                                                                                                          gTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAtt                                          <  1:610227/149‑1 (MQ=255)
                                                                                                            ccAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTga                                        <  2:1382803/149‑1 (MQ=255)
                                                                                                             cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                     >  2:1424683/1‑119 (MQ=255)
                                                                                                             cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                     <  1:1424683/119‑1 (MQ=255)
                                                                                                             cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                       >  2:172788/1‑149 (MQ=255)
                                                                                                             cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                       >  3:1221410/1‑149 (MQ=255)
                                                                                                               ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                     <  4:1472282/149‑1 (MQ=255)
                                                                                                               ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                     >  4:1142433/1‑149 (MQ=255)
                                                                                                                 aTCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGc                                                                                                           >  3:1246955/1‑77 (MQ=255)
                                                                                                                 aTCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGc                                                                                                           <  4:1246955/77‑1 (MQ=255)
                                                                                                                 aTCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGt                                   <  2:1544217/149‑1 (MQ=255)
                                                                                                                  tCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                  >  3:1444371/1‑149 (MQ=255)
                                                                                                                   cATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                                        >  1:973843/1‑94 (MQ=255)
                                                                                                                      gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGATATTAA‑CACCGCCAGAGATAACAACAATGTTGCCTTATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                              <  2:1095953/149‑1 (MQ=255)
                                                                                                                       aCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                             <  1:833911/149‑1 (MQ=255)
                                                                                                                            cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTAc                        <  4:622232/149‑1 (MQ=255)
                                                                                                                             aTAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACg                       >  3:1363960/1‑149 (MQ=255)
                                                                                                                                 tGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTaaa                   >  2:348951/1‑149 (MQ=255)
                                                                                                                                  ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                  >  4:90049/1‑149 (MQ=255)
                                                                                                                                  ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                  <  2:895847/149‑1 (MQ=255)
                                                                                                                                    aaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCa                >  4:652376/1‑149 (MQ=255)
                                                                                                                                        catcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                                    <  1:557813/106‑1 (MQ=255)
                                                                                                                                        catcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACttt            >  2:1405565/4‑149 (MQ=255)
                                                                                                                                         aCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa           <  2:430334/149‑1 (MQ=255)
                                                                                                                                         aCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa           <  2:804756/149‑1 (MQ=255)
                                                                                                                                           cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACg         >  3:983323/1‑149 (MQ=255)
                                                                                                                                           cTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                  <  2:356501/124‑1 (MQ=255)
                                                                                                                                            tGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                  >  3:516997/1‑59 (MQ=255)
                                                                                                                                            tGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATaa                                                                                                  <  4:516997/59‑1 (MQ=255)
                                                                                                                                            tGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGt        >  1:865597/1‑149 (MQ=255)
                                                                                                                                                tGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAg    <  2:105516/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGGATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTA  >  NC_000913/671964‑672258

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