Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 59,625 G→A 7.3% intergenic (+346/+62) yabP → / ← rluA pseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,6250GA 7.3% 258.2 / 7.6 95intergenic (+346/+62)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/50);  new base A (3/4);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAA  >  NC_000913/59478‑59767
                                                                                                                                                   |                                                                                                                                              
aGACCCGTTAAAATGAAATATAAGAGACGGTCACCGGGTGAAGAAAAAGGTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACCACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGa                                                                                                                                               >  1:650285/1‑149 (MQ=255)
aGACCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGGAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGa                                                                                                                                               >  1:87183/1‑149 (MQ=255)
  aCCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACt                                                                                                                                             >  1:842594/1‑149 (MQ=255)
     cctaaaaaTGAAAACTAAGAGACGGGCAACGGGTGAAGAAAAAGTTCAAAAAATCGCTGTGGAGCAGGAACGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCgg                                                                                                                                          <  2:694177/145‑1 (MQ=255)
       ttAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCgggg                                                                                                                                        >  1:160956/1‑149 (MQ=255)
           aaTTAAATATAAGAGACGGTCTACGTGTGAAGAAAAAGTTCAAAAATTCGCTGTGGTGCAGGAAGTGTATTACCGAATGGATAGCGTAGCCAAACGCAACAACTGATAGCAGTTTGGCATAAACAAAATATCCCCTGACTCGGGGAttt                                                                                                                                    <  2:961547/149‑1 (MQ=255)
           aaTGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAAATACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGAttt                                                                                                                                    <  1:183211/149‑1 (MQ=255)
             tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCCCACGCAACAACTGCAAGCAGTTTGGCAGAATCAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                  >  2:594149/1‑149 (MQ=255)
             tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                  <  2:109201/149‑1 (MQ=255)
             tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                  >  2:491887/1‑149 (MQ=255)
             tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGATCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                  <  1:477990/149‑1 (MQ=255)
              ggaaTATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                 <  2:135163/147‑4 (MQ=255)
              gAATTATAAGTGACGGTCCCCGGGCGAAGAAACAGTTCCACAAATCGCTGTGGAGCCGGAAGGGACTTACCGAAGGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                 <  2:587217/149‑4 (MQ=255)
              gAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTGCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                 <  2:531011/149‑4 (MQ=255)
              gAAATATAAGAGACGGTCAACGGGAGAAGAAAAAGTTCAAACATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                 <  1:883682/149‑4 (MQ=255)
               aaaTATAAGAGCCGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACACAAAATCCCCGGACTCGGGGATTTATGt                                                                                                                                <  1:773041/149‑1 (MQ=255)
               aaaTATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGt                                                                                                                                <  2:350793/149‑1 (MQ=255)
                aaTATAAGAGACTGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCATGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTa                                                                                                                               <  2:778842/149‑1 (MQ=255)
                    taAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGTAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAag                                                                                                                           <  1:899860/149‑1 (MQ=255)
                     aaGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAaga                                                                                                                          >  2:459910/1‑149 (MQ=255)
                       cagaCGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGCGCAGGAAGGGAATAACCGAATGGAAAGCGTAGCCACACGCAACAACAGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAgg                                                                                                                        <  1:83013/148‑1 (MQ=255)
                        cgaCGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGGGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGc                                                                                                                       <  2:392659/148‑1 (MQ=255)
                             gTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGccc                                                                                                                  >  1:55702/1‑149 (MQ=255)
                                 aCGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCCCCCGCAACAACTGAACGAAGTATGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGCGGCAGCCCTTAg                                                                                                              >  1:441737/1‑149 (MQ=255)
                                    ggTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAAGTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGGTTGGCAGAAACAAACAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGAt                                                                                                           <  2:601444/149‑1 (MQ=255)
                                      tGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGa                                                                                                         <  1:14009/149‑1 (MQ=255)
                                       gaagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTGTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAg                                                                                                        <  1:696338/149‑1 (MQ=255)
                                        aagaaAAAGTTCAAAAATTCGCTGTGGCGCAGGAAGGGAATTACCGAATTGAAAGCGTAGCGACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCTGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAgg                                                                                                       <  1:574288/149‑1 (MQ=255)
                                          gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGTATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGt                                                                                                     >  2:178023/1‑149 (MQ=255)
                                           aaaaaGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCAGACGCAACAACTGAAAGCAGTTTGGCATAAACAAAAAATCACCGGACTCGGGGATTTATGTACAATAGGCAGCCCTTAGGATGAGGGta                                                                                                    >  1:685359/1‑149 (MQ=255)
                                             aaaGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGtata                                                                                                  >  1:917990/1‑149 (MQ=255)
                                             aaaGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGtata                                                                                                  >  1:377730/1‑149 (MQ=255)
                                             aaaGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGtata                                                                                                  >  2:381175/1‑149 (MQ=255)
                                                   cAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTa                                                                                            <  1:722883/149‑1 (MQ=255)
                                                    aaaaaaTCGCTGTGTAGCAGGAAGGGAATGACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTAc                                                                                           <  2:349066/149‑1 (MQ=255)
                                                    aaaaaTTCGCTGTGGAGCAGGAAGGGATTTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAAGTACAAGAGGCAGCCCTTAGGATGAGGGTAGAAACGTAc                                                                                           >  2:842819/1‑149 (MQ=255)
                                                    aaaaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCGCACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTAc                                                                                           <  2:873544/149‑1 (MQ=255)
                                                     aaaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGACAGCGTAGCCACACGCAACAACTGACCGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGTGTATAAACGTACa                                                                                          <  2:880849/149‑1 (MQ=255)
                                                     aaaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCCCCCGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGCCGCGGGGATTAATGTACAAGAGGCAGAACTTAGGATGAGCGTATAAACGTACa                                                                                          >  2:211387/1‑149 (MQ=255)
                                                     aaaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAAc                                                                                               <  2:200327/144‑1 (MQ=255)
                                                     aaaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAAc                                                                                               >  1:200327/1‑144 (MQ=255)
                                                     aaaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACa                                                                                          <  1:1025142/149‑1 (MQ=255)
                                                      aaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAg                                                                                         <  2:19032/149‑1 (MQ=255)
                                                            gCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAgg                                                                                   <  2:173645/149‑1 (MQ=255)
                                                             cTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGt                                                                                  <  1:381175/149‑1 (MQ=255)
                                                              tgtgGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGtt                                                                                 <  1:275247/149‑1 (MQ=255)
                                                              tgtgGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGtt                                                                                 >  1:297314/1‑149 (MQ=255)
                                                                tgGAGCAGGAAGGGAATTACCGAATGCAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTaa                                                                               <  1:990964/149‑1 (MQ=255)
                                                                  gAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTaaaa                                                                             <  2:322195/149‑1 (MQ=255)
                                                                    tcAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGACACAAAAATTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAAt                                                                           <  2:430224/148‑1 (MQ=255)
                                                                        gAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAACCAAAAAATCCCCGGACTCGGGGCTTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGc                                                                       >  2:701538/1‑149 (MQ=255)
                                                                         aaGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAATGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGTATTTATGTACAAGAGGCAGCCGTTAGGATGAGGGTATAAACGTACAGGAATGGTTAAAAATCCGCt                                                                      <  2:387509/149‑1 (MQ=255)
                                                                             gagttACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGgcg                                                                  >  1:401939/4‑149 (MQ=255)
                                                                               aTTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTTGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAATGTTAAAAATCCGCTGGCGCt                                                                <  2:118591/149‑1 (MQ=255)
                                                                                 tACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCttt                                                              <  2:186133/149‑1 (MQ=255)
                                                                                  aCCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTa                                                             <  2:158027/149‑1 (MQ=255)
                                                                                     gaaaggaaagCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGAGTTATATACAACAGGAAGCCCTTAGGATTAGGGAATAAACGTACAGGATTGGTTAAAAATCCGCTGGCGCTTTAAAc                                                          >  2:595469/1‑149 (MQ=255)
                                                                                     gAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAAc                                                          >  1:973567/1‑149 (MQ=255)
                                                                                         ggAAAGCGTAGCCACACGCAAAAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCa                                                      <  1:1005462/149‑1 (MQ=255)
                                                                                           aaaGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATa                                                    <  2:296763/149‑1 (MQ=255)
                                                                                                gTAGCCACTCGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGTGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTTTAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTAtt                                               <  1:923192/149‑1 (MQ=255)
                                                                                                    ccACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCa                                           >  1:380692/1‑149 (MQ=255)
                                                                                                     cacaCGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCAt                                          >  2:194013/1‑149 (MQ=255)
                                                                                                      acacGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATa                                         >  2:902033/1‑149 (MQ=255)
                                                                                                      acacGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATa                                         >  2:384483/1‑149 (MQ=255)
                                                                                                          gCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATATACGTACAGGAAAGGTTAAAAATCCGCTGGCGCGTTAAACGTCATACTATTGCCATACGcc                                     >  2:615640/1‑149 (MQ=255)
                                                                                                           caacaaCTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCg                                    <  2:973567/149‑1 (MQ=255)
                                                                                                            aacaacTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCgg                                   >  1:805039/1‑149 (MQ=255)
                                                                                                               aacTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATg                                <  2:160956/149‑1 (MQ=255)
                                                                                                                acTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGATGGAATAAACGTACAGGAAATGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGACGGATgg                               >  1:818730/1‑149 (MQ=255)
                                                                                                                acTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                               <  1:178023/149‑1 (MQ=255)
                                                                                                                 cTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATggg                              <  2:377730/149‑1 (MQ=255)
                                                                                                                      aGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATc                         <  1:836616/149‑1 (MQ=255)
                                                                                                                       gCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCAGGATGGGTAATCg                        >  2:588509/1‑149 (MQ=255)
                                                                                                                         aGTTTGTCAGAAACAAAAAATGCCCGGACTCGGTGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTc                      <  2:736310/149‑1 (MQ=255)
                                                                                                                          gTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCa                     >  2:472113/1‑149 (MQ=255)
                                                                                                                           tttGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGCGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCaa                    <  2:660891/149‑1 (MQ=255)
                                                                                                                            ttGGCAGAAACAAAAAATCACCTGACTCGGGGATTTATGTACAAGAGGCAGCCCTTATGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAAc                   <  2:76707/149‑1 (MQ=255)
                                                                                                                                cAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtct               <  2:87183/149‑1 (MQ=255)
                                                                                                                                cAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtct               >  2:396751/1‑149 (MQ=255)
                                                                                                                                 aGAATCAAAAAATCCCCGGAATCGGGGATTTAAGTACAAGAGGCAGCCCTTAGTATGAGGGTATTAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAATCGTCATACTATTGCACTACGCCGGATGGGTAATCGTCCACAtctc              >  2:943112/1‑149 (MQ=255)
                                                                                                                                  gAAACAAAAAATCCCCGGCCTCGGGGATTTATGTACAAGCGGCCGCCCTTAGGATGAGGGTATAAACCTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATGCGCCGGAGGGGTAATCGGCAACATctct             >  1:262829/1‑149 (MQ=255)
                                                                                                                                  gAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATctct             >  2:736517/1‑149 (MQ=255)
                                                                                                                                  gAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGAAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATctct             >  2:677518/1‑149 (MQ=255)
                                                                                                                                    accAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTTGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGc           <  1:170658/147‑1 (MQ=255)
                                                                                                                                    aaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGc           <  1:1004315/149‑1 (MQ=255)
                                                                                                                                      cAAAAAATCCCCGGACTCGGGGATTTATGTACTAGAGGCAGCCATTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTGTAAACGTCATACTATTGCCATACGCAGGATGGGTTATCGTCAACATCTCTGCAt         >  1:616540/1‑149 (MQ=255)
                                                                                                                                      cAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACATGAAAGGTTAAAATTCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCAt         >  1:105759/1‑149 (MQ=255)
                                                                                                                                       cataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCGTTAGGATGAGGGTATTAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCTTg        >  1:814766/4‑149 (MQ=255)
                                                                                                                                        ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCgtgt       >  2:941268/3‑149 (MQ=255)
                                                                                                                                        aaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGt       <  1:372124/149‑1 (MQ=255)
                                                                                                                                         taaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTa      >  1:41255/2‑149 (MQ=255)
                                                                                                                                         aaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTa      <  1:999918/149‑1 (MQ=255)
                                                                                                                                           aaTCCCCGGACTCGGGGATTTATGTACTAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTCATCGTCATACTATTGCCTTACGCCGGTTGGGTATTCGTCAACATGTCTGCATGTAGc    <  2:194797/149‑1 (MQ=255)
                                                                                                                                             tCCACGGTCTGGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGcaa  <  1:68784/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGACCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAA  >  NC_000913/59478‑59767

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: