Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,449,020 A→T 7.9% R50R (CGT→CGA)  rplP ← 50S ribosomal subunit protein L16

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,449,0200AT 7.9% 257.4 / 10.8 76R50R (CGT→CGA) rplP50S ribosomal subunit protein L16
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base T (3/3);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ACCCGGCTGAATCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTA  >  NC_000913/3448873‑3449165
                                                                                                                                                   |                                                                                                                                                 
aCCCGGCTGAATCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCacgac                                                                                                                                                  <  4:459475/149‑1 (MQ=255)
aCCCGGCTGAATCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACCGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCacgac                                                                                                                                                  >  3:68384/1‑149 (MQ=255)
  ccGGCTGAATCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGt                                                                                                                                                <  2:629327/149‑1 (MQ=255)
    ggCTGAATCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGtgc                                                                                                                                              <  1:701325/149‑1 (MQ=255)
         aaTCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGTCGtgctatga                                                                                                                                         <  1:436016/149‑5 (MQ=255)
         aaTCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGTCGtgctatga                                                                                                                                         <  2:754800/149‑5 (MQ=255)
          aTCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTc                                                                                                                                        >  3:175505/1‑149 (MQ=255)
           tCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTTAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCg                                                                                                                                       >  1:5374/1‑149 (MQ=255)
            cAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATGCGCACTGCCAGCGGCTTTTCAGTGATCGGTGTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGa                                                                                                                                      <  2:125674/149‑1 (MQ=255)
              aGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATc                                                                                                                                    <  1:486125/149‑1 (MQ=255)
              aGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATc                                                                                                                                    >  2:298021/1‑149 (MQ=255)
                 cAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGa                                                                                                                                 <  2:706061/149‑1 (MQ=255)
                    cccAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGTATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACgg                                                                                                                              <  2:239565/149‑1 (MQ=255)
                       aaTACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCa                                                                                                                           <  1:568702/149‑1 (MQ=255)
                       aaTACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCa                                                                                                                           >  3:321111/1‑149 (MQ=255)
                         tACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGt                                                                                                                         >  3:650638/1‑149 (MQ=255)
                            tCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAg                                                                                                                      >  1:748670/1‑149 (MQ=255)
                            tCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAg                                                                                                                      >  4:300376/1‑149 (MQ=255)
                            tCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAg                                                                                                                      >  2:145667/1‑149 (MQ=255)
                             ccACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAga                                                                                                                     >  2:415893/1‑149 (MQ=255)
                              cACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAgac                                                                                                                    <  4:321111/149‑1 (MQ=255)
                                      ctttacctttacCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACgg                                                                                                            <  4:137721/149‑1 (MQ=255)
                                           cctttaccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAAc                                                                                                       <  3:467936/149‑1 (MQ=255)
                                            ctttaccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACa                                                                                                      >  2:746136/1‑149 (MQ=255)
                                            ctttaccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACa                                                                                                      >  1:635523/1‑149 (MQ=255)
                                            ctttaacCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCGAACa                                                                                                      <  4:472534/149‑1 (MQ=255)
                                             tttaccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACa                                                                                                      <  3:44738/148‑1 (MQ=255)
                                             tttaccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACa                                                                                                      >  4:44738/1‑148 (MQ=255)
                                               taccCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCt                                                                                                   <  1:516192/149‑1 (MQ=255)
                                                  ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGa                                                                                                                                 >  1:633442/1‑116 (MQ=255)
                                                  ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGa                                                                                                                                 <  2:633442/116‑1 (MQ=255)
                                                  ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTc                                                                                                >  4:680420/1‑149 (MQ=255)
                                                  ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTc                                                                                                >  4:141514/1‑149 (MQ=255)
                                                  ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTc                                                                                                >  2:256934/1‑149 (MQ=255)
                                                  ccATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATATTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTc                                                                                                >  4:286069/1‑149 (MQ=255)
                                                   cATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCa                                                                                               <  4:210541/149‑1 (MQ=255)
                                                      aCGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGAc                                                                                            <  3:549786/149‑1 (MQ=255)
                                                        gCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAgac                                                                                                                    <  4:205559/123‑1 (MQ=255)
                                                        gCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAgac                                                                                                                    >  3:205559/1‑123 (MQ=255)
                                                               cAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGc                                                                                   >  1:539376/1‑149 (MQ=255)
                                                               cAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGCCCAACAGCTTTCAGACCGAAGCTGc                                                                                   >  3:496139/1‑149 (MQ=255)
                                                                     cTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGc                                                                             <  2:272800/149‑1 (MQ=255)
                                                                              gATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCg                                                                    >  1:218520/1‑149 (MQ=255)
                                                                               aTCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTc                                                                                                                        >  3:654937/1‑96 (MQ=255)
                                                                               aTCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTc                                                                                                                        <  4:654937/96‑1 (MQ=255)
                                                                               aTCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGt                                                                   >  3:569510/1‑149 (MQ=255)
                                                                               aTCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGt                                                                   >  4:459677/1‑149 (MQ=255)
                                                                                tCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTa                                                                  <  2:539376/149‑1 (MQ=255)
                                                                                        tCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTgcgc                                                          <  1:135533/149‑1 (MQ=255)
                                                                                          cGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCa                                                        <  1:367600/149‑1 (MQ=255)
                                                                                            ggAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGa                                                      >  1:603957/1‑149 (MQ=255)
                                                                                                  cacGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCgg                                                >  2:584420/1‑149 (MQ=255)
                                                                                                   acGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGt                                               >  3:696854/1‑149 (MQ=255)
                                                                                                    cGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTACGCCAGACCGCGGtt                                              <  2:748670/149‑1 (MQ=255)
                                                                                                       aTCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACg                                           >  3:254780/1‑149 (MQ=255)
                                                                                                        tCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACgg                                          >  3:191089/1‑149 (MQ=255)
                                                                                                         ccAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCt                                                                                                                                   >  2:602997/1‑59 (MQ=255)
                                                                                                         ccAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCt                                                                                                                                   <  1:602997/59‑1 (MQ=255)
                                                                                                         ccAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGc                                         >  2:19346/1‑149 (MQ=255)
                                                                                                           aGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCACCAGACCGCGGTTACGGCCt                                       >  3:622263/1‑149 (MQ=255)
                                                                                                             aTCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCttt                                     <  2:635523/149‑1 (MQ=255)
                                                                                                             aTCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCttt                                     <  2:5374/149‑1 (MQ=255)
                                                                                                                 tACCTCGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGc                                 <  3:164077/149‑1 (MQ=255)
                                                                                                                       gACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTa                           <  1:59175/149‑1 (MQ=255)
                                                                                                                            ttAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTa                           <  3:588817/144‑1 (MQ=255)
                                                                                                                            ttAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTa                           >  4:588817/1‑144 (MQ=255)
                                                                                                                                   cACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGGAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTa               <  1:298021/149‑1 (MQ=255)
                                                                                                                                   cACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTa               <  4:191089/149‑1 (MQ=255)
                                                                                                                                        gaagcaGCACGTCGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCtt          >  1:229998/6‑149 (MQ=255)
                                                                                                                                        gTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCtt          >  4:156278/1‑149 (MQ=255)
                                                                                                                                          agcaGCACGTCGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTg        >  2:449232/4‑149 (MQ=255)
                                                                                                                                           aTAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACgg                                          <  3:392217/114‑1 (MQ=255)
                                                                                                                                           aTAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACgg                                          >  4:392217/1‑114 (MQ=255)
                                                                                                                                           aTAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTgg       <  2:141471/149‑1 (MQ=255)
                                                                                                                                             aGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGtt     >  3:731695/1‑149 (MQ=255)
                                                                                                                                              gCACGTCGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTtgt                                   <  1:113000/118‑1 (MQ=255)
                                                                                                                                              gCACGTCGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTtgt                                   >  2:113000/1‑118 (MQ=255)
                                                                                                                                              gCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTg    >  2:643337/1‑149 (MQ=255)
                                                                                                                                                acgacgTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTa  >  2:86504/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACCCGGCTGAATCAAGGCAACCCAATACTCCACGTTACCTTTACCTTTACCCATACGCACTGCCAGCGGCTTTTCAGTGATCGGTTTGTCCGGGAACACACGGATCCAGATCTTACCTTGACGCTTAACTGCACGGGTCATAGCACGACGTGCTGCTTCGATCTGACGGGCAGTCAGACGACCACGGCCAACAGCTTTCAGACCGAAGCTGCCGAAGCTAACATCCGTACCCTGCGCCAGACCGCGGTTACGGCCTTTGTGCATTTTACGGAATTTTGTACGCTTTGGTTGTA  >  NC_000913/3448873‑3449165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑