Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,231,009 G→C 8.1% G219A (GGG→GCG)  rluF → 23S rRNA pseudouridine(2604) and tRNA(Tyr) pseudouridine(35) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,231,0090GC 8.1% 463.7 / 26.0 148G219A (GGG→GCG) rluF23S rRNA pseudouridine(2604) and tRNA(Tyr) pseudouridine(35) synthase
Reads supporting (aligned to +/- strand):  ref base G (60/76);  new base C (5/7);  total (65/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-02

Only 100 of 156 total aligned reads displayed.

TAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAA  >  NC_000913/4230860‑4231157
                                                                                                                                                     |                                                                                                                                                    
 aaaaaaGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTgggg                                                                                                                                                     <  1:2022103/150‑1 (MQ=255)
         agCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccg                                                                                                                                             <  2:1957533/150‑1 (MQ=255)
         agCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCAGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAATCTCGAACGCACGCGAATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccg                                                                                                                                             <  2:2144538/150‑1 (MQ=255)
          gcgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGgcgc                                                                                                                                            <  2:733552/150‑1 (MQ=255)
          gcgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgc                                                                                                                                            <  1:1402204/150‑1 (MQ=255)
          gcgcCGTTTGGCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGgcgc                                                                                                                                            <  2:2485656/150‑1 (MQ=255)
           cgcCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgct                                                                                                                                           <  2:797683/150‑2 (MQ=255)
           cgcCGGTTGTCTTCCGCATTACCCTGGTGCACGGGCTGAACCGTCAGATCCGGCGCATGGGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGCCATTCCCCCAGCGGAATGccgct                                                                                                                                           <  2:2013602/150‑2 (MQ=255)
             ccGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGAt                                                                                                                                         <  1:240835/150‑1 (MQ=255)
                tttGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTa                                                                                                                                      <  2:2165395/150‑1 (MQ=255)
                  tGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAc                                                                                                                                    >  2:663814/1‑150 (MQ=255)
                  tGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAc                                                                                                                                    >  1:1331240/1‑150 (MQ=255)
                   gTCTTCCGCATTACCCTGGTTCAGTGTCTGAACCGTGTGATGCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGATCGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAcc                                                                                                                                   <  1:514251/150‑1 (MQ=255)
                   gTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAAcc                                                                                                                                   <  2:1755664/150‑1 (MQ=255)
                    tCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcg                                                                                                                                  <  1:1293644/150‑1 (MQ=255)
                      ttCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcgac                                                                                                                                <  2:2328247/150‑1 (MQ=255)
                      ttCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcgac                                                                                                                                <  1:96635/150‑1 (MQ=255)
                      ttCCGCATTACCCTGGCGCAGGGGCTGAACCGTCCGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTAAACcgac                                                                                                                                <  2:1352286/150‑1 (MQ=255)
                       tCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACcgacg                                                                                                                               <  1:973928/150‑1 (MQ=255)
                          gCATTACCATGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATg                                                                                                                            <  2:1439785/150‑1 (MQ=255)
                            aTTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAg                                                                                                                          >  2:1225069/1‑150 (MQ=255)
                            aTTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAg                                                                                                                          >  1:1005760/1‑150 (MQ=255)
                                cccTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTaa                                                                                                                      >  1:1750244/1‑150 (MQ=255)
                                cccTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTaa                                                                                                                      >  2:2016762/1‑150 (MQ=255)
                                  cTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATc                                                                                                                    <  2:2451681/150‑1 (MQ=255)
                                   tGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCg                                                                                                                   <  1:1981878/150‑1 (MQ=255)
                                      tGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGAcc                                                                                                                <  2:1878927/150‑1 (MQ=255)
                                      tGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGAcc                                                                                                                <  2:2400497/150‑1 (MQ=255)
                                       gCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACct                                                                                                               >  1:1787858/1‑150 (MQ=255)
                                          ggggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGATGTGAATAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGTCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACGTCtt                                                                                                            <  2:1532304/150‑1 (MQ=255)
                                          ggggCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCTCGATTTAACCGACGATGAGTTAATCGACCTCtt                                                                                                            <  2:541243/150‑1 (MQ=255)
                                              cTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAg                                                                                                        >  1:277232/1‑150 (MQ=255)
                                               tGAACCGTCAGATCCGGCGCATGTTCGAGCATTTCGGCTATGAATTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGc                                                                                                       <  2:2337555/150‑1 (MQ=255)
                                                  aCCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCa                                                                                                    <  1:1522709/150‑1 (MQ=255)
                                                    cGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                                  <  2:1707659/150‑1 (MQ=255)
                                                    cGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                                  >  1:439873/1‑150 (MQ=255)
                                                     gTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTg                                                                                                 <  1:1204757/150‑1 (MQ=255)
                                                       cAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaa                                                                                               >  1:523218/1‑150 (MQ=255)
                                                        aGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaa                                                                                              <  1:1144615/150‑1 (MQ=255)
                                                        aGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaa                                                                                              >  1:585436/1‑150 (MQ=255)
                                                         gATCCGGCGCATGTGCGAGCATTTCCGCTATGAAGTCAAAAAGCTGCAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCGGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGaaaa                                                                                             <  2:1670446/150‑1 (MQ=255)
                                                          aTCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAAt                                                                                            >  1:1314866/1‑150 (MQ=255)
                                                                       gCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAgg                                                                               <  1:246556/150‑1 (MQ=255)
                                                                        cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGt                                                                              >  1:858668/1‑150 (MQ=255)
                                                                        cGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCAGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCAtt                                                                                                  >  1:1154894/1‑130 (MQ=255)
                                                                         gAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTa                                                                             <  2:582261/150‑1 (MQ=255)
                                                                          aGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTaa                                                                            >  2:543699/1‑150 (MQ=255)
                                                                            cATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTaaaa                                                                          >  1:225283/1‑150 (MQ=255)
                                                                              tttCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAAcc                                                                        <  1:1464790/150‑1 (MQ=255)
                                                                                 cGGCTATGAAGTTAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTaa                                                                     <  1:1268894/150‑1 (MQ=255)
                                                                                    cTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGc                                                                  <  1:858811/150‑1 (MQ=255)
                                                                                       tGAAGTGAAAAAGCTGGAACGCACTAGCATCATGAACGTTAGCTTAAGCGGCTTTCCTGTGGGGGAATGGCGCGATTTAACCGCCGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGATGTAAAACCTAAAGCGaa                                                               <  2:1015252/150‑1 (MQ=255)
                                                                                        gAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAg                                                              >  1:325780/1‑150 (MQ=255)
                                                                                           gTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGcc                                                           >  1:600886/1‑150 (MQ=255)
                                                                                            tGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCa                                                          >  1:1038622/1‑150 (MQ=255)
                                                                                            tGAAAAAGCTGGAACGCACGAGCATCATGAACGTTAGCTTAAGCGGCATTCCGGTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACATAAAGCGAAGGCCa                                                          >  2:821901/1‑150 (MQ=255)
                                                                                             gAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaa                                                         >  1:575847/1‑150 (MQ=255)
                                                                                             gAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCaa                                                         >  1:1017184/1‑150 (MQ=255)
                                                                                               aaaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCCGCATTCCGCTGGGGGAATCGCGCGATTTAACCGAAGATGAGTTAATCGACATCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAc                                                       >  1:731538/1‑150 (MQ=255)
                                                                                                aaaGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAAcc                                                      <  1:1775274/150‑1 (MQ=255)
                                                                                                      ggAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAAc                                                >  2:951855/1‑150 (MQ=255)
                                                                                                           gcacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCgg                                           <  2:1331240/150‑1 (MQ=255)
                                                                                                            cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCggg                                          >  2:595350/1‑150 (MQ=255)
                                                                                                            cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCggg                                          <  1:2430786/150‑1 (MQ=255)
                                                                                                            cacgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCggg                                          >  2:1476287/1‑150 (MQ=255)
                                                                                                             acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGc                                         <  1:2411404/150‑1 (MQ=255)
                                                                                                             acgcGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGc                                         <  1:1599877/150‑1 (MQ=255)
                                                                                                                cgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATc                                      <  2:585436/150‑1 (MQ=255)
                                                                                                                cgcATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATc                                      <  1:2290537/150‑1 (MQ=255)
                                                                                                                    tcatGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAAc                                  >  1:2198896/1‑150 (MQ=255)
                                                                                                                     catGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACg                                 <  2:1012361/150‑1 (MQ=255)
                                                                                                                     catGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACg                                 <  1:257689/150‑1 (MQ=255)
                                                                                                                          aCGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAg                            <  2:1822180/150‑1 (MQ=255)
                                                                                                                            gTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTc                          <  2:1059897/150‑1 (MQ=255)
                                                                                                                                gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                      >  1:1681935/1‑150 (MQ=255)
                                                                                                                                gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                      >  1:698348/1‑150 (MQ=255)
                                                                                                                                gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                      >  2:762328/1‑150 (MQ=255)
                                                                                                                                gCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTa                      >  2:1419234/1‑150 (MQ=255)
                                                                                                                                 cTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTaa                     >  2:1097641/1‑150 (MQ=255)
                                                                                                                                 aaatcgcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCGAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTaa                     >  1:1599511/6‑150 (MQ=255)
                                                                                                                                  ttAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAg                    >  2:1707293/1‑150 (MQ=255)
                                                                                                                                  ttAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAg                    <  2:65024/150‑1 (MQ=255)
                                                                                                                                  aatcgcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAg                    >  1:391990/5‑150 (MQ=255)
                                                                                                                                    tcgcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGAt                  >  2:1095690/3‑150 (MQ=255)
                                                                                                                                      gcgcCATTCCCCCAGCGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATgg                >  2:2381331/1‑150 (MQ=255)
                                                                                                                                       cGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGa               >  1:1538374/1‑150 (MQ=255)
                                                                                                                                        ggCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaa              <  1:352714/150‑1 (MQ=255)
                                                                                                                                         gCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaa             >  1:1754819/1‑150 (MQ=255)
                                                                                                                                         gCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGaaa             >  1:2363308/1‑150 (MQ=255)
                                                                                                                                               cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAATCAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcgg       <  1:1001702/150‑1 (MQ=255)
                                                                                                                                               cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcgg       >  1:1069974/1‑150 (MQ=255)
                                                                                                                                               cGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGAAGGAAAAAAcgg       >  2:923113/1‑150 (MQ=255)
                                                                                                                                                gCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggc      >  1:1251623/1‑150 (MQ=255)
                                                                                                                                                 cTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggcg     <  2:271177/150‑1 (MQ=255)
                                                                                                                                                 cTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggcg     <  1:685133/150‑1 (MQ=255)
                                                                                                                                                  tGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggcgg    <  1:1080028/150‑1 (MQ=255)
                                                                                                                                                  tGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggcgg    <  2:1251623/150‑1 (MQ=255)
                                                                                                                                                  tGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAAcggcgg    <  2:600886/150‑1 (MQ=255)
                                                                                                                                                   gggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGa   <  2:1787858/150‑1 (MQ=255)
                                                                                                                                                    ggggAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGaa  <  1:1627633/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAAAAAAGAAGCGCCGTTTGTCTTCCGCATTACCCTGGTGCAGGGGCTGAACCGTCAGATCCGGCGCATGTGCGAGCATTTCGGCTATGAAGTGAAAAAGCTGGAACGCACGCGCATCATGAACGTTAGCTTAAGCGGCATTCCGCTGGGGGAATGGCGCGATTTAACCGACGATGAGTTAATCGACCTCTTTAAGCTCATTGAAAATTCCTCTTCCGAGGTAAAACCTAAAGCGAAGGCCAAACCGAAAACAGCGGGCATCAAACGTCCAGTCGTTAAGATGGAAAAAACGGCGGAA  >  NC_000913/4230860‑4231157

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