Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,682,899 C→A 6.2% I6I (ATC→ATA rstB → sensory histidine kinase of RstAB two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,682,8990CA 6.2% 230.6 / 2.8 76I6I (ATC→ATArstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base A (3/2);  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 = 1.00e+00

CTGGATCGTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGT  >  NC_000913/1682750‑1683046
                                                                                                                                                     |                                                                                                                                                   
cTGGATCGTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATc                                                                                                                                                     >  2:65753/1‑150 (MQ=255)
    aTCGTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCTTTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGt                                                                                                                                                 <  1:1346783/150‑1 (MQ=255)
    aTCGTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGt                                                                                                                                                 >  2:298670/1‑150 (MQ=255)
     tagTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGATCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGtt                                                                                                                                                <  1:490188/148‑1 (MQ=255)
      cGTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGttt                                                                                                                                               <  2:641556/150‑1 (MQ=255)
       gTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGtttt                                                                                                                                              <  1:589697/150‑1 (MQ=255)
         aGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGGATAAACAGtttttt                                                                                                                                            <  2:1208963/150‑3 (MQ=255)
           cGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCt                                                                                                                                          <  2:60414/150‑1 (MQ=255)
                aCGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAt                                                                                                                                     >  2:183994/1‑150 (MQ=255)
                 cGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAtt                                                                                                                                    >  2:352980/1‑150 (MQ=255)
                  gtagCTATTTCTCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAATACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAttg                                                                                                                                   <  2:62025/147‑1 (MQ=255)
                   tGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAttgt                                                                                                                                  >  2:181471/1‑150 (MQ=255)
                   gggCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAttgt                                                                                                                                  <  2:903671/149‑1 (MQ=255)
                    ggCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAttgtt                                                                                                                                 <  2:391807/150‑1 (MQ=255)
                     gCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAttgttt                                                                                                                                >  1:10267/1‑150 (MQ=255)
                      cTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTAttgtttg                                                                                                                               <  1:760650/150‑1 (MQ=255)
                        aTTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTc                                                                                                                             >  1:577650/1‑150 (MQ=255)
                        aTTTCGCGGTTAAGAAAAAAACTGCTCGATAACGACGCAGTACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATACATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTc                                                                                                                             <  1:1199407/150‑1 (MQ=255)
                          ttagcgGGTAAGAAAAAAACTGCTCGATAACGCCGCAGAAGCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGACTCATGCATGGGAATAAGCGATGAAAAAACTGTTAATCCAGTTTTACCTGTTATTGTTTGTCTg                                                                                                                           <  2:978036/147‑1 (MQ=255)
                                    aGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTtgtg                                                                                                                 >  2:1210309/1‑150 (MQ=255)
                                     gAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGa                                                                                                                >  2:670248/1‑150 (MQ=255)
                                     gAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGa                                                                                                                >  2:402042/1‑150 (MQ=255)
                                      aaaaaaaCTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGAt                                                                                                               <  1:65753/150‑1 (MQ=255)
                                            aCTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGGTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTctct                                                                                                         >  2:839446/1‑150 (MQ=255)
                                             cTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTctg                                                                                                        >  2:464853/1‑150 (MQ=255)
                                                  cGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCGATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGt                                                                                                   >  2:246564/1‑150 (MQ=255)
                                                    aTAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTg                                                                                                 >  2:51011/1‑150 (MQ=255)
                                                    aTAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTg                                                                                                 <  1:922572/150‑1 (MQ=255)
                                                      aaCGCCGCAGAACCTTATCGCATAAAAACTGTGGGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTggg                                                                                               <  1:1150834/150‑1 (MQ=255)
                                                          ccgcAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTgg                                                                                           <  1:744268/150‑1 (MQ=255)
                                                          ccgcAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTgg                                                                                           >  2:160043/1‑150 (MQ=255)
                                                             cAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGgtgt                                                                                        >  1:399700/1‑150 (MQ=255)
                                                                  ccTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACaaa                                                                                   >  1:958746/1‑150 (MQ=255)
                                                                  ccTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGATGGTGTACaaa                                                                                   >  2:248333/1‑150 (MQ=255)
                                                                    ttATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAAtt                                                                                 <  1:670248/150‑1 (MQ=255)
                                                                         gCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTAccg                                                                            >  1:1064495/1‑150 (MQ=255)
                                                                          cATTAAAACTGTGCGTAACAAAGGATATCTTTTTGAGCCTCATGAATGGTAATAAGAGATGAAAAAAATGTTTATCAAGTTTTACATTTTATTGTTTGTATGATTCCTTTTGATGTCTCTGATTTTTGTGCTGGTGTACAAATTTAccgc                                                                           >  2:213803/1‑150 (MQ=255)
                                                                           aTTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTAccgcc                                                                          >  2:331372/1‑150 (MQ=255)
                                                                               aaaCTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAAc                                                                      >  1:501083/1‑150 (MQ=255)
                                                                                    gagCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGGTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGGCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAAcgcgcg                                                                 <  2:908700/148‑1 (MQ=255)
                                                                                      gCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCggg                                                               <  2:784341/150‑1 (MQ=255)
                                                                                       cGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGATTCTCCAGTTTTACCTGTTTTTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACAGCCGAACGCGCGGGc                                                              >  1:281422/1‑150 (MQ=255)
                                                                                          aaCAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCaaa                                                           >  2:872376/1‑150 (MQ=255)
                                                                                            cAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACa                                                         <  1:402042/150‑1 (MQ=255)
                                                                                            acaataatATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACa                                                         >  2:439869/8‑150 (MQ=255)
                                                                                               aGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTc                                                      <  1:464853/150‑1 (MQ=255)
                                                                                                  cTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCt                                                   <  1:1187498/150‑1 (MQ=255)
                                                                                                    aTCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTgg                                                 <  1:181471/150‑1 (MQ=255)
                                                                                                      cTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGgat                                               >  2:933082/1‑150 (MQ=255)
                                                                                                         taagcgcCTCAAGCCTGGGAATACGCGATGAAAAAACTGCTAATCCAGTTTTACCTGTTATTGCTTGGCTGCTCCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTAACCGCCGAACGCGCGGGCAAACAGTCGCTGgatgat                                            <  2:1717616/147‑1 (MQ=255)
                                                                                                             cgcCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTtga                                        <  1:439869/150‑1 (MQ=255)
                                                                                                              gcCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTtgat                                       <  1:160043/150‑1 (MQ=255)
                                                                                                                cTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTtgatga                                     <  1:1243210/150‑1 (MQ=255)
                                                                                                                cTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTtgatga                                     <  2:399700/150‑1 (MQ=255)
                                                                                                                  catgcatgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAAc                                   >  2:1417888/1‑150 (MQ=255)
                                                                                                                  catgcatgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAAc                                   >  1:1702183/1‑150 (MQ=255)
                                                                                                                  catgcatgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAAc                                   <  1:352980/150‑1 (MQ=255)
                                                                                                                    tgcatgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACTCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAg                                 >  1:204900/1‑150 (MQ=255)
                                                                                                                      catgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGtt                               >  2:1530618/1‑150 (MQ=255)
                                                                                                                       atgTGAATAAGCTATGAAAAAACTTTTTATACAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTTTGATGTCTCTGCTGGTTGGTCTGGTTTACAAATTTACCGCCGAACGCTCTGGCAAACAGTCGCTGGATGATTTGATGAACAGTTc                              <  1:213803/150‑1 (MQ=255)
                                                                                                                        tgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAATCAGTCGCTGGATGATTTGATGAACAGTTCg                             >  1:1585625/1‑150 (MQ=255)
                                                                                                                        tgGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCg                             <  2:10267/150‑1 (MQ=255)
                                                                                                                          ggAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTCCAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCt                           >  1:1059802/1‑150 (MQ=255)
                                                                                                                           gAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTg                          <  1:298670/150‑1 (MQ=255)
                                                                                                                                  cGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGa                   >  1:868887/1‑150 (MQ=255)
                                                                                                                                  cGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGa                   >  1:1251855/1‑150 (MQ=255)
                                                                                                                                   gATGAAAAAACTCTTTAACCAGTTTTAACTGTTAAACTTTGAGAGCAACCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTAACCCCGAACGCGCGCCCAAACAGACGCTGCATGAATAGATGAACAGATCGCAGTATCTGAt                  >  2:588528/1‑150 (MQ=255)
                                                                                                                                    caggtaaaaCTGGATAAACAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATg                 >  1:1688512/6‑150 (MQ=255)
                                                                                                                                    aTGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATg                 >  1:398025/1‑150 (MQ=255)
                                                                                                                                     tGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATgc                >  1:98956/1‑150 (MQ=255)
                                                                                                                                      ggtaaaaCTGGATAAACAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATgcg               >  2:21274/4‑150 (MQ=255)
                                                                                                                                      ggtaaaaCTGGATAAACAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATgcg               >  2:71268/4‑150 (MQ=255)
                                                                                                                                            aCTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGa         <  1:51011/150‑1 (MQ=255)
                                                                                                                                                tttATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTATAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTg     <  1:999906/150‑1 (MQ=255)
                                                                                                                                                tttATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTg     <  1:331372/150‑1 (MQ=255)
                                                                                                                                                 ttATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGc    <  1:1193646/150‑1 (MQ=255)
                                                                                                                                                   aTCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGt  >  2:1464355/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGGATCGTAGCGTGGACGTGGCTATTTCGCGGTTAAGAAAAAAACTGCTCGATAACGCCGCAGAACCTTATCGCATTAAAACTGTGCGTAACAAAGGCTATCTTTTTGCGCCTCATGCATGGGAATAAGCGATGAAAAAACTGTTTATCCAGTTTTACCTGTTATTGTTTGTCTGCTTCCTTGTGATGTCTCTGCTGGTTGGGCTGGTGTACAAATTTACCGCCGAACGCGCGGGCAAACAGTCGCTGGATGATTTGATGAACAGTTCGCTGTATCTGATGCGCAGCGAATTGCGT  >  NC_000913/1682750‑1683046

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