Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 59,628 T→C 7.6% intergenic (+349/+59) 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,6280TC 7.6% 243.3 / 10.3 92intergenic (+349/+59)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (34/50);  new base C (3/4);  total (38/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.53e-01

ACCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACGTGG  >  NC_000913/59480‑59776
                                                                                                                                                    |                                                                                                                                                    
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:491887/1‑149 (MQ=255)
           tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                           <  2:109201/149‑1 (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:384483/1‑149 (MQ=255)
                                                                                                    acacGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATa                                                  >  2:902033/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)
                                                                                                                                                    tCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAACgtgg  >  1:798984/1‑149 (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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑