Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,272,542 T→A 5.4% C12* (TGT→TGA chaC → cation transport regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,272,5420TA 5.4% 265.4 / 2.9 89C12* (TGT→TGAchaCcation transport regulator
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base A (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGG  >  NC_000913/1272393‑1272684
                                                                                                                                                     |                                                                                                                                              
aaGTGGTCCGTGGATTGCATATTGTCCCGATAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTCGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGt                                                                                                                                                >  1:496625/1‑150 (MQ=255)
     gTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCCGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTTATGAATGCCGATTGTAAAAc                                                                                                                                           <  2:296994/150‑1 (MQ=255)
      tCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCTGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACg                                                                                                                                          >  1:880213/1‑150 (MQ=255)
       ccGTTGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACgg                                                                                                                                         <  1:573404/150‑1 (MQ=255)
       ccGTGGATTGCTTATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTACGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGTTGATAAATAGTGATATATTGATGAATGCCGATTGTAAAAAgg                                                                                                                                         >  1:142828/1‑150 (MQ=255)
       ccGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACgg                                                                                                                                         >  1:1762240/1‑150 (MQ=255)
        cGTGGATTGCATAATGTCCCGTTAGTGGGTTCAAAATGAGCAGTAAAAATGTCCGGAAGCCACCAAACAGGTGTCGCAGGGTAGTATGCAGTGGCGGAAGTGTAAGGGGATAACGCGGGATTTCTTGATGAATGCCGATTGTAAAACGGc                                                                                                                                        <  2:871410/150‑1 (MQ=255)
         gTGGATTGCATATTGTCCCGTTCGTGGTTTCAATATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCa                                                                                                                                       <  2:1314405/150‑1 (MQ=255)
          tGGATTGCATATTGTCCCTTTAGTGGTGTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAACAGTTGTCTCAGGGAAGTATGCAGTTGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGGTTTACAATCGGCAt                                                                                                                                      <  1:335584/150‑1 (MQ=255)
          tGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCAt                                                                                                                                      <  2:1807532/150‑1 (MQ=255)
                gCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGATGTATGCAGTGGCGGAAGTGTAAGGTGATAAAGCGTGATTTCTTGATGAATGCAGATTGTAAAACGGCATTTGGTg                                                                                                                                >  1:1188332/1‑150 (MQ=255)
                gCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGTTTTACAATCGGCAttcatca                                                                                                                                <  2:529290/150‑6 (MQ=255)
                 cATATTGTCCCGTTAGTGGTTGCAAAAAGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGc                                                                                                                               <  1:1647105/150‑1 (MQ=255)
                  atatTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGcc                                                                                                                              <  2:1717627/150‑1 (MQ=255)
                     ttGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCAtt                                                                                                                           >  2:322532/1‑150 (MQ=255)
                         cccGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTgaag                                                                                                                       >  1:1586501/1‑150 (MQ=255)
                             ttAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATc                                                                                                                   <  1:1270287/150‑1 (MQ=255)
                             gtAGTGGTTTAAAAATTAGCATTAAAAATGTCCGGAAGACAACAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAATGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATc                                                                                                                   <  2:283227/149‑1 (MQ=255)
                               aGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAc                                                                                                                 <  2:594065/150‑1 (MQ=255)
                                  ggTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCActct                                                                                                              >  2:1413733/1‑150 (MQ=255)
                                   gTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCTGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCAGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCtt                                                                                                             >  1:85200/1‑150 (MQ=255)
                                     ttCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTAt                                                                                                           <  1:1718381/150‑1 (MQ=255)
                                      tCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATg                                                                                                          >  1:1314259/1‑150 (MQ=255)
                                       cAAAATGAGCAGTAAAAATGTCCGGAAGACACCAATAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTTAGGTGAAAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTTCCATTGAAGAATCACTCTTATgg                                                                                                         >  2:283472/1‑150 (MQ=255)
                                       cAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATgg                                                                                                         <  1:669366/150‑1 (MQ=255)
                                        aaaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGt                                                                                                        >  2:623455/1‑150 (MQ=255)
                                         aaaTGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTc                                                                                                       <  1:561038/150‑1 (MQ=255)
                                               gCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGAATGCAGTGGCGGAAGTGTTAGGTGTTAACGCGTGATTACTTGATGAATGCCGATTGTAAAACGGCATTATGTGCCATTGAAGAATCACTCTTATGGACAGCAGa                                                                                                 >  1:1103159/1‑150 (MQ=255)
                                                cAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaa                                                                                                >  2:1086317/1‑150 (MQ=255)
                                                 aGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaac                                                                                               <  2:496625/150‑1 (MQ=255)
                                                 aGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGaac                                                                                               <  1:525985/150‑1 (MQ=255)
                                                      aaaTGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACgg                                                                                          >  2:354886/1‑150 (MQ=255)
                                                        aTGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggc                                                                                        <  1:757922/150‑1 (MQ=255)
                                                        aTGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGAAGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGggc                                                                                        <  2:132201/150‑1 (MQ=255)
                                                            ccGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCt                                                                                    >  2:1436221/1‑150 (MQ=255)
                                                             cGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCtt                                                                                   >  2:1439840/1‑150 (MQ=255)
                                                              ggAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTc                                                                                  <  1:517165/150‑1 (MQ=255)
                                                              ggAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTc                                                                                  >  2:669728/1‑150 (MQ=255)
                                                                aaGACACCACAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGc                                                                                <  1:412241/150‑1 (MQ=255)
                                                                    cacCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggc                                                                            >  2:1222943/1‑150 (MQ=255)
                                                                      ccAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcgg                                                                          >  1:1540568/1‑150 (MQ=255)
                                                                      ccAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcgg                                                                          >  1:1546014/1‑150 (MQ=255)
                                                                       cAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTggcggc                                                                         >  1:1394809/1‑150 (MQ=255)
                                                                        aaaaaGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGgcggcg                                                                        >  2:1223887/1‑150 (MQ=255)
                                                                         aaaaGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGa                                                                       >  1:1415762/1‑150 (MQ=255)
                                                                               tGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTgg                                                                 <  1:1818122/150‑1 (MQ=255)
                                                                                gTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCAATCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGc                                                                >  1:1291012/1‑150 (MQ=255)
                                                                                  cGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATAGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCCGCTTCGCTGTCGGCGACGCTGTCtt                                                              >  1:286412/1‑150 (MQ=255)
                                                                                   gCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGAGCCATTGAAGAATAACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTg                                                             >  1:994721/1‑150 (MQ=255)
                                                                                         aaGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGAcc                                                       <  1:452274/150‑1 (MQ=255)
                                                                                         aaGGATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAACACGGCATTTGGTGCCAGTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCGGGCTTGTCGAcc                                                       <  1:288376/150‑1 (MQ=255)
                                                                                             aTGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATTGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCtgat                                                   <  1:1791152/150‑1 (MQ=255)
                                                                                             aTGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCtgat                                                   <  2:1213728/150‑1 (MQ=255)
                                                                                                cAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGc                                                                <  2:516917/134‑1 (MQ=255)
                                                                                                cAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGc                                                                >  1:516917/1‑134 (MQ=255)
                                                                                                cAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGCCCTGATGAg                                                >  2:296303/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                   <  1:1086317/150‑1 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                   <  1:1554368/150‑1 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTCGa                                   >  2:547258/1‑150 (MQ=255)
                                                                                                             gtAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTAGGTGCCATTGAAGACTCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGa                                   <  1:280879/150‑1 (MQ=255)
                                                                                                               aaGCTGATAACGCGTGATTTCTTGATGAATGCCGCTTGAAAAACGGCATTGGGTGCCAGTGCAGAAGCACTCTTAGGGTCAGCAGATCAACGGGCGGCGTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATc                                 <  1:101066/150‑1 (MQ=255)
                                                                                                                  gTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTAGACCTGATGAGGGACCGGTGTGGATCTTc                              >  2:76652/1‑150 (MQ=255)
                                                                                                                   tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCg                             <  2:436796/150‑1 (MQ=255)
                                                                                                                   tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTCGATCTACg                             >  2:823226/1‑150 (MQ=255)
                                                                                                                   tGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCAGTGAAGAATCACTCTTATGGGCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCg                             <  1:740076/150‑1 (MQ=255)
                                                                                                                   ggATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCg                             <  2:1388957/149‑1 (MQ=255)
                                                                                                                    gATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTCGATCTTCgt                            >  2:1735617/1‑149 (MQ=255)
                                                                                                                     aTAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGc                           <  2:1613245/150‑1 (MQ=255)
                                                                                                                      tAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCt                          <  2:797779/150‑1 (MQ=255)
                                                                                                                         cgcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                       <  1:354886/150‑1 (MQ=255)
                                                                                                                         cgcgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATg                       >  1:685151/1‑150 (MQ=255)
                                                                                                                           cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGa                     >  1:1074462/1‑150 (MQ=255)
                                                                                                                           cgTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGa                     >  1:364732/1‑150 (MQ=255)
                                                                                                                               aTTTCTTGATGAATGCCGTTTGTACAACGGCATTTGGTGCAATTGAAGAAGCACTCTTATGGGCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGt                 <  1:468869/150‑1 (MQ=255)
                                                                                                                                 ttCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTg               <  2:652439/150‑1 (MQ=255)
                                                                                                                                  tCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGa              <  2:538571/150‑1 (MQ=255)
                                                                                                                                   cTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGCTGTGGATCATCGGCTATGGATCGTTGAt             >  1:308002/1‑150 (MQ=255)
                                                                                                                                    ttGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtg            >  2:1276355/1‑150 (MQ=255)
                                                                                                                                     tgatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt           <  2:1447557/150‑1 (MQ=255)
                                                                                                                                     tgatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt           >  2:540258/1‑150 (MQ=255)
                                                                                                                                     caccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgt           >  1:1700284/6‑150 (MQ=255)
                                                                                                                                      gatgaATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGGTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGAtgtg          <  2:729296/150‑1 (MQ=255)
                                                                                                                                       ccaaaTGCCGTTTTACAATCGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTgg         >  1:1641449/4‑150 (MQ=255)
                                                                                                                                        tgaATGCCGATTGTAACACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGTGCGGCTTCGCTGGCGGAGACGCTGGCTTGTCTACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGa        >  1:277418/1‑150 (MQ=255)
                                                                                                                                        tgaATGCCGATTGTAAAACGGCATTTGGTGTCATTGAAGAATAACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGAAAGTTGAttagga        >  1:734490/1‑145 (MQ=255)
                                                                                                                                          aaTGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAAt      >  1:497057/1‑150 (MQ=255)
                                                                                                                                           aTGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATc     >  2:358431/1‑150 (MQ=255)
                                                                                                                                             gCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCg   >  1:750332/1‑150 (MQ=255)
                                                                                                                                             gCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTACGGCAATGGATCGTTGATGTGGAATCCg   >  2:1795379/1‑150 (MQ=255)
                                                                                                                                              ccGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCgg  <  2:60340/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AAGTGGTCCGTGGATTGCATATTGTCCCGTTAGTGGTTTCAAAATGAGCAGTAAAAATGTCCGGAAGACACCAAAAAGTTGTCGCAGGGAAGTATGCAGTGGCGGAAGTGTAAGGTGATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGG  >  NC_000913/1272393‑1272684

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