Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,129 C→T 5.6% noncoding (67/84 nt) agrA → inactive antisense sRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,1290CT 5.6% 388.5 / 8.1 127noncoding (67/84 nt)agrAinactive antisense sRNA
Reads supporting (aligned to +/- strand):  ref base C (47/70);  new base T (4/3);  total (53/74)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 132 total aligned reads displayed.

AACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647980‑3648278
                                                                                                                                                      |                                                                                                                                                     
aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc                                                                                                                                                       >  2:1009628/1‑150 (MQ=255)
  cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaa                                                                                                                                                     >  2:408007/1‑150 (MQ=255)
     gATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt                                                                                                                                                  >  1:148183/1‑150 (MQ=255)
        aaTCATACAGTACATGCAGGTTATAAAACCAGCAAGTCCTTGCAATAGTTTCAGTTTGGTTTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCTGAGCCACCAACTCAGGGCTTGAAAGTAAAAAACCGACGCAAAGTCgg                                                                                                                                               <  1:1765188/150‑1 (MQ=38)
         aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                              >  1:725518/1‑150 (MQ=255)
            aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                                                                           <  1:948672/150‑1 (MQ=255)
              aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                         <  1:1017502/150‑1 (MQ=255)
              aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                         <  2:1660951/150‑1 (MQ=255)
               cAGTACATGCAGTTTATAAAACCAGCACGTCCTTGCAAAAGTTTCAGTATGTTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTGAGAGCCACCAACTCAGTGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  2:249278/150‑1 (MQ=37)
                aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGAAAAAAACCGACGCAAAGTCGGtttttatt                                                                                                                                       >  2:1878783/1‑147 (MQ=255)
                 gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                       <  2:990900/150‑3 (MQ=38)
                    cgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAAGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATc                                                                                                                                   >  2:1937197/3‑150 (MQ=37)
                     aTGCAGTTTATAAAACCAGCACGTCCTTGCAAAAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttactttc                                                                                                                                   <  1:2115433/150‑7 (MQ=37)
                      tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCg                                                                                                                                 >  1:477254/1‑150 (MQ=255)
                      tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCg                                                                                                                                 >  2:1266021/1‑150 (MQ=255)
                        cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGa                                                                                                                               >  2:1709258/1‑150 (MQ=255)
                        cAGGTTATAAAACCAGCACCTCCTTGCAATTCTATCAGTAAGGAATTAGCATTGAAGCGTAAGATTATGGCTAAAACAATCCAGACAAATCCCCCAACACAGCGATGGAAAGAAAAAAACCGACGAAAAGTCCATTTTTATAcgacctca                                                                                                                               >  2:422972/1‑143 (MQ=255)
                         aGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGAt                                                                                                                              <  1:994643/150‑1 (MQ=255)
                         aGGTTATAAAACCAGCACGTCCTTGCAAAAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCACAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGAt                                                                                                                              <  1:2002170/150‑1 (MQ=255)
                          ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACTAAAAGTCGTTTTTTTTACATCCGTAta                                                                                                                             >  2:341832/1‑149 (MQ=37)
                            ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCg                                                                                                                           <  1:1958687/150‑1 (MQ=255)
                              ataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGa                                                                                                                         <  1:156802/150‑1 (MQ=255)
                               taAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTAGATCCGGagtgggcc                                                                                                                        >  1:2075551/1‑143 (MQ=255)
                                aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACa                                                                                                                       >  2:2038618/1‑150 (MQ=255)
                                aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACa                                                                                                                       >  2:389371/1‑150 (MQ=255)
                                 aaaCTAGCAACTTCTTGCAATATTTTCCTTATGGTATTAGAATTGATGTGTTAGATTAGTTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcct                                                                                                                      <  2:480722/150‑3 (MQ=11)
                                     cAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCt                                                                                                                  <  2:86759/150‑1 (MQ=255)
                                          cGTCATTGCAATAGTATCAGTATGGGATTAGCATTGATGCGGTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCACCTCAGGGCTGGTAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                                                             <  2:326083/150‑1 (MQ=25)
                                           gTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAatat                                                                                                            <  2:470653/150‑1 (MQ=255)
                                             ccTTGCAATAGTTTCATTATGGTATTAGCATTGTTTCGTTAGATGATGGCTATCTCAATCCAGTCCGAGCCAACAACTCAGGGATTGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                          <  1:214649/150‑1 (MQ=37)
                                             ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGATCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                          <  1:1949781/150‑1 (MQ=255)
                                             ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                          >  2:890261/1‑150 (MQ=255)
                                              cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACACCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                         <  1:779244/150‑1 (MQ=255)
                                               ttGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                                                                        <  1:1397811/150‑1 (MQ=255)
                                                  cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACATCCATTAATATGACGaag                                                                                                     >  1:932870/1‑148 (MQ=18)
                                                  cAATAGTTTCAGTATGGTATTAGCATTGATGAGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGTCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATACGGAAAAGCTTTAAAATGACGATg                                                                                                     >  1:1077628/1‑150 (MQ=25)
                                                   aaTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                    <  1:408007/150‑1 (MQ=255)
                                                   aaTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                    <  2:822714/150‑1 (MQ=255)
                                                     tAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                                                                  <  1:1491912/150‑1 (MQ=255)
                                                             gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                                                                          >  1:1140954/1‑150 (MQ=255)
                                                              tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAACAAGTGaaat                                                                                         >  2:193249/1‑149 (MQ=255)
                                                                  gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                                                                     <  2:954791/150‑1 (MQ=255)
                                                                   tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                    <  2:537915/150‑1 (MQ=255)
                                                                   tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                    <  1:2114340/150‑1 (MQ=255)
                                                                    tttatctttGATGCGTTAGATTATCGCTATCTCACTGCAGTCTTAGCCACCATCTCTGGGGTGGAATGTTATAAACCGACGGAAAGTCTGTTTTTTTACATCCGGATTCGGACATGGCTTAATTTTACGATGACCCAGTGAATGtatata                                                                                   <  1:140210/143‑1 (MQ=255)
                                                                    aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTGTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                                                                   >  2:766101/1‑150 (MQ=255)
                                                                         cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                                                              <  2:864677/150‑1 (MQ=255)
                                                                           ttGATGCGTTAGATGATGGCTATCTCACACCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAAGCAAAGTCGGTTTTTTTACTTCCTGATTCGGAAAAGGCTTAATATTAATATGACACACTGAAAGAAAATAAATcagc                                                                            >  1:1606889/1‑147 (MQ=18)
                                                                            tGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                                                                           >  2:2116395/1‑150 (MQ=255)
                                                                             gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACATGGATTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                                                                          >  1:1035676/1‑150 (MQ=255)
                                                                              aTGCGTTAGATGATTTCTATCTCAATCCAGTCAGAGCCAACAACTCAGGTCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                                                                         <  1:341832/150‑1 (MQ=37)
                                                                              aTGCGTTAGATGATGTCTATCTCACTCCAGTCAGATCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                                                                         <  2:1038193/150‑1 (MQ=255)
                                                                                  gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                                                                     <  1:2045984/150‑1 (MQ=255)
                                                                                   ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  2:288415/150‑1 (MQ=255)
                                                                                   ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  1:353763/150‑1 (MQ=255)
                                                                                   ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  2:148183/150‑1 (MQ=255)
                                                                                   ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  1:2116395/150‑1 (MQ=255)
                                                                                    tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGTAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATCTGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                                                                   <  1:372398/150‑1 (MQ=255)
                                                                                     aGATGATGGCTATCTCCCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTATTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatat                                                                  >  1:1361417/1‑150 (MQ=255)
                                                                                       atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                                                                >  2:954297/1‑150 (MQ=255)
                                                                                          atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                                                             <  2:1838221/150‑1 (MQ=255)
                                                                                           tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCg                                                            <  1:1932389/150‑1 (MQ=255)
                                                                                             gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGaa                                                          >  2:1372036/1‑150 (MQ=255)
                                                                                              cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                         <  1:1266021/150‑1 (MQ=255)
                                                                                                 tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc                                                      >  2:1123719/1‑150 (MQ=255)
                                                                                                 tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc                                                      >  1:580517/1‑150 (MQ=255)
                                                                                                     aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                  <  1:1577907/150‑1 (MQ=255)
                                                                                                      cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcact                                                 >  2:741187/1‑150 (MQ=255)
                                                                                                         cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCg                                              <  1:1250380/150‑1 (MQ=255)
                                                                                                               gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCt                                        <  1:1803525/150‑1 (MQ=255)
                                                                                                                 gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACCTCCGGATTCGGACATGGCTTATTAAGACGATGACCGAGTGAAAGTATATAAATCGTGACTGAGATATATCCCGAAGTGCTCCCTCCGCCAGCTgg                                      >  2:1994428/1‑149 (MQ=25)
                                                                                                                   caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGTCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGTTATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    <  1:535166/150‑1 (MQ=255)
                                                                                                                   caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    >  1:1318130/1‑150 (MQ=255)
                                                                                                                   aactaaCTCAGGTATGGAAAGTATAAAACCGACTCAAAGTAGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    <  1:292728/146‑1 (MQ=255)
                                                                                                                     ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaagaa                                  <  2:123184/150‑1 (MQ=255)
                                                                                                                        aCTCAGGGCTTGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGTCCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATc                               <  1:801457/150‑1 (MQ=255)
                                                                                                                        aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACTAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGGGCTCCCTCCGCCAGCTGAAGAAATc                               <  2:794311/150‑1 (MQ=255)
                                                                                                                         cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCg                              >  2:737317/1‑150 (MQ=255)
                                                                                                                           cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCt                            <  1:1223837/150‑1 (MQ=255)
                                                                                                                            aGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAAGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTa                           <  2:1140954/150‑1 (MQ=255)
                                                                                                                             gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTaa                          <  1:2167986/150‑1 (MQ=255)
                                                                                                                             gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAACCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGGGCTCCCTCCGCCATCTGAAGAAATCGCTaa                          >  2:1059465/1‑150 (MQ=255)
                                                                                                                                cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTc                       <  1:1783464/150‑1 (MQ=255)
                                                                                                                                  ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                     >  2:1361174/1‑150 (MQ=255)
                                                                                                                                    gatgtaaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGAGTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCGCTGCGATTTATTCCGTAGTGCTCCCTCCGCCAGCTGAAGAAATGGCTAATTCTTg                    >  1:1308317/7‑150 (MQ=255)
                                                                                                                                    gatgtaaaaaaaCCGACTTTGCGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   >  1:1192199/7‑150 (MQ=255)
                                                                                                                                    aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   <  2:1216837/150‑1 (MQ=255)
                                                                                                                                    aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   <  1:1937197/150‑1 (MQ=255)
                                                                                                                                    aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   <  1:752152/150‑1 (MQ=255)
                                                                                                                                     aaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                  <  2:725518/150‑1 (MQ=255)
                                                                                                                                     aaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                  <  1:389371/150‑1 (MQ=255)
                                                                                                                                       gTAGAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                >  1:791803/1‑150 (MQ=255)
                                                                                                                                       gTAAAAAACCGACGCAAAGTCGGTTATTTTACATCCGGATTCGGACAAGCCTTAATATTACGATGACCCAGTGAAAGTATATAAATCGTCAATGCGATATATACCGAAGTGCTGCCTCCGCCAGCTGATGAAATAGCTAATTCTTGCAAt                >  2:1850721/1‑150 (MQ=255)
                                                                                                                                        tAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAtt               >  2:2199612/1‑149 (MQ=255)
                                                                                                                                          aaaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt              >  2:1771026/1‑150 (MQ=255)
                                                                                                                                          aaaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt              >  1:666297/1‑150 (MQ=255)
                                                                                                                                                ccGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCa        >  1:519485/1‑150 (MQ=255)
                                                                                                                                                      cAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt  <  1:1709258/150‑1 (MQ=255)
                                                                                                                                                      cAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt  <  2:1318130/150‑1 (MQ=255)
                                                                                                                                                      cAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGATTAATATGAAGATGACACAGTGAAAGGATATAAATCGTCACAGCGATATATAACGAAGTGCTCCCTCCGCAAAATGAAGAAATCGCTAATTATTGCAATGTTAGCCAATGGCt  >  2:1471767/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647980‑3648278

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: