Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,976,399 C→A 8.0% A41E (GCA→GAA)  wzxE → lipid III flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,976,3990CA 8.0% 500.9 / 20.6 150A41E (GCA→GAA) wzxElipid III flippase
Reads supporting (aligned to +/- strand):  ref base C (72/66);  new base A (8/4);  total (80/70)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 161 total aligned reads displayed.

GCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTC  >  NC_000913/3976250‑3976547
                                                                                                                                                     |                                                                                                                                                    
gCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcggc                                                                                                                                                      >  2:1980541/1‑150 (MQ=255)
 cGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                                                                     >  2:61159/1‑150 (MQ=255)
 cGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                                                                     >  2:70076/1‑150 (MQ=255)
       tgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTc                                                                                                                                               >  2:1549321/1‑150 (MQ=255)
        gttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTcc                                                                                                                                              <  1:181086/150‑1 (MQ=255)
         ttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCg                                                                                                                                             <  1:2426363/150‑1 (MQ=255)
           gAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGcc                                                                                                                                           >  1:2395704/1‑150 (MQ=255)
           gAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGcc                                                                                                                                           >  2:2289677/1‑150 (MQ=255)
             actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                                                                         >  2:838891/1‑150 (MQ=255)
              ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgc                                                                                                                                        <  1:740343/150‑2 (MQ=255)
              ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgc                                                                                                                                        <  1:1189632/150‑2 (MQ=255)
               tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGtt                                                                                                                                       <  2:1434188/150‑1 (MQ=255)
                actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgccc                                                                                                                                      <  2:668965/150‑4 (MQ=255)
                        tGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGt                                                                                                                              <  1:1378025/150‑1 (MQ=255)
                        tGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGt                                                                                                                              <  1:803864/150‑1 (MQ=255)
                           tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACAATGGTCAAGATTGGTGCCGGGTTACTGGTCGTTAAGTTGCTGGCGGTGTCATATGGTCCGGCGGCGCTTCGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCt                                                                                                                           >  1:593862/1‑150 (MQ=255)
                            atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                                                                          <  1:971999/150‑1 (MQ=255)
                             tGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCg                                                                                                                         >  1:1136377/1‑150 (MQ=255)
                             tGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCg                                                                                                                         <  2:154109/150‑1 (MQ=255)
                              gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCgg                                                                                                                        >  1:490443/1‑150 (MQ=255)
                                     gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                                                                                 >  2:981326/1‑150 (MQ=255)
                                     gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                                                                                 >  1:1310920/1‑150 (MQ=255)
                                      cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                                                                                >  2:1233985/1‑150 (MQ=255)
                                        aaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgc                                                                                                                                        <  1:1455366/124‑2 (MQ=255)
                                          aGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGTTTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCGGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                                                                                            <  1:653727/150‑1 (MQ=255)
                                           gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                                                                           >  2:469771/1‑150 (MQ=255)
                                                cTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTgg                                                                                                      <  2:29678/150‑1 (MQ=255)
                                                 ttGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTGGCCGGGGCTGGc                                                                                                     <  1:300873/150‑1 (MQ=255)
                                                   gtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCTGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                                                                                   <  2:1129708/150‑1 (MQ=255)
                                                    gggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                                                                                  <  2:1737093/149‑1 (MQ=255)
                                                     ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                                                                                 >  2:2452292/1‑150 (MQ=255)
                                                        cggcggCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTa                                                                                              >  2:1956171/1‑150 (MQ=255)
                                                          gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                                                                            <  1:32987/150‑1 (MQ=255)
                                                            ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                                                                          >  2:559798/1‑150 (MQ=255)
                                                             gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                                                                         >  1:401957/1‑150 (MQ=255)
                                                             gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                                                                         <  2:2206665/150‑1 (MQ=255)
                                                              ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                        <  1:1082497/150‑1 (MQ=255)
                                                              ccAGTACACTAGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                        >  2:258206/1‑150 (MQ=255)
                                                               ctgTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCGGGCTGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtgt                                                                                       <  1:1547438/148‑1 (MQ=255)
                                                                  tACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                                                                                    <  2:158919/150‑1 (MQ=255)
                                                                    cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                  >  1:2240406/1‑150 (MQ=255)
                                                                    cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                  <  1:1783013/150‑1 (MQ=255)
                                                                    cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                  >  1:287877/1‑150 (MQ=255)
                                                                       tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                                                                               >  1:70795/1‑150 (MQ=255)
                                                                        ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                              >  1:767036/1‑150 (MQ=255)
                                                                        ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                              >  2:366075/1‑150 (MQ=255)
                                                                        ggTCAAGATTGGGGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                              <  1:2487545/150‑1 (MQ=255)
                                                                          tCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACg                                                                            >  1:661496/1‑150 (MQ=255)
                                                                           cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                                                                           >  2:2100293/1‑150 (MQ=255)
                                                                           cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCGTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                                                                           <  2:1472509/150‑1 (MQ=255)
                                                                              gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                                                                        >  2:968842/1‑150 (MQ=255)
                                                                               ctttgTGCCGGGGTACTGGTCGGTAATTTGCTGGCGGGGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAATTCCGCCATTTGTTTACCGTGCTCGGCGTGCCTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGcc                                                                       <  2:1730913/146‑1 (MQ=255)
                                                                                  ggTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg                                                                    >  2:1345923/1‑150 (MQ=255)
                                                                                   gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt                                                                   >  1:432826/1‑150 (MQ=255)
                                                                                   gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt                                                                   <  2:2038964/150‑1 (MQ=255)
                                                                                    tGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa                                                                  <  2:2180691/150‑1 (MQ=255)
                                                                                    ggCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa                                                                  <  1:580202/149‑1 (MQ=255)
                                                                                     gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCTGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc                                                                 <  2:1776832/150‑1 (MQ=255)
                                                                                        gggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCAt                                                              >  1:565404/1‑150 (MQ=255)
                                                                                         ggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATg                                                             >  1:1574926/1‑150 (MQ=255)
                                                                                          gTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGGCCGGCGGGGCGTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGa                                                            <  1:1408753/150‑1 (MQ=255)
                                                                                           ttACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt                                                           <  1:258206/150‑1 (MQ=255)
                                                                                                ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                      <  2:1621162/150‑1 (MQ=255)
                                                                                                    ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCaa                                                  >  1:1631903/1‑150 (MQ=255)
                                                                                                      tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACa                                                >  1:962161/1‑150 (MQ=255)
                                                                                                        aGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGc                                              >  1:486061/1‑150 (MQ=255)
                                                                                                         gTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCt                                             >  2:2356232/1‑150 (MQ=255)
                                                                                                         gTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCt                                             >  2:2438860/1‑150 (MQ=255)
                                                                                                             cTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgcgc                                         >  2:1314901/1‑150 (MQ=255)
                                                                                                               ggCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGATGCCCAGTACCATGATAATCCGCAACAGCTGcgccg                                       >  1:628544/1‑150 (MQ=255)
                                                                                                                 cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCcgcg                                     >  1:856803/1‑150 (MQ=255)
                                                                                                                 cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCcgcg                                     <  2:1631903/150‑1 (MQ=255)
                                                                                                                  ggTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGt                                    >  2:1927261/1‑150 (MQ=255)
                                                                                                                      tCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTc                                <  1:2289677/150‑1 (MQ=255)
                                                                                                                      tCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTc                                <  2:1063405/150‑1 (MQ=255)
                                                                                                                       cATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCg                               <  2:1113920/150‑1 (MQ=255)
                                                                                                                        aTTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCgg                              >  2:2071659/1‑150 (MQ=255)
                                                                                                                           tGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCAc                           >  2:574940/1‑150 (MQ=255)
                                                                                                                             gTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACtt                         >  1:1046452/1‑150 (MQ=255)
                                                                                                                               ccGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTtca                       >  1:207476/1‑150 (MQ=255)
                                                                                                                                cggcggGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTtcat                      <  2:755102/150‑1 (MQ=255)
                                                                                                                                    ggggCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGc                  <  2:253194/150‑1 (MQ=255)
                                                                                                                                       gCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGAt               <  2:962161/150‑1 (MQ=255)
                                                                                                                                         ttGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATgg             <  1:1701131/150‑1 (MQ=255)
                                                                                                                                          tGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGt            <  1:2071659/150‑1 (MQ=255)
                                                                                                                                           gggCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTAGGCACTTCATCAGCGATGGTa           >  2:1006580/1‑150 (MQ=255)
                                                                                                                                            ggCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                    <  2:1755480/108‑1 (MQ=255)
                                                                                                                                            ggCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                    >  1:1755480/1‑108 (MQ=255)
                                                                                                                                            ggCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTAc          >  1:1687752/1‑150 (MQ=255)
                                                                                                                                             gCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt         >  2:1062296/1‑150 (MQ=255)
                                                                                                                                             acTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt         >  2:2073634/2‑150 (MQ=255)
                                                                                                                                                ggcggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTgg      <  1:1233985/150‑1 (MQ=255)
                                                                                                                                                 gcggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGt     <  1:1980541/150‑1 (MQ=255)
                                                                                                                                                 gcggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGt     >  1:2179669/1‑150 (MQ=255)
                                                                                                                                                 gcggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGt     <  1:2438860/150‑1 (MQ=255)
                                                                                                                                                  cggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGtt    >  1:1581487/1‑150 (MQ=255)
                                                                                                                                                  cggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGtt    <  2:1396694/150‑1 (MQ=255)
                                                                                                                                                  cggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGtt    <  2:767036/150‑1 (MQ=255)
                                                                                                                                                   ggcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGttt   <  2:2149303/150‑1 (MQ=255)
                                                                                                                                                    gcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCAATGGTACTTGGTTtc  <  1:1638351/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTC  >  NC_000913/3976250‑3976547

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