Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG45.2% 87.7 / 16.1 62noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/30);  new base G (28/0);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925756‑926050
                                                                                                                                                |                                                                                                                                                         
cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGa                                                                                                                                                      <  1:348165/150‑1 (MQ=255)
  cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                                                                                                                                    <  1:995670/150‑1 (MQ=255)
  cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAc                                                                                                                                                    <  1:316236/150‑1 (MQ=255)
    gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAccc                                                                                                                                                  <  1:39457/150‑1 (MQ=255)
      cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAacccc                                                                                                                                                >  1:962728/1‑150 (MQ=255)
      cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCg                                                                                                                                               <  2:227223/151‑1 (MQ=255)
       gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAccccc                                                                                                                                                >  2:94966/1‑149 (MQ=255)
               aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGtt                                                                                                                                        >  2:574514/1‑149 (MQ=255)
                    aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGta                                                                                                                                  <  2:1012998/150‑1 (MQ=255)
                            acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACtt                                                                                                                         <  1:303183/151‑1 (MQ=255)
                             caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAGTCGAACCCCCCCTACGTTGTCGTATACACACttt                                                                                                                        >  1:786189/1‑151 (MQ=255)
                               ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCCACGGTGCCGta                                                                                                                                  >  1:349756/1‑139 (MQ=25)
                               ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtat                                                                                                                                 <  2:349784/140‑1 (MQ=255)
                                gTAACATACTCCCACCGCCCGCCAATTTCACTGCCGTACAGACAGATAAAATCCGAAAAAAAAGCTCGCCCTTTCGTACGAGCTCTTCTTTAAATATCGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACTCTTTcc                                                                                                                      <  1:822978/150‑1 (MQ=18)
                                  aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCCCCACCTTGCCCGATACACACttt                                                                                                                        >  1:491989/1‑146 (MQ=25)
                                  aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAg                                                                                                                    <  2:492026/150‑1 (MQ=255)
                                     aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTGGAAACCCCACTAAGTTGCCGTTTTTACAATTTCCAGACg                                                                                                                 >  2:117966/1‑150 (MQ=21)
                                        cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCCAACGTTGCCGTCTTCACAACTTCCAcgcg                                                                                                                 >  1:666555/1‑147 (MQ=21)
                                         tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCCGTACGTTTCCGGATACACACTTTCCAGGc                                                                                                                  >  2:580442/1‑145 (MQ=21)
                                           cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCACACGTTTCCGTATAAACACTTTTCCGGCGTGCTcc                                                                                                           >  2:621260/1‑150 (MQ=21)
                                           cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGAGGGGGGTTTGAAACCCCCCCAAGTTGCCGTATAAACACCTTCCCAGCGTggccc                                                                                                           >  1:72015/1‑146 (MQ=21)
                                            aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGAACCCCCCAACAGGTGCCGGATACAAACTTTCCAGGCggg                                                                                                               >  1:874766/1‑143 (MQ=21)
                                              ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCa                                                                                                       <  1:538329/151‑1 (MQ=255)
                                              ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATACGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCa                                                                                                       <  2:954922/151‑1 (MQ=255)
                                                   caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTTCCGCATACACACCTTTCAGGCGGGCTCCTTCAGcccc                                                                                                  >  1:998971/1‑149 (MQ=18)
                                                   caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                                  <  2:893729/151‑1 (MQ=255)
                                                       aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCgg                                                                                             <  1:945987/151‑1 (MQ=255)
                                                         gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGGGGGGGGGGATCGTACCCCCCGAACGGTGTCCTGTACACACCTTTCCCGCGTGCTTCTTTTGCCAcaccg                                                                                              >  2:746494/1‑145 (MQ=12)
                                                          ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACACCCCATAAGTTGCCGTATACACACTTTCTCGGCGTGCTCCTTCATCCACTc                                                                                                >  1:627188/1‑146 (MQ=17)
                                                                gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTc                                                                                     <  2:702480/149‑1 (MQ=255)
                                                                  cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCCACACGTTGCCGTATAAACACTTTTCCGGGGGGCCCCCTCACCCCCCCGGGCAccccccc                                                                                  >  1:475732/1‑145 (MQ=17)
                                                                  cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGGTCGAAACCCCGCGACCTTGCCGTTTACACACCTTTCCCGCGTGCTCCTTCACACACACCGGCACCTCAcc                                                                                  >  2:307044/1‑150 (MQ=17)
                                                                   gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATCGAAACCCCGCTACGTTGACGTATACACACTTTTCAGGCGGGGTCCTTCCGCCACTCGGGCACCTCAcc                                                                                  >  2:174285/1‑149 (MQ=17)
                                                                     acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGGCACCTCACCcaa                                                                               >  2:824929/1‑150 (MQ=255)
                                                                      cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTCAAACCCCCCTACGTTGGCCGATAAACACTTTTCCGGCGGGCCCCCTCAGCCCCTCCGGCAccccccccaa                                                                               >  1:175772/1‑141 (MQ=9)
                                                                          cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGttt                                                                         <  1:206293/151‑1 (MQ=255)
                                                                             aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATAAGTTGCCGGATACACACTTTTCAGGCGTGCTCCTTCAGCCACTCCGACAACTCACCCAAATGtttt                                                                        >  2:806050/1‑149 (MQ=14)
                                                                               aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTa                                                                    <  2:1019004/151‑1 (MQ=255)
                                                                                aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCGCTACGGTGCCGTTTACACACTTTTCAGGCGGGCTCCCTTAGCCACCCCGGCACCCCCCCAAATTTTTTTTCt                                                                     >  2:198915/1‑149 (MQ=14)
                                                                                aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGAGTGGTACCCCCAATACACTGGCGAATCAACACTTTTCAATGGGGGTCCTTCCACCACTCGGGAACCCCACCAAAATAttttt                                                                       >  2:437662/1‑146 (MQ=9)
                                                                                    gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCGTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACTCTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCaaa                                                               <  1:265854/151‑1 (MQ=255)
                                                                                       aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCATACGTTGCCGTATACACACCTTCCAGGCGGGCTCCTTCAGCCACTCCGACACCTCCCCAAATTGTTTTTCTACCAAAcc                                                             >  2:578465/1‑150 (MQ=14)
                                                                                       aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCATAAGTTGCCGTATATACACTTTTCAGGCGTGGTCCTTTAGCCACCCGGGCACCCCACCAAAATGTTTTTTTACCAAAcc                                                             >  2:321870/1‑150 (MQ=14)
                                                                                         aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                                                          <  2:666580/151‑1 (MQ=255)
                                                                                              aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGt                                                    <  2:365280/151‑1 (MQ=255)
                                                                                                cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATAAGTTTCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCCCCCCGACACCTCACCAAAATGTTTTTCTACCAAACCTCATGGGTg                                                   >  2:370207/1‑150 (MQ=17)
                                                                                                 tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTTCACCCCCCACACGTTTTCGTATACAAACTTTTCAGGGGTGGTCCCTCAGCCCCCCCGGCACCCCACCAAAATTTTTTTTTACCAAAACCCCTgggggg                                                  >  1:334814/1‑150 (MQ=2)
                                                                                                   gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCCTACGTTGCCGTATAAACACTTTCCAGGCGGGGTCCCTCAGCCCCTCGGGCACCCCCCCCAATTTTTTTTCTTCCCAAACCCATGGGTgg                                                  >  2:426081/1‑148 (MQ=2)
                                                                                                   gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCaa                                               <  2:311867/151‑1 (MQ=255)
                                                                                                    cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAAc                                              <  2:17905/151‑1 (MQ=255)
                                                                                                        ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACggg                                           <  2:168354/150‑1 (MQ=255)
                                                                                                           gTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTGAAACCCCCCTAATTTTAGTTCTACACACATTCTCGGCGGGCTCCTTCTGCCACCCGGGACACCCACCCAAATGTTTTTATAACTAACCTCACAGGGTGGGACGGGgcg                                        >  1:797750/1‑150 (MQ=2)
                                                                                                             aCGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGAAACCCCGATACGGTGCCGGATAAACACCTTTCAGGCGTGTTCCTTCAGCCCCCCCGGCACCCCCCCCAAATGTTTTTTTTACCAAACCCCATGGGGGGAAAgggg                                           >  2:411206/1‑146 (MQ=2)
                                                                                                               gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctact                                   <  1:425341/151‑1 (MQ=255)
                                                                                                                 gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctacta                                  <  1:715104/150‑1 (MQ=255)
                                                                                                                 gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtat                                 <  2:572573/151‑1 (MQ=255)
                                                                                                                 gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtat                                 <  2:275183/151‑1 (MQ=255)
                                                                                                                           aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTg                        <  1:621208/150‑1 (MQ=255)
                                                                                                                           aaaTATGGCGGTGAGGGGGGGAGTCGAACCCCCCCAACCGTTCCGTATACACACTTTCCCGGGGGGCTTCTTCCGCCCCCCCGACACCCCCCCCAAATTTTTTTCTTCCCAACCCCATTGGGGGGGAAGGGGCGCCCCCATAAGGAGGTgg                       >  1:344366/1‑151 (MQ=2)
                                                                                                                             atatGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                     <  1:805990/151‑1 (MQ=255)
                                                                                                                                   cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc               <  1:450649/151‑1 (MQ=255)
                                                                                                                                           gggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  2:148001/1‑148 (MQ=255)
                                                                                                                                                aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:344393/151‑1 (MQ=255)
                                                                                                                                                aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:558658/151‑1 (MQ=255)
                                                                                                                                                aTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAAGTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  2:882857/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                         
CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925756‑926050

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