Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,083 A→C 7.6% M217L (ATG→CTG)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,0830AC 7.6% 548.7 / 28.8 157M217L (ATG→CTG) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base A (70/75);  new base C (7/5);  total (77/80)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

AAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCA  >  NC_000913/1788935‑1789230
                                                                                                                                                    |                                                                                                                                                   
aaGTCAGATCCAACGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGAATGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGGCAAAAATGGTCAGa                                                                                                                                                     <  4:319093/149‑1 (MQ=255)
 aGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                                                                                                                    <  1:605883/149‑1 (MQ=255)
 aGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                                                                                                                    >  1:939708/1‑149 (MQ=255)
   tCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGa                                                                                                                                                  >  2:1692273/1‑149 (MQ=255)
    cAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAc                                                                                                                                                 <  1:126125/149‑1 (MQ=255)
     aGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                                                                                                                                >  2:1884541/1‑149 (MQ=255)
         ccACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCt                                                                                                                                            >  4:439792/1‑149 (MQ=255)
          cACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTa                                                                                                                                           <  3:439792/149‑1 (MQ=255)
           aCCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                          >  4:1341151/1‑149 (MQ=255)
            ccGCGAAATGGCTTCGGCACTCGCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATc                                                                                                                                         <  1:1610014/149‑1 (MQ=255)
             cgcgAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                                                                                                                        <  1:384436/149‑1 (MQ=255)
                  aaTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTTTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                                                                                                                   <  4:1632189/149‑1 (MQ=255)
                   aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                                                                                                                  >  2:636733/1‑149 (MQ=255)
                        ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                                                                                                                      >  2:845021/1‑140 (MQ=255)
                          cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCTGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                                                                                                                           >  2:1609178/1‑149 (MQ=255)
                              aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGt                                                                                                                       <  4:864079/149‑1 (MQ=255)
                               ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                      <  1:1630413/149‑1 (MQ=255)
                               ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                      <  1:1609178/149‑1 (MQ=255)
                               ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                      <  3:427390/149‑1 (MQ=255)
                                tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                          >  2:921007/1‑128 (MQ=255)
                                tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                          <  1:921007/128‑1 (MQ=255)
                                     tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                                                                                                                  <  3:1033799/131‑1 (MQ=255)
                                     tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                                                                                                                  >  4:1033799/1‑131 (MQ=255)
                                       tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                     >  4:906259/1‑142 (MQ=255)
                                       tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                     >  2:1271966/1‑142 (MQ=255)
                                       tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                     <  3:906259/142‑1 (MQ=255)
                                       tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                     <  1:1271966/142‑1 (MQ=255)
                                       tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca                                                                                                              >  2:266724/1‑149 (MQ=255)
                                         cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatt                                                                                                            >  4:175438/1‑149 (MQ=255)
                                          ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatta                                                                                                           <  1:1884541/149‑1 (MQ=255)
                                           gTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                                                                                                                        <  2:929972/119‑1 (MQ=255)
                                                ttttAAAAATGGTACCGATGGCAATACGCGGATTGTTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCATCTGACCAttt                                                                                                                                            >  4:752059/1‑108 (MQ=255)
                                                ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                     >  1:140156/1‑101 (MQ=255)
                                                ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                     <  2:140156/101‑1 (MQ=255)
                                                ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc                                                                                                     >  3:1304253/1‑149 (MQ=255)
                                                ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc                                                                                                     >  2:1533372/1‑149 (MQ=255)
                                                              ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaaa                                                                                       >  1:480769/1‑149 (MQ=255)
                                                                      aaTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTa                                                                               >  2:244372/1‑149 (MQ=255)
                                                                       aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTAt                                                                              <  4:1835753/149‑1 (MQ=255)
                                                                       aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTAt                                                                              >  3:861803/1‑149 (MQ=255)
                                                                         aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCa                                                                            <  2:1620581/149‑1 (MQ=255)
                                                                           gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                          >  2:1172899/1‑149 (MQ=255)
                                                                               aTTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCga                                                                      <  1:822872/149‑1 (MQ=255)
                                                                                ttGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgat                                                                     <  4:861803/149‑1 (MQ=255)
                                                                                      tgGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                               <  1:1692273/149‑1 (MQ=255)
                                                                                        gATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGa                                                                                   >  3:1304627/1‑127 (MQ=255)
                                                                                        gATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGa                                                                                   <  4:1304627/127‑1 (MQ=255)
                                                                                           gCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                     <  1:1840729/90‑1 (MQ=255)
                                                                                              aTCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATc                                                                             >  2:428786/1‑127 (MQ=255)
                                                                                              aTCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATc                                                                             <  1:428786/127‑1 (MQ=255)
                                                                                               tCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc                                                                                                     <  1:1064558/102‑1 (MQ=255)
                                                                                               tCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc                                                                                                     >  2:1064558/1‑102 (MQ=255)
                                                                                                ccGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCt                                                     >  2:1380875/1‑149 (MQ=255)
                                                                                                       gCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATAc                                              <  2:179773/149‑1 (MQ=255)
                                                                                                       gCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATAc                                              >  2:392850/1‑149 (MQ=255)
                                                                                                         ccGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGc                                            >  1:1691472/1‑149 (MQ=255)
                                                                                                               cagccaTATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                                                                                                                    <  1:1376710/39‑1 (MQ=255)
                                                                                                               cagccaTATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                                                                                                                    >  2:1376710/1‑39 (MQ=255)
                                                                                                               cagccaTATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACg                                      >  2:453991/1‑149 (MQ=255)
                                                                                                                agccaTATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGa                                     <  4:222673/149‑1 (MQ=255)
                                                                                                                agccaTATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGa                                     <  3:1361666/149‑1 (MQ=255)
                                                                                                                      atGTTCCTCTCGCCAGACAAAAATGGTCAGCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCTCAGCCTGCGATa                                               >  3:692513/1‑133 (MQ=255)
                                                                                                                      atGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATa                                               <  2:800758/133‑1 (MQ=255)
                                                                                                                      atGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATa                                               <  4:692513/133‑1 (MQ=255)
                                                                                                                      atGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGa                               <  3:720764/149‑1 (MQ=255)
                                                                                                                      atGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGa                               <  2:480769/149‑1 (MQ=255)
                                                                                                                         ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca                                                                                                              <  2:329311/67‑1 (MQ=255)
                                                                                                                         ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca                                                                                                              >  1:329311/1‑67 (MQ=255)
                                                                                                                           ccTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          <  3:1885021/149‑1 (MQ=255)
                                                                                                                            ctctCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt                                                                                                   <  2:1481894/75‑1 (MQ=255)
                                                                                                                            ctctCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt                                                                                                   >  1:1481894/1‑75 (MQ=255)
                                                                                                                            ctctCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          >  3:110785/1‑148 (MQ=255)
                                                                                                                            ctctCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          <  2:822393/148‑1 (MQ=255)
                                                                                                                            ctctCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          >  1:822393/1‑148 (MQ=255)
                                                                                                                               tcGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGa                      >  4:1069843/1‑149 (MQ=255)
                                                                                                                                 gCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGa                               >  2:1590767/1‑138 (MQ=255)
                                                                                                                                 gCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGa                               <  1:1590767/138‑1 (MQ=255)
                                                                                                                                 gCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCt                        >  2:384830/1‑145 (MQ=255)
                                                                                                                                   cAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc                                                                                                     <  1:1520601/66‑1 (MQ=255)
                                                                                                                                    aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATc                 >  4:673976/1‑149 (MQ=255)
                                                                                                                                      tgatagaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtg               >  3:79488/7‑149 (MQ=255)
                                                                                                                                      aCAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTGTTATGCGTGGCGGCAAATAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtg               >  3:1587669/1‑149 (MQ=255)
                                                                                                                                       cAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtgg              <  1:392850/149‑1 (MQ=255)
                                                                                                                                        atagaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggt             >  1:722693/5‑149 (MQ=255)
                                                                                                                                          aaaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtgg           >  4:839140/1‑149 (MQ=255)
                                                                                                                                            aTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATcg                                                              >  1:916981/1‑96 (MQ=255)
                                                                                                                                              ggTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGa       <  4:79488/149‑1 (MQ=255)
                                                                                                                                                tCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGttt                                 >  3:1724294/1‑121 (MQ=255)
                                                                                                                                                tCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGttt                                 <  4:592017/121‑1 (MQ=255)
                                                                                                                                                tCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGttt                                 <  4:1724294/121‑1 (MQ=255)
                                                                                                                                                tCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGAtt     >  2:153762/1‑149 (MQ=255)
                                                                                                                                                tCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGAtt     >  4:1371590/1‑149 (MQ=255)
                                                                                                                                                  aGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          >  2:1435277/1‑126 (MQ=255)
                                                                                                                                                  aGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          <  1:1435277/126‑1 (MQ=255)
                                                                                                                                                   gATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt          <  1:589857/141‑1 (MQ=255)
                                                                                                                                                   gATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt          >  2:589857/1‑141 (MQ=255)
                                                                                                                                                   gATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCa  <  1:244372/149‑1 (MQ=255)
                                                                                                                                                   gATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCa  <  1:1172899/149‑1 (MQ=255)
                                                                                                                                                    aTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTGTGATTTAc                          >  4:1249764/1‑124 (MQ=255)
                                                                                                                                                    aTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACACTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAc                          <  3:1249764/124‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCA  >  NC_000913/1788935‑1789230

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