Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,441,061 G→C 5.8% S205S (TCC→TCG)  rpsD ← 30S ribosomal subunit protein S4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,441,0610GC 5.8% 201.9 / 2.6 75S205S (TCC→TCG) rpsD30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base C (3/2);  total (38/37)
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

CAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTT  >  NC_000913/3440913‑3441209
                                                                                                                                                    |                                                                                                                                                    
cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg                                                                                                                                                      <  1:840739/149‑1 (MQ=255)
     tATGGACAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTGCTCTCTTTTGTACTAAGCTTTACTTGGAGTa                                                                                                                                                 >  2:680913/1‑149 (MQ=255)
        ggCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAg                                                                                                                                              <  1:839949/149‑1 (MQ=255)
         gCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc                                                                                                                                             >  2:536356/1‑149 (MQ=255)
         gCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc                                                                                                                                             >  1:319280/1‑149 (MQ=255)
          ccAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt                                                                                                                                            >  1:224757/1‑149 (MQ=255)
           cAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTc                                                                                                                                           >  1:98943/1‑149 (MQ=255)
             aaGCCACGCTCTAAAGGCGCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcga                                                                                                                                         <  2:1006448/149‑1 (MQ=255)
             aaGCCACGATCTAAAGGCTCAAGGGTCACGTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcga                                                                                                                                         <  2:886823/149‑1 (MQ=255)
              aGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgcc                                                                                                                                        >  1:231642/1‑147 (MQ=255)
                ccACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgacga                                                                                                                                      >  1:522376/1‑149 (MQ=255)
                 cACGCTCTAAAGTCTCAAGGGTCACCTTGGCGTGCGTAGAAGTCAGTTGCTCGATATCAAGCAGTCGCGGTTTTAGAAACTCTGTCACAGATCCCTGCATTGTGTCCTCTCTTTGGTACTTAGCTTTACTTGGAGTAAAGCTCGACGAt                                                                                                                                     <  2:929541/149‑1 (MQ=255)
                  aCGCTCTAAAGGCTCAAGGGGCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGAGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATc                                                                                                                                    <  2:811907/149‑1 (MQ=255)
                    gCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAg                                                                                                                                  <  1:965903/149‑1 (MQ=255)
                     ctctAAAGGCGGAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCCCGGTTTTAGATACTCTGTCACAGAACCCTGCAATGTGTCCTCTCTTTGGTACTAAGCTTAACTGGGAGTAAAGCTCGACGATCAgg                                                                                                                                 <  1:535859/149‑1 (MQ=255)
                      tctAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCATCCAGGCGCGGGTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAATCTTTTCTTGGAGTAAAGCTCGACGATCAGgt                                                                                                                                >  2:102391/1‑149 (MQ=255)
                          aaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAGTTGCTCGATATCAACCAGGCGCGGGTTTAGAAACTCTGTCACAGACCCCGGCATTGTGTGCTCTGTTTGGTACTAAGCTTTACTTGGAGTAACGCGCGCCGATCAGGTGTTc                                                                                                                            >  1:820651/1‑149 (MQ=255)
                          aaGGCTCAAGGGGCACCTTGGCGTGCGGCGAACTCACTTGCTCGATCTCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGGACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTc                                                                                                                            <  1:682282/149‑1 (MQ=255)
                           aGGCTCCAGGGTCACCTTGGCGGGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGGCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCg                                                                                                                           <  1:1016450/149‑1 (MQ=255)
                           aGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCg                                                                                                                           >  2:615501/1‑149 (MQ=255)
                            ggCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGt                                                                                                                          <  1:160242/149‑1 (MQ=255)
                                cAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAAt                                                                                                                      <  2:133276/149‑1 (MQ=255)
                                      gcACCTTGGCGTGCGTCTACCTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGc                                                                                                                <  2:30864/148‑1 (MQ=255)
                                        aCCTTGGCGTGCGTCGAACTCACTAGCTCGAAATCAACCAGGCGCGGTTTTAGAAACACTGTCAAAGAACCGAGCATTGTGAGCTCTGTTTGGTCCTAAGCTTTACTTGGAGTAAAGCTCGCCGATCAGGAGTTCGTAAAAGTCCGcag                                                                                                              >  2:657613/1‑149 (MQ=255)
                                         ccTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGCACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGTAGTTAAGCTCGACGCTCAGGTGTTCGTTAATGTCCGcaga                                                                                                             >  2:161501/1‑149 (MQ=255)
                                         ccTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcaga                                                                                                             >  1:745560/1‑149 (MQ=255)
                                          cTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagac                                                                                                            >  2:522835/1‑149 (MQ=255)
                                                gtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCa                                                                                                      <  1:932275/149‑1 (MQ=255)
                                                 tgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATTTCCGCAGACAGATCAg                                                                                                     >  2:94663/1‑149 (MQ=255)
                                                     tCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACg                                                                                                 <  1:537436/149‑1 (MQ=255)
                                                      cGAACTCACTTGCTCGATATCAACCCGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTATAGCTCGACGATCAGGTGTTCGTTAAAGTCCGCAGACAGATCAGAACGc                                                                                                >  2:990443/1‑149 (MQ=255)
                                                      cGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGc                                                                                                >  2:341337/1‑149 (MQ=255)
                                                        aaCTCACTTGCTCGATATCAACCAGGCGCGGGTTTAGAAACGCTGTCACAGACCCCTGCATGGTGTCCTCTCTTTGGTACTAAGCTTTACTTCGAGTAAAGCTCGACGAGCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTc                                                                                              <  2:601675/149‑1 (MQ=255)
                                                          cTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCg                                                                                            <  1:103921/149‑1 (MQ=255)
                                                             aCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCt                                                                                         >  2:422226/1‑149 (MQ=255)
                                                                  cTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGc                                                                                    <  2:224757/149‑1 (MQ=255)
                                                                       tatCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCTATGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTaaa                                                                               <  1:707584/149‑1 (MQ=255)
                                                                       tatCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTaaa                                                                               <  2:54249/149‑1 (MQ=255)
                                                                        atCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAAc                                                                              <  1:356770/149‑1 (MQ=255)
                                                                            aCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          >  1:292988/1‑149 (MQ=255)
                                                                             ccATGCGCGGATTTAGAAACTCTGTCACAGAACACTGCATTGTGTCGTCTCATTGGTACTAAGCTTTACTTCTTGTAAAGCTCGACGATCAGGTGTTCGTTAATGTGCGCAGACAGATCAGAACGCTCCGGCTTACTCTTAAACGTAcc                                                                         <  2:929825/149‑1 (MQ=255)
                                                                             ccAGGCGCGGTTTTAGAAACTCAGGCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTGGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCTCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAcc                                                                         <  2:950598/149‑1 (MQ=255)
                                                                                 gcgcGGTTTTAGAAACTCTGTCACAGAACCCTGTATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTACGGCTTACGCTTAAACGTACCTTcc                                                                     >  1:174931/1‑149 (MQ=255)
                                                                                  cgcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCa                                                                    <  2:364958/149‑1 (MQ=255)
                                                                                  cgcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTGCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCa                                                                    <  1:928452/149‑1 (MQ=255)
                                                                                      gTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCtt                                                                <  1:70957/149‑1 (MQ=255)
                                                                                         ttAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGcc                                                             >  2:965225/1‑149 (MQ=255)
                                                                                           aGAATCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTTCATGGAGTAAAGCTGGACGATCAGGTGTTCGTTAATGTCCGGAGACAGATCAGAACGGGCCGGCTTACGCTTAAACGAACCTTCCATCTTgacgg                                                           >  2:569481/1‑145 (MQ=255)
                                                                                            gAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGc                                                          <  1:407237/149‑1 (MQ=255)
                                                                                                 tctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGTCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCaa                                                     <  2:552431/149‑1 (MQ=255)
                                                                                                       acaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCTCAGACAGATCAGAACGCTCCGGCTTACTCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcca                                               <  1:536356/149‑1 (MQ=255)
                                                                                                       acaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccc                                               >  2:775271/1‑148 (MQ=255)
                                                                                                       acaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcca                                               <  1:675303/149‑1 (MQ=255)
                                                                                                       acaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcca                                               >  1:29882/1‑149 (MQ=255)
                                                                                                         aGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccagc                                             >  1:731423/1‑149 (MQ=255)
                                                                                                           aaCCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccagcca                                           >  2:239199/1‑149 (MQ=255)
                                                                                                               cttaaTTGTGTCCTCTCTTTCGTACTATGCTTTACTTGGAGTAATGCTCGACTTTCAGGTGTTCGTTAATGTCCGCAGTCAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGtt                                       <  2:535096/145‑1 (MQ=255)
                                                                                                                 gctttGTGTCCTCTCTTTTGTACTAAGATTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTgg                                     <  1:680913/146‑1 (MQ=255)
                                                                                                                  cATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccagccaacctggc                                    <  1:615501/149‑8 (MQ=255)
                                                                                                                     tgtgTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCAGCAGACAGATCATAACGCTCCGGCTTACGCTTAAACGTACCTTACATCTTGCCAGCATCAACTTACAGCCAGGTTGGCttt                                 >  1:160116/1‑149 (MQ=255)
                                                                                                                       ggtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGGTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTc                               <  1:385433/148‑1 (MQ=255)
                                                                                                                           ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                           >  2:636189/1‑149 (MQ=255)
                                                                                                                              tctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctgc                        >  2:750866/1‑149 (MQ=255)
                                                                                                                                tttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTc                      >  2:143359/1‑149 (MQ=255)
                                                                                                                                tttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCGTACGCTTAAACGTTCCTTCCATCGTGCGAGCATCAACTTCCCGCCTGGTTGGCTTTTCACGCTGCTc                      >  2:852075/1‑149 (MQ=255)
                                                                                                                                  tGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcag                    <  2:849353/149‑1 (MQ=255)
                                                                                                                                     tACTAAGCTTTACTTGGAGTAAAGCTCGCCGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAACCGTACCTGCCATCTTGCCAGCAGCAACTTCCAGCCAGGTTGGCTTTGCACGCTGCTcagcca                 >  2:628916/1‑149 (MQ=255)
                                                                                                                                      aCTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccag                >  2:984029/1‑149 (MQ=255)
                                                                                                                                       tcgaGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTACGCAGAAAGATCAGAACGCTCCGGCTTACGCTTAAACGTACATTCCATCTTGCCAGCATCAACTTCCAGCCAGGTAGGCTTTTCACGCAGCTcagcccgc               >  2:737501/4‑149 (MQ=255)
                                                                                                                                       cTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccagc               <  1:520168/149‑1 (MQ=255)
                                                                                                                                          aGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTcc            >  1:1015525/1‑149 (MQ=255)
                                                                                                                                          aGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAAAGCTCCGGCTTACGCTTAAAAGTACCTTCCATCTTGCCAGCATCCACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTcc            >  1:256508/1‑149 (MQ=255)
                                                                                                                                             tttACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCagc         >  2:823586/1‑148 (MQ=255)
                                                                                                                                                 cTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCgg     <  1:960488/149‑1 (MQ=255)
                                                                                                                                                  ttGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCTGCCAGGTTTGCTTGTCACGCTGCTCCGCCAGCTCCAGCGCGGc    >  1:570919/1‑149 (MQ=255)
                                                                                                                                                  ttGGAGTAAAGCTCGACGATCAGGGGGGCGTTAATGTCCGCCGACAGATCAGAACGCTCCGGCTTACGCTTAAATGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGc    <  2:562907/149‑1 (MQ=255)
                                                                                                                                                    ggAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCtt  >  2:380686/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTT  >  NC_000913/3440913‑3441209

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: ‑