Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,536,877 C→G 5.6% A145P (GCA→CCA)  narW ← nitrate reductase 2 (NRZ), delta subunit (assembly subunit)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,536,8770CG 5.6% 471.4 / 14.1 142A145P (GCA→CCA) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base C (63/71);  new base G (3/5);  total (66/76)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 143 total aligned reads displayed.

AAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGG  >  NC_000913/1536728‑1537025
                                                                                                                                                     |                                                                                                                                                     
        ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgtt                                                                                                                                               <  1:1056346/150‑1 (MQ=255)
        ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgtt                                                                                                                                               >  2:1552737/1‑150 (MQ=255)
           gTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGGCATCGCGCGCTTCGCGGTTTACTTGTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGttt                                                                                                                                           <  2:1185022/150‑1 (MQ=255)
           gTTCCGCTTCCCACACTTCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                                                                                                                                            <  1:1212955/150‑1 (MQ=255)
            ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                                           <  1:1266496/150‑1 (MQ=255)
            ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGttt                                                                                                                                           <  2:1796288/150‑1 (MQ=255)
             tCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttg                                                                                                                                          <  1:863982/150‑2 (MQ=255)
                cctTCCCACACTGCATCCAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAg                                                                                                                                       <  2:1250308/149‑1 (MQ=255)
                  ttCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgatgct                                                                                                                                     <  2:753387/150‑7 (MQ=38)
                      cacaCTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGTGACCg                                                                                                                                 >  1:767660/1‑150 (MQ=255)
                      cacaCTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCg                                                                                                                                 >  2:12216/1‑150 (MQ=255)
                      cacaCTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCg                                                                                                                                 <  2:891678/150‑1 (MQ=255)
                            gCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc                                                                                                                           >  2:1273061/1‑150 (MQ=255)
                            gCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc                                                                                                                           >  1:1115342/1‑150 (MQ=255)
                            gCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc                                                                                                                           <  1:1627835/150‑1 (MQ=255)
                              aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                         <  2:1400740/150‑1 (MQ=255)
                              aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                         >  1:1672879/1‑150 (MQ=255)
                              aTCAAGCGCCTGGCGGGTGTCAGCGCGCTCTACGCTGTTTAATTGTTTTGGGACACTGTCACGTGAAAGGGAGCGTCCCGCCAGTGGCAGCAGAGCAGCAACCAACCCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                         <  1:81595/150‑1 (MQ=255)
                                 aaGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGc                                                                                                                      >  2:511094/1‑150 (MQ=255)
                                  aGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCg                                                                                                                     >  1:1223/1‑150 (MQ=255)
                                     gcCTTGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGAAACTGTCACTTGAAAGGGAGCGTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                                  <  1:122579/150‑1 (MQ=255)
                                      ccTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa                                                                                                                 >  1:939700/1‑150 (MQ=255)
                                         ggCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCg                                                                                                              <  2:98635/150‑1 (MQ=255)
                                         ggCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCg                                                                                                              <  1:1208831/150‑1 (MQ=255)
                                          gCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACTCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCgg                                                                                                             <  2:1143121/150‑1 (MQ=255)
                                            gggTGTCACCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcg                                                                                                           <  1:962853/150‑1 (MQ=255)
                                             ggTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcgc                                                                                                          <  1:714462/150‑1 (MQ=255)
                                             ggTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcgc                                                                                                          >  2:1605978/1‑150 (MQ=255)
                                                 tCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACa                                                                                                      >  1:740742/1‑150 (MQ=255)
                                                  cATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACAt                                                                                                     >  1:938837/1‑150 (MQ=255)
                                                      gcgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACTGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAg                                                                                                 <  2:1262432/150‑1 (MQ=255)
                                                      gcgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAg                                                                                                 <  2:1115342/150‑1 (MQ=255)
                                                      gcgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAg                                                                                                 <  1:723227/150‑1 (MQ=255)
                                                       cgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                                                                                <  1:1339590/150‑1 (MQ=255)
                                                           ctTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATc                                                                                            >  2:628494/1‑150 (MQ=255)
                                                              cGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCctt                                                                                         <  1:1422668/150‑1 (MQ=255)
                                                                cTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttct                                                                                       <  1:1348333/150‑1 (MQ=255)
                                                                      aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                                 <  1:1826699/150‑1 (MQ=255)
                                                                       cTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCt                                                                                >  1:1516763/1‑150 (MQ=255)
                                                                         tGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGa                                                                              >  2:382685/1‑150 (MQ=255)
                                                                         tGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGa                                                                              <  2:1401199/150‑1 (MQ=255)
                                                                          ttcaagtgACAGTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGAt                                                                             >  1:1005198/6‑150 (MQ=255)
                                                                          gTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGAt                                                                             >  2:1077019/1‑150 (MQ=255)
                                                                           ttttGTGACAGTGTCACTTGAATGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAAGAACGCGTACCTCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGTTCGGCGCAACTTTGAGCAATCCTTCTTTCGCCTgttt                                                                            <  1:245008/150‑4 (MQ=255)
                                                                           ttttGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATc                                                                            <  1:958133/150‑1 (MQ=255)
                                                                              tgtgACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGTGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTc                                                                         <  1:501116/150‑1 (MQ=255)
                                                                              tgtgACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTc                                                                         <  2:1388397/150‑1 (MQ=255)
                                                                                      tggcACTTGAACGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCGTCCCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTa                                                                 <  2:778186/147‑1 (MQ=255)
                                                                                        tCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACg                                                               <  1:792363/150‑1 (MQ=255)
                                                                                        tCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACg                                                               <  2:1516763/150‑1 (MQ=255)
                                                                                         cACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGc                                                              <  1:1408948/150‑1 (MQ=255)
                                                                                          aCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                             >  2:1306077/1‑150 (MQ=255)
                                                                                          aCTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                             >  2:806260/1‑150 (MQ=255)
                                                                                            ttGAAAGGTAGCTTCCCGCAAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTTAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTa                                                           <  2:1366296/150‑1 (MQ=255)
                                                                                            ttGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTCTTTTAAGCGACCGCCAAGCAGGGCGAGGCTCGGCGCCACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTTCGCGTa                                                           >  1:1686079/1‑150 (MQ=255)
                                                                                              gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                         >  2:1226398/1‑150 (MQ=255)
                                                                                               aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGTATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt                                                        <  1:1767108/150‑1 (MQ=255)
                                                                                               aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt                                                        <  1:1504241/150‑1 (MQ=255)
                                                                                               aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt                                                        >  2:1486566/1‑150 (MQ=255)
                                                                                               aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCGTAAAt                                                        >  2:535882/1‑150 (MQ=255)
                                                                                                aaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATa                                                       >  1:655415/1‑150 (MQ=255)
                                                                                                  gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                                     >  2:1529578/1‑150 (MQ=255)
                                                                                                       cTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGa                                                >  1:528097/1‑150 (MQ=255)
                                                                                                        ttCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGTTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGAt                                               <  1:385890/150‑1 (MQ=255)
                                                                                                        ttCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGCTCGTCAGGCAGTACGCTTAAATACTCCAGAt                                               >  1:54549/1‑150 (MQ=255)
                                                                                                        ttCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGAt                                               <  1:592665/150‑1 (MQ=255)
                                                                                                              ccAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCg                                         <  2:30583/150‑1 (MQ=255)
                                                                                                                 gTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTa                                      <  1:173627/150‑1 (MQ=255)
                                                                                                                 gTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTa                                      >  1:814090/1‑150 (MQ=255)
                                                                                                                 gTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTa                                      <  1:579904/150‑1 (MQ=255)
                                                                                                                  ttGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGTGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAg                                     <  1:540233/150‑1 (MQ=255)
                                                                                                                   tGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCGTTCGCCTGCTCGTCAGGCAGTCCGCTTAAATACGCCAGATACAGCGGTAGa                                   >  2:182464/1‑150 (MQ=255)
                                                                                                                     cagcagAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACAGTGAGCACTCCTTGTTTCGCCTGATCGGCAGGCAGTACGCTTAAATACTCCAGATAGAGCGGTAGATa                                 >  2:1174623/1‑150 (MQ=255)
                                                                                                                         agagCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCg                             >  1:507935/1‑150 (MQ=255)
                                                                                                                         agagCATCAAACAACGCGTACCACGGTGCCACGAGTAGTTTAAAGCGCCCGCCCAGCAGCGCGAGGATCGGCGCAACAGTGAGCAATCCTTATTTCGCCTGATCGACAGGCAGGACGCTTAAATACTCCAGATGACGCGGTAGATAATCg                             >  2:1150128/1‑150 (MQ=255)
                                                                                                                          gagCATCTAACAACGCGTACCTCGGTGGCTCGCGTTGTTTTAAGCGTCCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTAGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCgg                            <  1:365427/150‑1 (MQ=255)
                                                                                                                           agCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGTATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCggg                           <  2:1672879/150‑1 (MQ=255)
                                                                                                                            gCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGTCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGc                          >  1:351303/1‑150 (MQ=255)
                                                                                                                             cttaaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                         >  2:1765362/5‑150 (MQ=38)
                                                                                                                             cttaaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                         >  2:1428509/5‑150 (MQ=38)
                                                                                                                             cATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCATTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                         >  1:415950/1‑150 (MQ=255)
                                                                                                                             cATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                         >  2:1813731/1‑150 (MQ=255)
                                                                                                                              agcAAACAACGCGTACACCGGTGCCTCTCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAg                        <  1:1146147/148‑1 (MQ=255)
                                                                                                                              aTCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGCTACAGCGGTAGATAATCGGGCAg                        >  2:505054/1‑150 (MQ=255)
                                                                                                                              aTCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAg                        >  2:1279588/1‑150 (MQ=255)
                                                                                                                               taaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGCGCAATCCTTCTTTCGCCTGATCGTAAGGCAGGGCGCTTAAATACTCCTGATTCAGCGGTAGATAATCGGGCCGt                       >  2:145290/3‑150 (MQ=255)
                                                                                                                                cAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGtt                      >  1:803209/1‑150 (MQ=255)
                                                                                                                                 aaaCAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTc                     <  2:740742/150‑1 (MQ=255)
                                                                                                                                  aacaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTcc                    >  2:661125/1‑150 (MQ=255)
                                                                                                                                   acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc                   <  1:1022823/150‑1 (MQ=255)
                                                                                                                                      acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc                <  1:152731/150‑1 (MQ=255)
                                                                                                                                         cgTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAAt             <  1:148780/150‑1 (MQ=255)
                                                                                                                                             ccACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAg         <  1:1605978/150‑1 (MQ=255)
                                                                                                                                             ccACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAg         >  1:1761143/1‑150 (MQ=255)
                                                                                                                                              cACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGCCAATCCAGc        >  2:635079/1‑150 (MQ=255)
                                                                                                                                              cACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGc        <  1:365814/150‑1 (MQ=255)
                                                                                                                                              cACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGc        <  1:769038/150‑1 (MQ=255)
                                                                                                                                               aCGGTGCCTCGAGTTGTTTTAACCGTCCGCCAAGCAGGGCGAGGATCGGCGCAACATGGAGAAAACCTTCGTTCGCCGGATCGTCTGGCAGGACGCTGAAAGACTGCAGCTTCAGCGGTAGCATAACGGGCAGTTCCGGACAGTCCAGcg       >  2:263438/1‑149 (MQ=255)
                                                                                                                                                 ggTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCTGGGCGAGGAGCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCTGTTCCCGACACTCCAGCTGc     >  1:73075/1‑150 (MQ=255)
                                                                                                                                                    gCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAgg  >  1:1394655/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGG  >  NC_000913/1536728‑1537025

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