Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,129 C→T 9.7% noncoding (67/84 nt) agrA → small RNA AgrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,1290CT 9.7% 497.7 / 29.5 177noncoding (67/84 nt)agrAsmall RNA AgrA
Reads supporting (aligned to +/- strand):  ref base C (65/91);  new base T (6/11);  total (74/103)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

AACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647980‑3648278
                                                                                                                                                     |                                                                                                                                                     
aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc                                                                                                                                                       >  1:1829042/1‑150 (MQ=255)
 aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCa                                                                                                                                                      <  2:1095532/150‑1 (MQ=255)
   tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaaa                                                                                                                                                    >  2:1781824/1‑150 (MQ=255)
   tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaaa                                                                                                                                                    >  2:202807/1‑150 (MQ=255)
    ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAg                                                                                                                                                   >  1:493129/1‑150 (MQ=255)
         aTCATACAGAACATGCACGTTATAAAACCAGCAATTCCTTGCAATTTTTTCAGTTTGGTATTAGTATTGATGTGTTAGATTATTTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                                                              <  2:193920/150‑1 (MQ=21)
              aCAGCAGATGCACGCTATAAAAACAGCACCTCCTTGCAATCTTTTCCCCTTGGTAATAGCATTGATGCGTTAGATCATCCCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttt                                                                                                                                         <  1:2359052/150‑1 (MQ=12)
               cAGTACTTGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  2:353689/150‑1 (MQ=255)
               cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                                                                        <  1:316105/150‑1 (MQ=255)
               cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                         >  1:773938/1‑149 (MQ=255)
               cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                                                                         <  2:773938/149‑1 (MQ=255)
                aGTACATGCCCGTTATAAAACCAGCTCTTCCTTGCAATATTTTCATCATGGTATTAGCATTGATGCGTTAGATTATCCCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                       <  1:81887/150‑3 (MQ=21)
                aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttac                                                                                                                                       <  1:248589/150‑3 (MQ=38)
                aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                                                                       <  2:307241/150‑1 (MQ=255)
                 gTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTAc                                                                                                                                      <  2:1332479/150‑1 (MQ=255)
                  tACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACa                                                                                                                                     <  2:1582609/150‑1 (MQ=255)
                    cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttactttc                                                                                                                                   <  1:1075458/150‑7 (MQ=38)
                     gtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATcc                                                                                                                                  >  2:1291244/2‑150 (MQ=255)
                          ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGAATTGATGCGTTAGATGATTGCTATCTCACTCCAGTCAGAGCCACAAAATCAGTGCTGGAAAGTAAAAAACCTACGCAAAGTCGGTTTTTTTACATCCgcctt                                                                                                                             >  2:2022465/1‑146 (MQ=37)
                             tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCgg                                                                                                                          <  1:159095/150‑1 (MQ=255)
                              ataAAACCCGCACGCCCTTGCACTCGTTCCAGTATGGTATTAGCATTCATGCGCTAGATGATGGCTATCTCCCTCCAGTCCGCGGCACCAACTCAGGGCGGTAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGa                                                                                                                         <  2:2439432/150‑1 (MQ=25)
                               taAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGAc                                                                                                                        <  1:1752739/150‑1 (MQ=255)
                                aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAACTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACa                                                                                                                       <  1:1214046/150‑1 (MQ=255)
                                    acAGCACCTCCTTGCAATTTTTTCCCTTTGGTATTAGCATTGATGCGTTAGATGATTGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcctcct                                                                                                                   <  2:1593590/149‑6 (MQ=17)
                                        cACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                                                               >  1:692975/1‑150 (MQ=255)
                                        cACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGATTaa                                                                                                               >  1:582169/1‑150 (MQ=255)
                                              cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                         >  1:2158299/1‑150 (MQ=255)
                                                tGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                                                                       >  2:365807/1‑150 (MQ=255)
                                                  cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                                                                     >  2:1228814/1‑150 (MQ=255)
                                                    aTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                                                                   >  2:2474497/1‑150 (MQ=255)
                                                      aGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                                                                 <  2:1815546/150‑1 (MQ=255)
                                                      aGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                                                                 <  1:450653/150‑1 (MQ=255)
                                                      aGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                                                                 <  2:2473960/150‑1 (MQ=255)
                                                       gTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                                                                >  2:92606/1‑150 (MQ=255)
                                                        tttCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCATGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                                                               <  1:1943153/150‑1 (MQ=255)
                                                        tttCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                                                               <  1:1291244/150‑1 (MQ=255)
                                                         ttCAGTATGGTATAAGCATTGATGCGTTTGATGAGGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                                                                              <  2:1440388/150‑1 (MQ=255)
                                                          tCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                                                             <  2:84750/150‑1 (MQ=255)
                                                           cAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGa                                                                                            >  1:1166229/1‑150 (MQ=255)
                                                             gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                                                                          >  2:1270463/1‑150 (MQ=255)
                                                             gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                                                                          >  1:1743140/1‑150 (MQ=255)
                                                              tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg                                                                                                                                               <  1:943598/96‑1 (MQ=255)
                                                               aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGt                                                                                        <  1:661648/150‑1 (MQ=255)
                                                                tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  2:2219359/150‑1 (MQ=255)
                                                                tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  2:1829042/150‑1 (MQ=255)
                                                                tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  1:1571858/150‑1 (MQ=255)
                                                                tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  2:1342346/150‑1 (MQ=255)
                                                                tGGTATTAGCATTGATGCGTTAGATGATGGATATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                       <  2:225079/150‑1 (MQ=255)
                                                                  gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                                                                     <  2:2299403/150‑1 (MQ=255)
                                                                   tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                    >  1:172476/1‑150 (MQ=255)
                                                                        gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                               <  1:365807/150‑1 (MQ=255)
                                                                        gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                               <  1:1765490/150‑1 (MQ=255)
                                                                         cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                                                              >  2:224232/1‑150 (MQ=255)
                                                                         cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                                                              <  2:493129/150‑1 (MQ=255)
                                                                                gCGTTAGATGATGGCCATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                                                                       <  1:997979/150‑1 (MQ=255)
                                                                                   ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                                                                    <  1:481056/150‑1 (MQ=255)
                                                                                    tAGATTATTTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                                                                   <  2:2374008/150‑1 (MQ=255)
                                                                                    tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAGCTCAGAACAGCAAAGTAAAAAAAAAAAGAAAAGTCGGGTTTATAACAGCCGGATTCGCAGAACGATTAATATGAAGATGACCCAGTGAAAGTATATAAATCGTCGCTGCggtt                                                                   >  2:1461595/1‑147 (MQ=255)
                                                                                    tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                                                                   >  1:2095060/1‑150 (MQ=255)
                                                                                      gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                                                                 >  1:2446898/1‑150 (MQ=255)
                                                                                         gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                                                              <  1:1962919/150‑1 (MQ=255)
                                                                                              cTATCTAACTCCAGTCAGAGCCACCATCTCATGGATTTAAAGTAAAAATCCGACGCAAAGTAGGTTTTTTTACATGCGGAGTCGTACATGGCTTAATATGACGATGACCGTGTGTAAGTATATAAATCGTCACTGGGATTTATTCCGAAg                                                         <  2:2088657/150‑1 (MQ=255)
                                                                                                 tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc                                                      <  2:596858/150‑1 (MQ=255)
                                                                                                  ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTTCGGTATATATCGCAGTgac                                                     <  1:2378929/150‑3 (MQ=255)
                                                                                                   tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTCACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc                                                    >  2:441316/1‑150 (MQ=255)
                                                                                                    cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc                                                   <  2:2395397/150‑1 (MQ=255)
                                                                                                    cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGAGATATATACCGAAGTGctcc                                                   >  2:785056/1‑150 (MQ=255)
                                                                                                     aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                  <  2:1743140/150‑1 (MQ=255)
                                                                                                       tCCAGTCAGAGCCACCAACTCAGAACTGGAAAGTAAAAAAACGACGAAAATTCGGTTTTTTTACATCCGGATTCGGAGAAGGATTAATATAACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctc                                                >  2:1428287/1‑150 (MQ=255)
                                                                                                        ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATACACCGAAGTGctccctcc                                               >  1:168565/1‑150 (MQ=255)
                                                                                                         cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                         >  2:1412270/1‑139 (MQ=255)
                                                                                                         cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                         <  1:1412270/139‑1 (MQ=255)
                                                                                                         cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCg                                              <  1:2220655/150‑1 (MQ=255)
                                                                                                           gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGcc                                            <  1:232863/150‑1 (MQ=255)
                                                                                                               gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCt                                        <  1:1012626/150‑1 (MQ=255)
                                                                                                                agCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTg                                       <  2:16390/150‑1 (MQ=255)
                                                                                                                   caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                    <  1:2474497/150‑1 (MQ=255)
                                                                                                                   caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATCAAGATGACCCAGTCAATGAATATAAATCGTCACTGAGATATATACAGAAGTGCTCCCTCCCCAAGCTgaat                                    >  1:871690/1‑149 (MQ=25)
                                                                                                                    accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaaga                                   <  2:738090/150‑1 (MQ=255)
                                                                                                                    accaACTATTGGCTGGAAATTAAAAAAGCGACGCAATGTCGGTTTTTTTACATCCTGATTAGGATAAGGCTTAATATGACGATGACGAAGTTAAAGTATATAAATCGTCACTTCGATATATACGGAATTGCTCCCTCAGCCAGCTgaagt                                   <  2:2040899/150‑2 (MQ=255)
                                                                                                                          tCAGGGCTGGAAAGTAAAAAACCGACTCAAAGTCTGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGc                             <  1:1270463/150‑1 (MQ=255)
                                                                                                                            aGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTa                           <  1:673162/150‑1 (MQ=255)
                                                                                                                             gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTaa                          <  1:2111923/150‑1 (MQ=255)
                                                                                                                             gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTaa                          >  1:1223771/1‑150 (MQ=255)
                                                                                                                             gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCAGACAACGATTAAAATGACGATCAACCAGAAAAAGTATAGAAATAAACAATGCGATATATAACCAAGTGGTACCTCAGCAAGCTGAAGAAATCGATaa                          >  2:1250391/1‑150 (MQ=16)
                                                                                                                                  ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                     >  1:208398/1‑150 (MQ=255)
                                                                                                                                   gatgtaaaaaaaCCGACTTTGGGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTg                    >  1:2415806/7‑150 (MQ=255)
                                                                                                                                    aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                   <  1:2022465/150‑1 (MQ=255)
                                                                                                                                     tgtaaaaaaaCCGACTTTGAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATAGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                  >  2:82466/5‑150 (MQ=255)
                                                                                                                                      aGTAAAAAACCGACGCAAAGTCGTTTTTTTTACATCCTGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCaa                 <  1:457097/150‑1 (MQ=255)
                                                                                                                                       taaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                >  2:1975791/3‑150 (MQ=255)
                                                                                                                                       taaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                >  1:671741/3‑150 (MQ=255)
                                                                                                                                        tAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATg               <  2:682661/150‑1 (MQ=255)
                                                                                                                                        aaaaaaaCCGACTTTGCGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATg               >  1:112624/2‑150 (MQ=255)
                                                                                                                                         aaaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTCAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt              >  2:1879445/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGtt             >  1:85378/1‑150 (MQ=255)
                                                                                                                                             aaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGGAATCGCTAATTCTTGCAATGTTAGc          <  2:1205472/150‑1 (MQ=255)
                                                                                                                                             aaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGc          >  1:298158/1‑150 (MQ=255)
                                                                                                                                                     cAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt  <  1:1042670/150‑1 (MQ=255)
                                                                                                                                                     cAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGAGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt  <  2:978271/150‑1 (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 < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: