Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,361 C→T 7.1% H231Y (CAC→TAC)  guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,872,3610CT 7.1% 331.8 / 10.6 112H231Y (CAC→TAC) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base C (60/44);  new base T (6/2);  total (66/46)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 113 total aligned reads displayed.

TATGCCGAAAGCAAGGCACTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGAAGA  >  NC_002947/4872214‑4872509
                                                                                                                                                   |                                                                                                                                                    
tATGCCGAAAGCAAGGCACTGATCGAGCGCTGGCATCGCAAGGGACGCCTGCACTCAGCAGCTCCCGCTCGCCTCGCTGCAACCCGCACAGCTGAACGAGATGCCCTTGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGcac                                                                                                                                                    <  2:1477649/150‑1 (MQ=255)
      gAAAGCATGGCACTGATCGAGCGCTGGCATGGCAAGGGCCGTCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAGCCAGCACACCGGAACAATTGACCCTGGTCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCAc                                                                                                                                              <  2:771796/150‑1 (MQ=255)
      gAAAGCAAGGCACTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCAc                                                                                                                                              >  2:1357890/1‑150 (MQ=255)
       aaaGCAAGGCACTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCAcc                                                                                                                                             <  1:1848015/150‑1 (MQ=255)
             aGGCACTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCgg                                                                                                                                       >  1:1791447/1‑150 (MQ=255)
                 aCTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGaaaa                                                                                                                                   >  2:164126/1‑150 (MQ=255)
                  cTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAAc                                                                                                                                  <  1:302751/150‑1 (MQ=255)
                              tggcatggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATc                                                                                                                      >  1:1754113/1‑150 (MQ=255)
                                gcatggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTGCATGTACACGccggg                                                                                                                                            >  2:1801802/1‑122 (MQ=255)
                                gcatggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTGCATGTACACGccggg                                                                                                                                            <  1:1801802/126‑5 (MQ=255)
                                  atggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACt                                                                                                                  >  2:751566/1‑150 (MQ=255)
                                  atggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACt                                                                                                                  >  2:89775/1‑150 (MQ=255)
                                  atggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACt                                                                                                                  >  2:2222311/1‑150 (MQ=255)
                                    ggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTgg                                                                                                                <  1:2059582/150‑1 (MQ=255)
                                    ggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTgg                                                                                                                <  1:1357890/150‑1 (MQ=255)
                                    ggcaAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTgg                                                                                                                <  1:164126/150‑1 (MQ=255)
                                           gccgccTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGt                                                                                                         <  1:751566/150‑1 (MQ=255)
                                            ccgccTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTc                                                                                                        <  1:1303329/150‑1 (MQ=255)
                                             cgccTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGCCCCTGGCCTGCCAATTGCCCAAGGAACACCCCGGCGTGTACCTGCCCACGCCCCTGTCGGACAACCTCAAGGAAATCGACTGCGTCAAGTcc                                                                                                       >  2:117145/1‑149 (MQ=255)
                                                cTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTg                                                                                                    >  1:2258833/1‑150 (MQ=255)
                                                 tGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCaa                                                                                                                             >  2:1505764/1‑124 (MQ=255)
                                                   cACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTc                                                                                                 <  2:1754113/150‑1 (MQ=255)
                                                       aTGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCg                                                                                             >  2:1408086/1‑150 (MQ=255)
                                                        tGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGa                                                                                            <  1:1542994/150‑1 (MQ=255)
                                                               aCCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCg                                                                                                       <  2:1747281/132‑1 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGt                                                                                                                                          >  2:1649604/1‑95 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGt                                                                                                                                          <  1:1649604/95‑1 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAAc                                                                                                                                  >  1:2179009/1‑103 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAAc                                                                                                                                  <  2:2179009/103‑1 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCt                                                                                                     >  1:289276/1‑132 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCt                                                                                                     <  2:289276/132‑1 (MQ=255)
                                                                 cccGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAggg                                                                                   >  1:429974/1‑150 (MQ=255)
                                                                   cgcgCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTc                                                                                                                                         >  2:418951/1‑94 (MQ=255)
                                                                   cgcgCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGtt                                                                                 >  2:1568223/1‑150 (MQ=255)
                                                                        ttCGCCCCAAACAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTTCATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAACTCGCTGTTCCCCGAGCAGAAGTGTGATCTg                                                                            >  2:202479/1‑150 (MQ=255)
                                                                         tCGCCCCAACCAGCACCCCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGa                                                                                       <  1:922613/138‑1 (MQ=255)
                                                                         tCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAACCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGa                                                                                       >  2:922613/1‑138 (MQ=255)
                                                                            ccccAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTg                                                                                                                 <  2:1586220/109‑1 (MQ=255)
                                                                            ccccAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTg                                                                                                                 >  1:1586220/1‑109 (MQ=255)
                                                                             cccAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAg                                                                                                                            >  2:754345/1‑97 (MQ=255)
                                                                              ccAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTc                                                                      >  1:1452013/1‑150 (MQ=255)
                                                                                     gCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGCAGGGTTATCTGGACGTCTACGATc                                                               >  2:1207464/1‑150 (MQ=255)
                                                                                       acacCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCATGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCAc                                                             >  2:1150613/1‑150 (MQ=255)
                                                                                         acCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGt                                                                                                              >  2:1204304/1‑99 (MQ=255)
                                                                                         acCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGt                                                                                                              <  1:1204304/99‑1 (MQ=255)
                                                                                         acCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACtt                                                           >  1:2094050/1‑150 (MQ=255)
                                                                                           cGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCg                                                         >  2:897400/1‑150 (MQ=255)
                                                                                             gAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAg                                                       <  1:89775/150‑1 (MQ=255)
                                                                                              aaCAATTGACCCTGGCCGGCCAATTGCTCAAGGAAGACCCCGGCGCGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGACGGGTTATCTGGACGTCTCCGATCACTGCGAgc                                                      >  1:112894/1‑150 (MQ=255)
                                                                                                 aaTTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGa                                                                                            <  1:709312/109‑1 (MQ=255)
                                                                                                 aaTTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGa                                                                                            >  2:709312/1‑109 (MQ=255)
                                                                                                 aaTTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAgctgc                                                   >  1:877281/1‑150 (MQ=255)
                                                                                                     gACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAg                                                                                        >  1:1575/1‑109 (MQ=255)
                                                                                                     gACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAg                                                                                        <  2:1575/109‑1 (MQ=255)
                                                                                                        ccTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTc                                                                                                             <  1:1677271/85‑1 (MQ=255)
                                                                                                        ccTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTc                                                                                                             >  2:1677271/1‑85 (MQ=255)
                                                                                                               ggccAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCg                                     >  2:1012225/1‑150 (MQ=255)
                                                                                                                 ccAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGgt                                   <  2:1152116/150‑1 (MQ=255)
                                                                                                                    aTTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAgg                                                                                    <  2:254887/98‑1 (MQ=255)
                                                                                                                    aTTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAgg                                                                                    >  1:254887/1‑98 (MQ=255)
                                                                                                                     ttGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGttt                               <  1:897400/150‑1 (MQ=255)
                                                                                                                     ttGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGCACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGttt                               <  1:646613/150‑1 (MQ=255)
                                                                                                                         tCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGccc                           >  1:1721695/1‑150 (MQ=255)
                                                                                                                          cAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCa                          >  1:934563/1‑150 (MQ=255)
                                                                                                                           aaGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACGGCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCAc                         <  1:870789/150‑1 (MQ=255)
                                                                                                                           aaGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGAc                                                                         >  1:366352/1‑102 (MQ=255)
                                                                                                                           aaGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGAc                                                                         <  2:366352/102‑1 (MQ=255)
                                                                                                                           aaGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCAc                         >  1:363941/1‑150 (MQ=255)
                                                                                                                           aaGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCACCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCAc                         >  1:770366/1‑150 (MQ=255)
                                                                                                                            aGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAACACTAAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACg                        >  1:1253487/1‑150 (MQ=255)
                                                                                                                            aGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACg                        >  1:198277/1‑150 (MQ=255)
                                                                                                                             ggAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACgg                       <  2:1452013/150‑1 (MQ=255)
                                                                                                                             ggAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCCACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACgg                       >  1:323545/1‑150 (MQ=255)
                                                                                                                              gAACCCCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGc                      <  2:112894/150‑1 (MQ=255)
                                                                                                                              gAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGc                                                                                                      >  1:349522/1‑70 (MQ=255)
                                                                                                                              gAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGc                                                                                                      <  2:349522/70‑1 (MQ=255)
                                                                                                                                 gacaggtgCGTGTGCATGTACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTa                   >  1:2057479/8‑150 (MQ=255)
                                                                                                                                  acCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGaa                                                                                      >  2:568019/1‑82 (MQ=255)
                                                                                                                                  acCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGaa                                                                                      <  1:568019/82‑1 (MQ=255)
                                                                                                                                    cccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTc                                                          <  2:600558/108‑1 (MQ=255)
                                                                                                                                    cccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTc                                                          >  1:600558/1‑108 (MQ=255)
                                                                                                                                    cccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTacac                >  1:1474824/1‑150 (MQ=255)
                                                                                                                                    aggtgCGTGTGCATGTACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTacac                <  2:2057479/146‑1 (MQ=255)
                                                                                                                                     ccGGCGTGTACATGCACCCGCCCCTGTCGGAAAACCTCAAGGAAATCGACTGGGGCAAGTCGCTGTTCCCCGAGCAGAAGGGGTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGATCGGTGCTTGCCCCCGGCGTACAcc               >  2:1106522/1‑150 (MQ=255)
                                                                                                                                     ccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAgctgct                                                  <  1:493618/115‑1 (MQ=255)
                                                                                                                                     ccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAgctgct                                                  >  2:493618/1‑115 (MQ=255)
                                                                                                                                     ccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAg                                                       <  1:1905641/110‑1 (MQ=255)
                                                                                                                                     ccGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAg                                                       >  2:1905641/1‑110 (MQ=255)
                                                                                                                                      gtgCGTGTGCATGTACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGGTATCTGTACGCCGACTATCACTTCGAGCTGCTGGGCAAGCGCTCGGGGTGGGCCCACGGCGTACACCt              >  2:1706166/3‑150 (MQ=255)
                                                                                                                                      cGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCt              >  1:1585049/1‑150 (MQ=255)
                                                                                                                                      cGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCt              >  1:17909/1‑150 (MQ=255)
                                                                                                                                          gtgtACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGc          >  1:307635/1‑150 (MQ=255)
                                                                                                                                              gcATGTACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGcgacg      >  2:2112064/2‑150 (MQ=255)
                                                                                                                                                aTGTACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGaa    >  1:2135746/1‑150 (MQ=255)
                                                                                                                                                aTGTACACGCACCTGTCGGAAAACATCACGTAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGAGGaa    >  1:64032/1‑150 (MQ=255)
                                                                                                                                                aTGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGaa    >  2:1627883/1‑150 (MQ=255)
                                                                                                                                                 tGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGGCGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGttt                               >  1:1765675/1‑122 (MQ=255)
                                                                                                                                                 tGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGAag   >  2:2058444/1‑150 (MQ=255)
                                                                                                                                                  gCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGAaga  <  1:1264243/150‑1 (MQ=255)
                                                                                                                                                  gCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCGGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGAaga  <  1:202479/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATGCCGAAAGCAAGGCACTGATCGAGCGCTGGCATGGCAAGGGCCGCCTGCACTATGCAGTCACCCCGCGCTTCGCCCCAACCAGCACACCGGAACAATTGACCCTGGCCGGCCAATTGCTCAAGGAACACCCCGGCGTGTACATGCACACGCACCTGTCGGAAAACCTCAAGGAAATCGACTGGGTCAAGTCGCTGTTCCCCGAGCAGAAGGGTTATCTGGACGTCTACGATCACTTCGAGCTGCTGGGCGAGCGCTCGGTGTTTGCCCACGGCGTACACCTGTGCGACGAAGA  >  NC_002947/4872214‑4872509

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