Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,531 T→G 5.3% noncoding (12/120 nt) rrfF ← 5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,5310TG 5.3% 293.7 / 2.0 105noncoding (12/120 nt)rrfF5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (11/88);  new base G (4/2);  total (15/90)
Fisher's exact test for biased strand distribution p-value = 3.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGT  >  NC_000913/3423382‑3423679
                                                                                                                                                     |                                                                                                                                                    
cAAACGGCCCAGTCATCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt                                                                                                                                                      <  2:1620966/150‑1 (MQ=255)
cAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAATTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt                                                                                                                                                      <  2:1134319/150‑1 (MQ=38)
    aGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACAATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGCCCACCGCGCTACGgccg                                                                                                                                                  >  1:1832212/1‑150 (MQ=255)
     ggCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGTCAGTTCCCTACTCTCTCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTgccgc                                                                                                                                                 <  2:1310677/150‑1 (MQ=255)
      gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc                                                                                                                                                <  2:1306504/150‑1 (MQ=37)
      gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc                                                                                                                                                <  1:475525/150‑1 (MQ=37)
      gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc                                                                                                                                                >  1:830757/1‑150 (MQ=37)
      gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgcc                                                                                                                                                <  1:374248/150‑1 (MQ=37)
      gCCCAGTCTTCCGACTGAGCATTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc                                                                                                                                                    <  2:1852731/146‑1 (MQ=255)
       cccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCa                                                                                                                                               <  1:451532/150‑1 (MQ=37)
          aGTCTTCCGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGCGTTGCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGAc                                                                                                                                            <  2:1321826/150‑1 (MQ=255)
          aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGAc                                                                                                                                            <  1:591292/150‑1 (MQ=38)
           gTCTTCCGCCTGAGCCTTTTGTTTTATTTGATGGCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACa                                                                                                                                           <  2:350199/150‑1 (MQ=255)
                ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTc                                                                                                                                      >  2:1126612/1‑150 (MQ=17)
                ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTc                                                                                                                                      >  2:877511/1‑150 (MQ=17)
                ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTc                                                                                                                                      <  1:787806/150‑1 (MQ=17)
                   ccTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                                                                                                                                   <  2:1131255/149‑1 (MQ=38)
                   ccTGAGCCTTTTGTTTTATTTGATGGCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                                                                                                                                   <  1:302012/149‑1 (MQ=25)
                     tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                                                                                                                                 <  1:1430906/150‑1 (MQ=40)
                     tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                                                                                                                                 <  2:1826804/150‑1 (MQ=40)
                     tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTACCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                                                                                                                                 <  1:1601497/150‑1 (MQ=255)
                      gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                                                                                                                                <  2:1085661/150‑1 (MQ=40)
                       aGCCTTTTGTTTTATTAGATGTCTGGCAGTCCCATACTCTCGAATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGGGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                                                                                                                               <  1:439228/150‑1 (MQ=14)
                        gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                                                                                                                              <  2:1810420/150‑1 (MQ=255)
                             ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCTTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCCGCATGGGGTCAGGTGTGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc                                                                                                                         <  1:2019013/150‑1 (MQ=255)
                             ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCTCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc                                                                                                                         <  2:843490/150‑1 (MQ=255)
                              tGTTTTATTTGATGTCTGGCAGTACCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc                                                                                                                        >  1:15006/1‑150 (MQ=255)
                                     tttGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACttt                                                                                                                 <  2:2049893/150‑1 (MQ=255)
                                        gATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAc                                                                                                              <  1:629088/150‑1 (MQ=255)
                                         aTGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAAcc                                                                                                             <  2:1937486/150‑1 (MQ=255)
                                          tGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCt                                                                                                            <  1:638142/150‑1 (MQ=255)
                                             cTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACAACCATCGGCGCTACGGCGTTTCACTTCTGGGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTaaa                                                                                                         <  2:1886375/150‑1 (MQ=255)
                                               ggCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTaaaaa                                                                                                       >  1:689998/1‑150 (MQ=255)
                                                gCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAg                                                                                                      <  1:319839/150‑1 (MQ=255)
                                                gCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAg                                                                                                      <  1:1382543/150‑1 (MQ=255)
                                                 cAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAgt                                                                                                     <  2:1986342/150‑1 (MQ=255)
                                                    tgcccTACTCTCGCAGGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAgtggt                                                                                                  <  2:1221491/148‑1 (MQ=255)
                                                     tCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAgtggtg                                                                                                 <  1:164405/150‑1 (MQ=255)
                                                       ccTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCt                                                                                               <  2:1043061/150‑1 (MQ=255)
                                                        cTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTg                                                                                              <  2:2018965/150‑1 (MQ=255)
                                                         tACTCTCGCAGGGGTAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGAATGTGGTCAGGTGGTACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGa                                                                                             <  1:363925/150‑1 (MQ=18)
                                                           ctctCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACAACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATa                                                                                           <  1:221855/150‑1 (MQ=255)
                                                             ctcGCATGGGGAGACCCCACACTACCATCGGCGCTACGGTGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGTACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATAcc                                                                                         <  1:884482/150‑1 (MQ=255)
                                                              tcGCATGCGGAGACCCCACACTTCCATCGTCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACAACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAGTTTAACCTAAAAAGTGGTGCTGATAccc                                                                                        <  1:1781037/150‑1 (MQ=18)
                                                                  aTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCagag                                                                                    <  1:905437/150‑1 (MQ=255)
                                                                        agaCCCCACACTACCATCTGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGGCGAAc                                                                              <  1:135052/150‑1 (MQ=255)
                                                                        agaCCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAAc                                                                              <  2:934609/150‑1 (MQ=255)
                                                                          aCCCCACACTACCATCGGGGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTg                                                                            <  2:1638626/150‑1 (MQ=255)
                                                                            cccACACTACCATCGGCGCTTCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTggg                                                                          <  1:314832/150‑1 (MQ=255)
                                                                            cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTggg                                                                          <  1:479811/150‑1 (MQ=255)
                                                                            cccACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTggg                                                                          <  1:1229939/150‑1 (MQ=255)
                                                                             ccACACTATCATCTGCGCGACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTgggg                                                                         <  1:1276438/150‑1 (MQ=255)
                                                                                cacTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGAcc                                                                      <  1:172248/150‑1 (MQ=255)
                                                                                    aCCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCAc                                                                  <  1:1192538/150‑1 (MQ=255)
                                                                                      cATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTAAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCAccc                                                                <  1:425138/150‑1 (MQ=255)
                                                                                         cGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTa                                                             <  1:641390/150‑1 (MQ=255)
                                                                                              cTACGGCGTGTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCacg                                                        <  2:1310706/150‑3 (MQ=255)
                                                                                               tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCAAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAgg                                                       <  1:978039/150‑1 (MQ=255)
                                                                                                tcGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAggg                                                      <  1:1251560/149‑1 (MQ=255)
                                                                                                  ggCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTg                                                    <  1:1415520/150‑1 (MQ=255)
                                                                                                  ggCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGGGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTg                                                    <  2:342714/150‑1 (MQ=255)
                                                                                                   gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTAACCAAGGGTgc                                                   <  2:689998/150‑1 (MQ=255)
                                                                                                    cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgag                                                  <  2:1154718/150‑3 (MQ=255)
                                                                                                    cGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgag                                                  <  1:348287/150‑3 (MQ=255)
                                                                                                     gTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTTCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACGCAGAGTCGAACTGGGGACCTCACCCTTGGTAAGGGTgagg                                                 <  2:1489593/150‑4 (MQ=37)
                                                                                                      ttgcACTTCTGAGTTCGGCCTGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGct                                                <  2:262149/147‑1 (MQ=255)
                                                                                                        tCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGctct                                              <  1:1284488/150‑1 (MQ=255)
                                                                                                          ccTTCTGCGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAACAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                            <  2:1520208/149‑1 (MQ=255)
                                                                                                          aCTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                            <  2:1687073/150‑1 (MQ=255)
                                                                                                          aCTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAc                                            <  1:1776171/150‑1 (MQ=255)
                                                                                                           cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                           <  2:574981/150‑1 (MQ=255)
                                                                                                           cTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTAcc                                           <  1:1528094/150‑1 (MQ=255)
                                                                                                              cTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAAc                                        <  1:198894/150‑1 (MQ=255)
                                                                                                               tGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTTCTGCCGCCAGACAAATTCTTTTCTATTCTGCCTAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCATGGGTGCGCTCTACCAACt                                       <  2:1796059/150‑1 (MQ=255)
                                                                                                                gAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTGTAATCTGCCGAACTTTACCCTAAAAAGTGGTGCTGAAACCCAGAGTCGAAATGGGGACCTCACCCTTACCAAGGGTGCGATGTACCAACTg                                      >  1:1163037/1‑150 (MQ=38)
                                                                                                                  gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                    <  2:1847894/150‑1 (MQ=255)
                                                                                                                  gTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAg                                    >  1:847611/1‑150 (MQ=255)
                                                                                                                    tCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGcc                                  <  1:176246/150‑1 (MQ=255)
                                                                                                                     gggCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCa                                 <  1:990708/149‑1 (MQ=255)
                                                                                                                     cGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCa                                 <  2:200233/150‑1 (MQ=255)
                                                                                                                       gCAGGGGGTCAGGTGGGACCACCGCGCTACTGCCGCAAGACAATTTCTTGTCTACTCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTCCCAAGGGTGCGCTCTACCAACTGAGCCata                               <  1:998361/150‑1 (MQ=37)
                                                                                                                          gggggTCAGTTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGGCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCa                            <  1:1510707/149‑1 (MQ=255)
                                                                                                                           ggggTCAGGTGGGACCACCGCGCTACTGCCGCCAGTCAAATTCTTTTCTAATCTGCCGAACTTTAACCTAATATGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCGTATCAg                           <  1:1738815/150‑1 (MQ=255)
                                                                                                                               tCAGGTGGGACCACCGCTCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACg                       <  2:778746/150‑1 (MQ=255)
                                                                                                                                cAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAGCCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCGCCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGc                      <  1:1112786/150‑1 (MQ=255)
                                                                                                                                 aGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCt                     >  1:491546/1‑150 (MQ=255)
                                                                                                                                     gggACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAAt                 >  2:10823/1‑150 (MQ=255)
                                                                                                                                     gggACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAAt                 <  1:505877/150‑1 (MQ=255)
                                                                                                                                      ggACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAAtt                <  2:1733231/150‑1 (MQ=255)
                                                                                                                                        accaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTTCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTg              <  1:1295903/150‑1 (MQ=255)
                                                                                                                                         ccaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGa             <  2:541153/150‑1 (MQ=255)
                                                                                                                                         ccaccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGa             <  1:1608485/150‑1 (MQ=255)
                                                                                                                                          caccGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGAt            >  2:1859862/1‑150 (MQ=255)
                                                                                                                                           accGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATg           <  2:468721/150‑1 (MQ=255)
                                                                                                                                             cgcgCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGcc         <  2:550038/150‑1 (MQ=255)
                                                                                                                                              gcgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCt        >  2:1880500/1‑150 (MQ=255)
                                                                                                                                               ggcTACTGCCGCCAGACAAATTCTTATCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg       <  1:1409529/149‑1 (MQ=255)
                                                                                                                                               cgcTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTg       <  1:2061318/150‑1 (MQ=255)
                                                                                                                                                    cTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGTTGCTGATACCCAGAGGCGAACTGGGGACCTCACCCTTACCAAGGGTGCTCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGt  <  2:15006/150‑1 (MQ=255)
                                                                                                                                                    cTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGTTGCGCTCTACCAACTGAGACATATCAGCACGCTAAATTTGATGCCTGGCAGt  >  1:383851/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAGCACGCTAAATTTGATGCCTGGCAGT  >  NC_000913/3423382‑3423679

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: