Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,536,877 C→G 6.0% 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 6.0% 544.9 / 19.7 167A145P (GCA→CCA) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base C (77/79);  new base G (3/7);  total (80/87)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 170 total aligned reads displayed.

AAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGC  >  NC_000913/1536728‑1537026
                                                                                                                                                     |                                                                                                                                                     
aaaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGc                                                                                                                                                       <  2:1408064/150‑1 (MQ=255)
aaaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGc                                                                                                                                                       <  2:2110069/150‑1 (MQ=255)
    ttCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcg                                                                                                                                                   <  1:58893/150‑1 (MQ=255)
       aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGAGACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                                                                                <  1:883344/150‑1 (MQ=255)
       aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                                                                                                                                <  2:1627370/150‑1 (MQ=255)
        ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgtt                                                                                                                                               >  2:867662/1‑150 (MQ=255)
          tGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgt                                                                                                                                             >  1:187167/1‑150 (MQ=255)
           gTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                                                                                                                                            >  2:129783/1‑150 (MQ=255)
            ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                                           <  2:64369/150‑1 (MQ=255)
            ttCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                                                                                                                           <  1:828772/150‑1 (MQ=255)
              ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                                                                                                                                         >  2:1894577/1‑150 (MQ=255)
              ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                                                                                                                                         >  1:830724/1‑150 (MQ=255)
              ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttga                                                                                                                                         <  2:2154983/150‑3 (MQ=255)
               ctctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgat                                                                                                                                        <  1:13180/150‑4 (MQ=38)
               ctctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTATTGTCACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGtttgat                                                                                                                                        <  1:1054942/150‑4 (MQ=38)
                   tCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                                                                                                                                    >  2:673634/1‑150 (MQ=255)
                   tCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                                                                                                                                    >  2:757537/1‑150 (MQ=255)
                    cccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAc                                                                                                                                   >  2:1820437/1‑150 (MQ=255)
                     ccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAcc                                                                                                                                  >  1:1242909/1‑150 (MQ=255)
                     ccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGGTTTACGCGAcc                                                                                                                                  >  1:1463107/1‑150 (MQ=255)
                          cTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCaa                                                                                                                             >  1:1856606/1‑150 (MQ=255)
                           tGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                                                                                                                            <  1:2115663/150‑1 (MQ=255)
                           tGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTCCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                                                                                                                            >  1:428677/1‑150 (MQ=255)
                              aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAACACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                         >  1:298255/1‑150 (MQ=255)
                              aTCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTCTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg                                                                                                                         <  1:170849/150‑1 (MQ=255)
                               tCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg                                                                                                                        <  2:1529276/150‑1 (MQ=255)
                               tCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg                                                                                                                        <  1:1736245/150‑1 (MQ=255)
                                cAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAggg                                                                                                                       >  2:1020722/1‑150 (MQ=255)
                                cAAGCGCCTGGCGGGTGTCAGCGCGCTCTTCGCTGTTTCCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAggg                                                                                                                       <  1:1413529/150‑1 (MQ=255)
                                    cgcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAg                                                                                                                   >  1:1031404/1‑150 (MQ=255)
                                     gcCTGGCGTGTGACATCGCGATCTTCGCTTTTTAATTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCATCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                                  <  2:377078/150‑1 (MQ=255)
                                     gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAACCAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                                  <  1:1460625/150‑1 (MQ=255)
                                     gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                                  <  1:1646900/150‑1 (MQ=255)
                                     gcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAgg                                                                                                                  >  2:901143/1‑150 (MQ=255)
                                            gggTGGCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGgcg                                                                                                           <  1:1309764/150‑1 (MQ=255)
                                              gtgtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCa                                                                                                         <  1:983841/150‑1 (MQ=255)
                                               tgtCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCaa                                                                                                        <  1:1383397/150‑1 (MQ=255)
                                                gtCATCGCGCTCTTCGCTGATTACTTGTTTTGTGACACTGTCACTTGACAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAAc                                                                                                       <  1:1771486/150‑1 (MQ=255)
                                                gtCATCGCGCGCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAAc                                                                                                       <  1:1787053/150‑1 (MQ=255)
                                                    tCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTg                                                                                                   >  2:1631713/1‑150 (MQ=255)
                                                    tCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTACCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCACAACATTg                                                                                                   <  2:428677/150‑1 (MQ=255)
                                                      gcgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAg                                                                                                 >  2:1697542/1‑150 (MQ=255)
                                                       cgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                                                                                <  2:2004869/150‑1 (MQ=255)
                                                       cgcTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGc                                                                                                >  2:425889/1‑150 (MQ=255)
                                                         ctctTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCaa                                                                                              >  2:87568/1‑150 (MQ=255)
                                                              cGCTGTATACTTGTTTTGTGACCCTGTCACTTGAAAGGGCGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCCTACCACGGTGCCTCGCGTTGTTTTATGCGACCGCCCAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCctt                                                                                         >  2:1071623/1‑150 (MQ=255)
                                                               gCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttc                                                                                        <  2:269585/150‑1 (MQ=255)
                                                                 tGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCACGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCcttctt                                                                                      <  1:1478103/150‑1 (MQ=255)
                                                                   tttACTTGGTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTc                                                                                    <  2:144869/150‑1 (MQ=255)
                                                                     tACTTGTTTTGTGACACTGTCACTCGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGCTCGGCGCAACATTGCGCAATCCGTCTTTCGc                                                                                  >  2:631160/1‑150 (MQ=255)
                                                                      aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                                 <  1:1131566/150‑1 (MQ=255)
                                                                      aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                                 <  2:1286907/150‑1 (MQ=255)
                                                                      aCTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGAGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGcc                                                                                 >  1:1791751/1‑150 (MQ=255)
                                                                         tGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGa                                                                              >  2:608738/1‑150 (MQ=255)
                                                                         tGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGa                                                                              <  1:706456/150‑1 (MQ=255)
                                                                             ttGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGt                                                                          <  1:1358683/150‑1 (MQ=255)
                                                                                 gACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGTCTGATCGTCAgg                                                                      >  2:1524979/1‑150 (MQ=255)
                                                                                 gACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAgg                                                                      >  1:2203809/1‑150 (MQ=255)
                                                                                   cacTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCa                                                                    <  1:1100849/150‑1 (MQ=255)
                                                                                    acTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAg                                                                   <  1:2143641/150‑1 (MQ=255)
                                                                                     cTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCAAAGCAGGGCGAGGATCGGAGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGt                                                                  >  1:1611780/1‑150 (MQ=255)
                                                                                      tGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCCGGCAGTa                                                                 >  2:2002742/1‑150 (MQ=255)
                                                                                      tGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGACGATCAGGcgaaa                                                                 <  1:373263/150‑5 (MQ=255)
                                                                                          aCTTGAAAGGGAGCTTCCCGCCAGTTGCGGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCt                                                             <  2:246725/150‑1 (MQ=255)
                                                                                            ttGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTa                                                           <  1:485967/150‑1 (MQ=255)
                                                                                              gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                         <  1:2024438/150‑1 (MQ=255)
                                                                                              gAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTaaa                                                         >  2:1928230/1‑150 (MQ=255)
                                                                                               aaaGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAAt                                                        >  2:119120/1‑150 (MQ=255)
                                                                                                  gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                                     >  2:52353/1‑150 (MQ=255)
                                                                                                  gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                                     >  2:505398/1‑150 (MQ=255)
                                                                                                  gggAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACt                                                     <  2:57156/150‑1 (MQ=255)
                                                                                                      gCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCTGATCGTCAGGCAGTACGCTTAAATACTCCAGa                                                <  2:1031404/150‑1 (MQ=255)
                                                                                                      gCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCCTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGCCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAg                                                 >  2:1925738/1‑150 (MQ=255)
                                                                                                       cTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGCCCGCCAAGAAGTGCGAGGATCGGCGCATCATTGAGCAATACTTCTTTCGCCAGATCGTCAGGCACTACGCTTAAATACTCCAGa                                                >  2:1217111/1‑150 (MQ=255)
                                                                                                        ttCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGAt                                               >  1:508438/1‑150 (MQ=255)
                                                                                                         tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATa                                              >  2:266839/1‑150 (MQ=255)
                                                                                                         tCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATa                                              <  2:1727786/150‑1 (MQ=255)
                                                                                                            cGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAg                                           >  1:288980/1‑150 (MQ=255)
                                                                                                             gCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATAGAGCAATCCTTCATTCGCCAGATCGTCAGCCACTACGCTTAAATAATCCAGATACAGc                                          >  1:263623/1‑150 (MQ=255)
                                                                                                               cAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCgg                                        <  1:688896/150‑1 (MQ=255)
                                                                                                                    gcagcaGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                                   <  1:867662/150‑1 (MQ=255)
                                                                                                                    gcagcaGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                                   <  2:1084290/150‑1 (MQ=255)
                                                                                                                    gcagcaGAGCATCAAACAACGCGGACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGAt                                   <  1:1639230/150‑1 (MQ=255)
                                                                                                                     cagcagAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCATGGCGAGGATAGGCGCAACATTGAGCAATCCTTCTGTCGCCTGATCGTCAGGCAGTACGCTTAAATACTACAGATACAGCGCTAGATa                                  >  2:1108277/1‑150 (MQ=255)
                                                                                                                       gcagAGCATCAAACAACGAGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACTCTTAAATACTCCAGATACAGCGGTAGATAAt                                >  2:2120557/1‑150 (MQ=255)
                                                                                                                            gCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGc                           <  1:1820437/150‑1 (MQ=255)
                                                                                                                             cttaaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                          >  1:747002/5‑150 (MQ=38)
                                                                                                                             cttaaaaCAACGCGAGGCACCGTGGTACGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCa                          >  1:1907705/5‑150 (MQ=38)
                                                                                                                                 aaaCAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTc                      >  1:1014066/1‑150 (MQ=255)
                                                                                                                                  aacaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGAAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTcc                     >  2:1396611/1‑150 (MQ=255)
                                                                                                                                   acaacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTccc                    <  2:953832/150‑1 (MQ=255)
                                                                                                                                     aacGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGa                  <  2:1687695/150‑1 (MQ=255)
                                                                                                                                      acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc                 <  2:1966225/150‑1 (MQ=255)
                                                                                                                                      acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc                 <  1:169827/150‑1 (MQ=255)
                                                                                                                                      acGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGAc                 <  2:1639616/150‑1 (MQ=255)
                                                                                                                                       cgcgTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACa                >  1:1256585/1‑150 (MQ=255)
                                                                                                                                         cgTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGTGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAAt              >  2:2202217/1‑150 (MQ=255)
                                                                                                                                              cACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGc         <  1:443135/150‑1 (MQ=255)
                                                                                                                                               aCGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCt        <  1:1279255/150‑1 (MQ=255)
                                                                                                                                               aCGGTGCCTCGAGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCTACATTGAGCTATCCTTCTTTCGCCTGTTCGTCAGGCAGTCCGCTTAAATACTCCAGATTGCGCGGTAGATTAGCGGGCAGTTCCCGACTCTCCAgat        >  1:381335/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGGCGCAACATTGAGCAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGC  >  NC_000913/1536728‑1537026

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: