Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,132 A→C 5.4% noncoding (70/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,1320AC 5.4% 396.7 / 8.1 130noncoding (70/84 nt)agrAinactive antisense sRNA
Reads supporting (aligned to +/- strand):  ref base A (49/74);  new base C (4/3);  total (53/77)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

Only 100 of 132 total aligned reads displayed.

GATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAA  >  NC_000913/3647985‑3648280
                                                                                                                                                    |                                                                                                                                                    
gATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGt                                                                                                                                                    >  1:148183/1‑150 (MQ=255)
   aaTCATACAGTACATGCAGGTTATAAAACCAGCAAGTCCTTGCAATAGTTTCAGTTTGGTTTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCTGAGCCACCAACTCAGGGCTTGAAAGT‑AAAAAACCGACGCAAAGTCgg                                                                                                                                                 <  1:1765188/150‑1 (MQ=38)
    aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGt                                                                                                                                                >  1:725518/1‑150 (MQ=255)
       aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttt                                                                                                                                             <  1:948672/150‑1 (MQ=255)
         aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                           <  1:1017502/150‑1 (MQ=255)
         aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                           <  2:1660951/150‑1 (MQ=255)
          cAGTACATGCAGTTTATAAAACCAGCACGTCCTTGCAAAAGTTTCAGTATGTTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTGAGAGCCACCAACTCAGTGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                          <  2:249278/150‑1 (MQ=37)
          cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                          >  1:2124538/1‑150 (MQ=255)
           aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGA‑AAAAAACCGACGCAAAGTCGGtttttatt                                                                                                                                         >  2:1878783/1‑147 (MQ=255)
           aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGAGGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑TAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                         <  1:653630/150‑3 (MQ=37)
            gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                         <  2:990900/150‑3 (MQ=38)
               cgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAAGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATc                                                                                                                                     >  2:1937197/3‑150 (MQ=37)
                aTGCAGTTTATAAAACCAGCACGTCCTTGCAAAAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttactttc                                                                                                                                     <  1:2115433/150‑7 (MQ=37)
                 tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCg                                                                                                                                   >  2:1266021/1‑150 (MQ=255)
                   cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGa                                                                                                                                 >  2:1709258/1‑150 (MQ=255)
                   cAGGTTATAAAACCAGCACCTCCTTGCAATTCTATCAGTAAGGAATTAGCATTGAAGCGTAAGATTATGGCTAAAACAATCCAGACAAATCCCCCAACACAGCGATGGAAAGA‑AAAAAACCGACGAAAAGTCCATTTTTATAcgacctca                                                                                                                                 >  2:422972/1‑143 (MQ=255)
                    aGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGAt                                                                                                                                <  1:994643/150‑1 (MQ=255)
                    aGGTTATAAAACCAGCACGTCCTTGCAAAAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCACAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGAt                                                                                                                                <  1:2002170/150‑1 (MQ=255)
                     ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACAGACTAAAAGTCGTTTTTTTTACATCCGTAta                                                                                                                               >  2:341832/1‑149 (MQ=37)
                       ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCg                                                                                                                             <  1:1958687/150‑1 (MQ=255)
                         ataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGa                                                                                                                           <  1:156802/150‑1 (MQ=255)
                           aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACa                                                                                                                         >  2:2038618/1‑150 (MQ=255)
                           aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACa                                                                                                                         >  2:389371/1‑150 (MQ=255)
                            aaaCTAGCAACTTCTTGCAATATTTTCCTTATGGTATTAGAATTGATGTGTTAGATTAGTTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcct                                                                                                                        <  2:480722/150‑3 (MQ=11)
                                cAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCt                                                                                                                    <  2:86759/150‑1 (MQ=255)
                                     cGTCATTGCAATAGTATCAGTATGGGATTAGCATTGATGCGGTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCACCTCAGGGCTGGTAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                                                               <  2:326083/150‑1 (MQ=25)
                                        ccTTGCAATAGTTTCATTATGGTATTAGCATTGTTTCGTTAGATGATGGCTATCTCAATCCAGTCCGAGCCAACAACTCAGGGATTGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                            <  1:214649/150‑1 (MQ=37)
                                        ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                            >  2:890261/1‑150 (MQ=255)
                                         cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACACCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                           <  1:779244/150‑1 (MQ=255)
                                          ttGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                                                                          <  1:1397811/150‑1 (MQ=255)
                                             cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACATCCATTAATATGACGaag                                                                                                       >  1:932870/1‑148 (MQ=18)
                                             cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                                                                       <  2:42412/150‑1 (MQ=255)
                                             cAATAGTTTCAGTATGGTATTAGCATTGATGAGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGTCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATACGGAAAAGCTTTAAAATGACGATg                                                                                                       >  1:1077628/1‑150 (MQ=25)
                                              aaTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                      <  2:822714/150‑1 (MQ=255)
                                              aaTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                      <  1:408007/150‑1 (MQ=255)
                                                tAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                                                                    <  1:1491912/150‑1 (MQ=255)
                                                     tCAGGATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                                                               <  1:374794/150‑1 (MQ=255)
                                                        gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                                                                            >  1:1140954/1‑150 (MQ=255)
                                                         tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                                                                           <  1:511063/150‑1 (MQ=255)
                                                         tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAACAAGTGaaat                                                                                           >  2:193249/1‑149 (MQ=255)
                                                             gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                                                                       <  2:954791/150‑1 (MQ=255)
                                                              tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                      <  2:537915/150‑1 (MQ=255)
                                                              tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                      <  1:2114340/150‑1 (MQ=255)
                                                               tttatctttGATGCGTTAGATTATCGCTATCTCACTGCAGTCTTAGCCACCATCTCTGGGGTGGAATGT‑TATAAACCGACGGAAAGTCTGTTTTTTTACATCCGGATTCGGACATGGCTTAATTTTACGATGACCCAGTGAATGtatata                                                                                     <  1:140210/143‑1 (MQ=255)
                                                               aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTGTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                                                                     >  2:766101/1‑150 (MQ=255)
                                                                    cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                                                                <  2:864677/150‑1 (MQ=255)
                                                                      ttGATGCGTTAGATGATGGCTATCTCACACCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGAAGCAAAGTCGGTTTTTTTACTTCCTGATTCGGAAAAGGCTTAATATTAATATGACACACTGAAAGAAAATAAATcagc                                                                              >  1:1606889/1‑147 (MQ=18)
                                                                       tGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                                                                             >  2:2116395/1‑150 (MQ=255)
                                                                        gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACAGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACATGGATTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                                                                            >  1:1035676/1‑150 (MQ=255)
                                                                         aTGCGTTAGATGATTTCTATCTCAATCCAGTCAGAGCCAACAACTCAGGTCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                                                                           <  1:341832/150‑1 (MQ=37)
                                                                         aTGCGTTAGATGATGTCTATCTCACTCCAGTCAGATCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                                                                           <  2:1038193/150‑1 (MQ=255)
                                                                             gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                                                                       <  1:2045984/150‑1 (MQ=255)
                                                                              ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                      <  1:353763/150‑1 (MQ=255)
                                                                              ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                      <  2:148183/150‑1 (MQ=255)
                                                                              ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                      <  2:288415/150‑1 (MQ=255)
                                                                              ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCACGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                      <  1:301838/150‑1 (MQ=255)
                                                                                aGATGATGGCTATCTCCCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTATTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatat                                                                    >  1:1361417/1‑150 (MQ=255)
                                                                                  atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                                                                  >  2:954297/1‑150 (MQ=255)
                                                                                     atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                                                               <  2:1838221/150‑1 (MQ=255)
                                                                                      tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCg                                                              <  1:1932389/150‑1 (MQ=255)
                                                                                        gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGaa                                                            >  2:1372036/1‑150 (MQ=255)
                                                                                         cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                           <  1:1266021/150‑1 (MQ=255)
                                                                                            tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc                                                        >  2:1123719/1‑150 (MQ=255)
                                                                                                aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                    <  1:1577907/150‑1 (MQ=255)
                                                                                                 cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcact                                                   >  2:741187/1‑150 (MQ=255)
                                                                                                   ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc                                                 <  1:1950489/150‑1 (MQ=255)
                                                                                                    cAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCg                                                <  1:1250380/150‑1 (MQ=255)
                                                                                                      gTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGcc                                              <  1:825954/150‑1 (MQ=255)
                                                                                                          gagCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCt                                          <  1:1803525/150‑1 (MQ=255)
                                                                                                            gCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACCTCCGGATTCGGACATGGCTTATTAAGACGATGACCGAGTGAAAGTATATAAATCGTGACTGAGATATATCCCGAAGTGCTCCCTCCGCCAGCTgg                                        >  2:1994428/1‑149 (MQ=25)
                                                                                                              caccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGTCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGTTATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                      <  1:535166/150‑1 (MQ=255)
                                                                                                              caccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                      >  1:1318130/1‑150 (MQ=255)
                                                                                                                ccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaagaa                                    <  2:123184/150‑1 (MQ=255)
                                                                                                                   aCTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACTAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGGGCTCCCTCCGCCAGCTGAAGAAATc                                 <  2:794311/150‑1 (MQ=255)
                                                                                                                    cTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCg                                >  2:737317/1‑150 (MQ=255)
                                                                                                                      cAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCt                              <  1:1223837/150‑1 (MQ=255)
                                                                                                                       aGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAAGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTa                             <  2:1140954/150‑1 (MQ=255)
                                                                                                                        gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTaa                            <  1:2167986/150‑1 (MQ=255)
                                                                                                                        gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAACCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGGGCTCCCTCCGCCATCTGAAGAAATCGCTaa                            >  2:1059465/1‑150 (MQ=255)
                                                                                                                           cTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTc                         <  1:1783464/150‑1 (MQ=255)
                                                                                                                             ggAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                       >  2:1361174/1‑150 (MQ=255)
                                                                                                                               gatgtaaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGAGTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCGCTGCGATTTATTCCGTAGTGCTCCCTCCGCCAGCTGAAGAAATGGCTAATTCTTg                      >  1:1308317/7‑150 (MQ=255)
                                                                                                                               gatgtaaaaaaaCCGACTTTGCGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                     >  1:1192199/7‑150 (MQ=255)
                                                                                                                               aaaGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                     <  1:1937197/150‑1 (MQ=255)
                                                                                                                               aaaGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                     <  2:1216837/150‑1 (MQ=255)
                                                                                                                                aaGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                    <  1:389371/150‑1 (MQ=255)
                                                                                                                                aaGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                    <  2:725518/150‑1 (MQ=255)
                                                                                                                                  gT‑AGAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                  >  1:791803/1‑150 (MQ=255)
                                                                                                                                  gT‑AAAAAACCGACGCAAAGTCGGTTATTTTACATCCGGATTCGGACAAGCCTTAATATTACGATGACCCAGTGAAAGTATATAAATCGTCAATGCGATATATACCGAAGTGCTGCCTCCGCCAGCTGATGAAATAGCTAATTCTTGCAAt                  >  2:1850721/1‑150 (MQ=255)
                                                                                                                                   taaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGCCCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                  >  1:1899136/3‑150 (MQ=255)
                                                                                                                                   taaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                  >  1:2106100/3‑150 (MQ=255)
                                                                                                                                   t‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAtt                 >  2:2199612/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt                >  1:666297/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt                >  2:1771026/1‑150 (MQ=255)
                                                                                                                                        aaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAg             >  1:75936/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)
                                                                                                                                                   aaGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTaa  >  2:58551/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAA  >  NC_000913/3647985‑3648280

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: ‑