Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,508,221 (G)5→4 coding (1441/1962 nt) rlhA → 23S rRNA 5‑hydroxycytidine C2501 synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,508,2170G.98.2% 230.4 / ‑1.8 55coding (1437/1962 nt)rlhA23S rRNA 5‑hydroxycytidine C2501 synthase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (29/25);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

TGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCGGTTGCTGATCCTGCGCCGGTTTATCCG  >  NC_000913/1508079‑1508352
                                                                                                                                          |                                                                                                                                       
tgtCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCgg                                                                                                                                        >  1:225516/1‑140 (MQ=255)
tgtCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCgg                                                                                                                                        >  1:752871/1‑140 (MQ=255)
 gtCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggg                                                                                                                                       >  2:649359/1‑140 (MQ=255)
   cAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggggc                                                                                                                                     >  1:80900/1‑139 (MQ=255)
    aGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggggcg                                                                                                                                    <  1:649359/140‑3 (MQ=255)
    aGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggggcg                                                                                                                                    <  2:807213/140‑3 (MQ=255)
      cATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggggcgct                                                                                                                                  <  1:1077619/140‑5 (MQ=255)
        tCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGt                                                                                                                               <  1:266452/140‑1 (MQ=255)
        tCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGt                                                                                                                               >  2:161699/1‑140 (MQ=255)
          aCGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGttt                                                                                                                             <  1:175348/140‑1 (MQ=255)
          aCGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGttt                                                                                                                             >  1:737209/1‑140 (MQ=255)
            gCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGt                                                                                                                           <  2:225516/140‑1 (MQ=255)
                aCGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCa                                                                                                                       >  1:743727/1‑140 (MQ=255)
                  gCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCaaa                                                                                                                     <  1:574418/140‑1 (MQ=255)
                  gCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCaaa                                                                                                                     >  2:384942/1‑140 (MQ=255)
                      gACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGt                                                                                                                 >  2:770225/1‑140 (MQ=255)
                          ggCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCTAACAGTCTGt                                                                                                             >  1:794860/1‑140 (MQ=255)
                          ggCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGt                                                                                                             >  2:205264/1‑140 (MQ=255)
                                tCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAAcc                                                                                                       >  1:140353/1‑140 (MQ=255)
                                tCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAAcc                                                                                                       >  2:940865/1‑140 (MQ=255)
                                tCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAAcc                                                                                                       >  1:1069895/1‑140 (MQ=255)
                                         ccAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccgcc                                                                                              >  2:126924/1‑140 (MQ=255)
                                         ccAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccgcc                                                                                              >  1:601899/1‑140 (MQ=255)
                                         ccAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGAAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccgcc                                                                                              >  2:758318/1‑140 (MQ=255)
                                             tAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGa                                                                                          <  1:940865/140‑1 (MQ=255)
                                                 gCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAgct                                                                                      >  2:9635/1‑140 (MQ=255)
                                                  ccGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAgctg                                                                                     >  2:513849/1‑140 (MQ=255)
                                                   cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAgctgc                                                                                    >  1:218928/1‑140 (MQ=255)
                                                   cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAgc                                                                                       <  2:666011/137‑1 (MQ=255)
                                                   cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAgc                                                                                       >  1:666011/1‑137 (MQ=255)
                                                   cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAACCAGTTCCGCCGTGAAgctgct                                                                                   <  1:161699/140‑1 (MQ=255)
                                                    gAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAgctgct                                                                                   <  1:384942/140‑1 (MQ=255)
                                                      aaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGa                                                                                 <  1:770225/140‑1 (MQ=255)
                                                           cAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGc                                                                            >  1:199271/1‑140 (MQ=255)
                                                                      cTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTgcg                                                                 <  2:80900/140‑1 (MQ=255)
                                                                      cTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTgcg                                                                 <  2:737209/140‑1 (MQ=255)
                                                                          aGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTc                                                             >  2:996709/1‑140 (MQ=255)
                                                                           ggATGGTCTGGCAAAACTGGGGCAACCCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCt                                                            <  1:126924/140‑1 (MQ=255)
                                                                           ggATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCt                                                            <  2:752871/140‑1 (MQ=255)
                                                                              tGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGc                                                         <  2:199271/140‑1 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTgcg                                                                 <  2:422945/131‑1 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTgcg                                                                 >  1:422945/1‑131 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGcgcg                                                               <  1:104316/133‑1 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGcgcg                                                               >  2:104316/1‑133 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGcc                                                        >  1:222977/1‑140 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGcc                                                        >  2:566789/1‑140 (MQ=255)
                                                                               ggTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGcc                                                        >  2:721563/1‑140 (MQ=255)
                                                                                          aCTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGtt                                                    <  2:203819/133‑1 (MQ=255)
                                                                                          aCTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGtt                                                    >  1:203819/1‑133 (MQ=255)
                                                                                             ggggCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCgg                                          <  1:377233/140‑1 (MQ=255)
                                                                                             ggggCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCgg                                          <  2:1069895/140‑1 (MQ=255)
                                                                                             ggggCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAATTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCgg                                          <  2:802735/140‑1 (MQ=255)
                                                                                                     cccTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTa                                  >  1:285914/1‑140 (MQ=255)
                                                                                                        tctATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAAc                               >  2:112314/1‑140 (MQ=255)
                                                                                                         ctATTACGCCGGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAGCAGTCTGTTAAGCCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAAcc                              <  2:1155304/140‑1 (MQ=255)
                                                                                                         ctATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAAcc                              <  2:236495/140‑1 (MQ=255)
                                                                                                          tATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCg                             <  1:566789/140‑1 (MQ=255)
                                                                                                                    gcgATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCGGTTGCTGATc                   <  1:170110/140‑1 (MQ=255)
                                                                                                                     cgATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAAc                               <  2:745799/127‑1 (MQ=255)
                                                                                                                     cgATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAAc                               >  1:745799/1‑127 (MQ=255)
                                                                                                                      gATGTGCAAATTAATTTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCGGTTGCTGATCCt                 >  2:529277/1‑140 (MQ=255)
                                                                                                                                ttaattTGCCGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCGGTTGCTGATCCTGCGCCGGttt       <  1:205264/140‑1 (MQ=255)
                                                                                                                                      tGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCGGTTGCTGATCCTGCGCCGGTTTATCCg  >  2:656848/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTGAAGCTGCTGACATGCTGGATGCTGCGCGTCTTGCCAGTTACCAGCGCGGCAGCCGTAAACCGGTTGCTGATCCTGCGCCGGTTTATCCG  >  NC_000913/1508079‑1508352

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