Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 59,627 C→T 9.5% intergenic (+348/+60) yabQ → / ← rluA protein YabQ/23S rRNA pseudouridine(746) and tRNA pseudouridine(32) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,6270CT 9.5% 333.3 / 18.9 116intergenic (+348/+60)yabQ/rluAprotein YabQ/23S rRNA pseudouridine(746) and tRNA pseudouridine(32) synthase
Reads supporting (aligned to +/- strand):  ref base C (52/52);  new base T (8/3);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 117 total aligned reads displayed.

ACCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGG  >  NC_000913/59480‑59776
                                                                                                                                                   |                                                                                                                                                     
 cccGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCg                                                                                                                                                    >  1:1287523/1‑150 (MQ=255)
   cGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCggg                                                                                                                                                  <  2:715707/150‑1 (MQ=255)
   cGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCggg                                                                                                                                                  >  1:1424281/1‑150 (MQ=255)
    gTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCgggg                                                                                                                                                 >  1:372134/1‑150 (MQ=255)
      tAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGTGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGAt                                                                                                                                               >  1:589029/1‑150 (MQ=255)
         aaTGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTa                                                                                                                                            <  1:924092/150‑1 (MQ=255)
          aTGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                           <  2:591918/150‑1 (MQ=255)
          aTGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                           >  2:468355/1‑150 (MQ=255)
           tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATg                                                                                                                                          >  1:109636/1‑150 (MQ=255)
           tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATg                                                                                                                                          <  2:1055509/150‑1 (MQ=255)
           tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATg                                                                                                                                          <  1:106984/150‑1 (MQ=255)
            gAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAttttttg                                                                                                                                         <  2:1268952/150‑5 (MQ=255)
                tataAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACaa                                                                                                                                     >  2:695366/1‑150 (MQ=255)
                  taAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAaga                                                                                                                                   >  1:523831/1‑150 (MQ=255)
                  taAGAGACGCTCAACGGGTGAAGAAAAAGTACAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACACGGGGATTTATGTACAaga                                                                                                                                   >  2:1473287/1‑150 (MQ=255)
                      agaCGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCa                                                                                                                               <  1:124978/150‑1 (MQ=255)
                         cGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGcc                                                                                                                            >  1:68721/1‑150 (MQ=255)
                         cGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAttttttgtt                                                                                                                                       <  2:263224/139‑7 (MQ=255)
                         cGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAttttttgtt                                                                                                                                       >  1:263224/1‑133 (MQ=255)
                          ggTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGccc                                                                                                                           <  2:702089/150‑1 (MQ=255)
                          ggTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGccc                                                                                                                           >  1:1247912/1‑150 (MQ=255)
                             cAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTa                                                                                                                        >  2:1055523/1‑150 (MQ=255)
                              aaCGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAg                                                                                                                       <  1:1254943/150‑1 (MQ=255)
                                    agaagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGCAATTACCGAAAGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAg                                                                                                                 >  2:66460/2‑150 (MQ=255)
                                     gaagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAgg                                                                                                                >  2:166631/1‑150 (MQ=255)
                                      aagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAggg                                                                                                               <  2:1379111/150‑1 (MQ=255)
                                       agaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGt                                                                                                              >  1:90982/1‑150 (MQ=255)
                                                 cAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTAc                                                                                                    >  2:1505913/1‑150 (MQ=255)
                                                     aaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGa                                                                                                <  2:1497691/150‑1 (MQ=255)
                                                      aTTCGCTGTGGAGCAGGAAGGGAATAACCGAATGGAAAGCGTAGCCACACGCGACAACTGAAAGCAGTTAGGCAGAAACAAAAAATCCCCCGACTCGGGGATTTATGAACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGaa                                                                                               >  1:1108705/1‑150 (MQ=255)
                                                        tCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAg                                                                                             <  1:887571/150‑1 (MQ=255)
                                                          gCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGt                                                                                           <  1:1036956/150‑1 (MQ=255)
                                                          gCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGt                                                                                           >  1:1001661/1‑150 (MQ=255)
                                                             gtgGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGAACAGGAAAGGTTaa                                                                                        >  2:1520791/1‑150 (MQ=255)
                                                              tgGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTaaa                                                                                       <  1:467718/150‑1 (MQ=255)
                                                               ggAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTaaaa                                                                                      >  1:1451716/1‑150 (MQ=255)
                                                                  gCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATc                                                                                   <  1:749119/150‑1 (MQ=255)
                                                                  gCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATc                                                                                   <  1:695366/150‑1 (MQ=255)
                                                                   cAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATcc                                                                                  <  2:68721/150‑1 (MQ=255)
                                                                    aGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCg                                                                                 >  2:1199538/1‑150 (MQ=255)
                                                                    aGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCg                                                                                 <  2:1287523/150‑1 (MQ=255)
                                                                       aaGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTg                                                                              <  1:723068/150‑1 (MQ=255)
                                                                          ggAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGgcg                                                                           >  2:1078133/1‑150 (MQ=255)
                                                                           gAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGgcgc                                                                          >  1:957538/1‑150 (MQ=255)
                                                                            aaTTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATTTACAATAGGCAGCCCTTAGTATGAGTGTATAAACTTACAGGAAAGTTTAAAAATCCGATGGCGCt                                                                         <  2:1252297/150‑1 (MQ=255)
                                                                             aTTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCtt                                                                        <  1:468355/150‑1 (MQ=255)
                                                                              ttACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCttt                                                                       <  1:166631/150‑1 (MQ=255)
                                                                                 ccGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTaaa                                                                    <  1:504438/150‑1 (MQ=255)
                                                                                   gAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATc                                                                                   <  2:1295867/133‑1 (MQ=255)
                                                                                   gAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACg                                                                  <  1:1505913/150‑1 (MQ=255)
                                                                                     aTGGAAAGCGTAGCGACACGAAACAACTGAAAGCAGTTTGGCAGAAACAATAATTCCCCTTACTCGGGGATTTATGTACAAGAGGCAGCCGTTAGGTTGATGGTATAATGGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTc                                                                <  2:240615/150‑1 (MQ=255)
                                                                                      tGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCa                                                               <  1:985593/150‑1 (MQ=255)
                                                                                        gAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATa                                                             >  2:714659/1‑150 (MQ=255)
                                                                                         aaaGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATAc                                                            >  1:776044/1‑150 (MQ=255)
                                                                                         aaaGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATAc                                                            >  2:419537/1‑150 (MQ=255)
                                                                                            gCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTAt                                                         >  2:975818/1‑150 (MQ=255)
                                                                                                aGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGcc                                                     >  1:624063/1‑150 (MQ=255)
                                                                                                aGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGcc                                                     <  2:586985/150‑1 (MQ=255)
                                                                                                   cacaCGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATa                                                  >  2:1494099/1‑150 (MQ=255)
                                                                                                    acacGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATAc                                                 <  1:258179/150‑1 (MQ=255)
                                                                                                         caacaaCTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCgg                                            >  1:847891/1‑150 (MQ=255)
                                                                                                         caacaaCTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCgg                                            >  1:83642/1‑150 (MQ=255)
                                                                                                            caacTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATg                                         >  2:1451057/1‑150 (MQ=255)
                                                                                                             aacTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                                        >  1:249878/1‑150 (MQ=255)
                                                                                                               cTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGt                                      >  1:169553/1‑150 (MQ=255)
                                                                                                               cTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGt                                      <  1:1078133/150‑1 (MQ=255)
                                                                                                                 gAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTaa                                    <  2:1247912/150‑1 (MQ=255)
                                                                                                                    aGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGTATGAGTGTATAAACGTACAGTAAATGTTAAATATCCTCTGGCGCTTTAAACGTGATACTATTGCCATACGCCTGATGGGTAATCg                                 <  2:1504260/150‑1 (MQ=255)
                                                                                                                    aGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGAGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCg                                 <  1:172877/150‑1 (MQ=255)
                                                                                                                            ggCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtc                         >  1:133659/1‑150 (MQ=255)
                                                                                                                              cAGAAACAAAAAATCGCCGGACTCGGGGATTTATGTACAAGCGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGCTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtctc                       <  2:819361/150‑1 (MQ=255)
                                                                                                                              cAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtctc                       >  1:1181172/1‑150 (MQ=255)
                                                                                                                              cAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtctc                       <  2:516679/150‑1 (MQ=255)
                                                                                                                                gAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTg                     >  2:908451/1‑150 (MQ=255)
                                                                                                                                gAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTg                     <  2:433906/150‑1 (MQ=255)
                                                                                                                                 aaaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGc                    <  2:83642/150‑1 (MQ=255)
                                                                                                                                  aaCAAAAAATTCCCGGAATCGGGGTTTTAAGTAAAATATGGATCCCTTATGTAGTTGGGTTTAACGTTCATGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCa                   <  2:625748/150‑1 (MQ=255)
                                                                                                                                  aaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCa                   <  1:155246/150‑1 (MQ=255)
                                                                                                                                  aaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCa                   <  2:1001661/150‑1 (MQ=255)
                                                                                                                                  aaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCa                   <  1:836580/150‑1 (MQ=255)
                                                                                                                                  aaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCa                   >  1:1352764/1‑150 (MQ=255)
                                                                                                                                   tacataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCAt                  >  1:320202/6‑150 (MQ=255)
                                                                                                                                    cAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATg                 >  1:582207/1‑150 (MQ=255)
                                                                                                                                    cAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATg                 >  2:982052/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATg                 <  1:1421894/149‑1 (MQ=255)
                                                                                                                                     cataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGt                >  1:482059/4‑150 (MQ=255)
                                                                                                                                     cataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGt                >  1:1389425/4‑150 (MQ=255)
                                                                                                                                     cataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGt                >  1:187252/4‑150 (MQ=255)
                                                                                                                                     aaaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCCACATCTCTGCATGt                >  2:576353/1‑150 (MQ=255)
                                                                                                                                     aaaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACATACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGt                <  2:387862/150‑1 (MQ=255)
                                                                                                                                      ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTa               >  1:323272/3‑150 (MQ=255)
                                                                                                                                      aaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTa               <  2:918067/150‑1 (MQ=255)
                                                                                                                                      aaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTa               >  1:543001/1‑150 (MQ=255)
                                                                                                                                       taaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAg              >  1:1089743/2‑150 (MQ=255)
                                                                                                                                       aaagTCCCCGGACTCGGGTATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAg              <  1:675450/150‑1 (MQ=255)
                                                                                                                                        aaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGc             >  1:476525/1‑150 (MQ=255)
                                                                                                                                             cccGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGcaacaa        >  2:490760/1‑150 (MQ=255)
                                                                                                                                             cccGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGcaacaa        >  2:91576/1‑150 (MQ=255)
                                                                                                                                               cGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAAc      <  2:1451716/150‑1 (MQ=255)
                                                                                                                                                  aCTCGGGGATTTATGTACAAGAGGGAGCAGTTAGGATGAGGGTATAAAAGGACATGAAAGGTTAAATATCCGCTGGCTCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGGATTTAGCAACAAACgt  <  2:1479070/150‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGG  >  NC_000913/59480‑59776

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: