Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,923,719 G→C 5.6% intergenic (‑39/+25) rsmG ← / ← mnmG 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,7190GC 5.6% 1413.1 / 57.9 391intergenic (‑39/+25)rsmG/mnmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base G (168/201);  new base C (7/15);  total (175/216)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 399 total aligned reads displayed.

CGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTT  >  NC_000913/3923570‑3923868
                                                                                                                                                     |                                                                                                                                                     
      gTTCCATTTATGCAGCCTATTCCCGTCGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTa                                                                                                                                                 <  2:3858626/150‑1 (MQ=255)
       ttCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTAc                                                                                                                                                <  2:1806465/150‑1 (MQ=255)
         ccATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCa                                                                                                                                              <  1:185041/150‑1 (MQ=255)
         ccATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCa                                                                                                                                              <  1:256203/150‑1 (MQ=255)
         ccATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCa                                                                                                                                              <  1:1289445/150‑1 (MQ=255)
           aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCAgg                                                                                                                                            <  2:4218679/150‑1 (MQ=255)
               aTGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAtt                                                                                                                                        <  1:558701/150‑1 (MQ=255)
                   agcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTAGCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTTCCAGGCATTTTTa                                                                                                                                    >  1:1774085/1‑150 (MQ=255)
                   agcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTa                                                                                                                                    <  1:5649776/150‑1 (MQ=255)
                   agcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTa                                                                                                                                    <  1:3381854/150‑1 (MQ=255)
                    gcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTaa                                                                                                                                   >  1:3047970/1‑150 (MQ=255)
                         ttCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGt                                                                                                                              <  1:2453236/150‑1 (MQ=255)
                          tCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGtt                                                                                                                             <  1:1627308/150‑1 (MQ=255)
                          tCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGtt                                                                                                                             >  2:3960276/1‑150 (MQ=255)
                           cACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTa                                                                                                                            <  2:2628732/150‑1 (MQ=255)
                            aCGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTAt                                                                                                                           >  1:4661850/1‑150 (MQ=255)
                             cGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGTAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATg                                                                                                                          <  2:1285553/150‑1 (MQ=255)
                             cGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATg                                                                                                                          <  1:1708661/150‑1 (MQ=255)
                              gTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGACATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATgc                                                                                                                         <  2:3562011/150‑1 (MQ=255)
                              gTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATgc                                                                                                                         <  2:3655369/150‑1 (MQ=255)
                               tAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATgcg                                                                                                                        >  2:5596716/1‑150 (MQ=255)
                               tAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATgcg                                                                                                                        <  2:3718394/150‑1 (MQ=255)
                               tAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAAAGCGTTATgcg                                                                                                                        >  1:3189764/1‑150 (MQ=255)
                                    aaTAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTacg                                                                                                                   <  1:2023665/150‑1 (MQ=255)
                                       aaGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTacgacg                                                                                                                <  2:4013618/150‑1 (MQ=255)
                                        aGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgc                                                                                                               <  1:2346265/150‑1 (MQ=255)
                                        aGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgc                                                                                                               <  2:722070/150‑1 (MQ=255)
                                         gCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgca                                                                                                              <  2:1489447/150‑1 (MQ=255)
                                         gCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgca                                                                                                              <  1:5894089/150‑1 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGACATACCTGCGTCGTTCAGCAGTCAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTGACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             <  2:5407324/150‑1 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             <  1:4791493/150‑1 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             >  2:929610/1‑150 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             >  2:614954/1‑150 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             >  1:2741253/1‑150 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             <  1:3867161/150‑1 (MQ=255)
                                          cTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAATTTTAATGCGTTATGCGCTACGACgcag                                                                                                             >  1:38122/1‑150 (MQ=255)
                                             gTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCAt                                                                                                          >  2:3512973/1‑150 (MQ=255)
                                             gTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCAt                                                                                                          >  2:5130455/1‑150 (MQ=255)
                                             gTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTAAGACGCAGCAt                                                                                                          >  1:536617/1‑150 (MQ=255)
                                              tttttCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATa                                                                                                         >  1:4433989/1‑150 (MQ=255)
                                                    tGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGt                                                                                                   <  1:1020139/150‑1 (MQ=255)
                                                         aTCTGTAATCGAAAAAACTGAGTCATACAGCACTAAGGAGAGTGTGTTGAGAAAGGATAATACAAGATCTTGATGCGTTGCCTGGTAAGAGGGTGCTTACCAGGCATTTTTAATGCGTTATGCTCTACGACGCAGCAAACCAGGtttatt                                                                                              >  2:4308061/1‑150 (MQ=255)
                                                         aTCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtttttt                                                                                              >  2:1604946/1‑150 (MQ=255)
                                                         aTCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGtttttt                                                                                              >  1:4288291/1‑150 (MQ=255)
                                                            ggTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCCTTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAg                                                                                           >  1:3620216/1‑150 (MQ=255)
                                                            ggTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAg                                                                                           <  2:4651188/150‑1 (MQ=255)
                                                             gTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGc                                                                                          <  2:498463/150‑1 (MQ=255)
                                                             gTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGc                                                                                          >  2:5292706/1‑150 (MQ=255)
                                                             gTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGc                                                                                          <  2:4150909/150‑1 (MQ=255)
                                                                    aaaTACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACca                                                                                   <  2:5819041/150‑1 (MQ=255)
                                                                    aaaTACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACca                                                                                   >  1:5012573/1‑150 (MQ=255)
                                                                            gCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATg                                                                           >  2:4794169/1‑150 (MQ=255)
                                                                             cGTCTGTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATgg                                                                          <  2:1691712/150‑1 (MQ=255)
                                                                                cTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGgaga                                                                       <  2:4455435/150‑1 (MQ=255)
                                                                                   tCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATgg                                                                    <  1:5298364/150‑1 (MQ=255)
                                                                                   tCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTAACAGGCATTTTTAATTCGTTATTCGCTACGACGCAGCATACCCAGTTTTTTCAGCCACACCAGCAGAATGGAGATgg                                                                    >  1:4907268/1‑150 (MQ=255)
                                                                                    cagcagTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTTTGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGc                                                                   <  1:3542276/150‑1 (MQ=255)
                                                                                       cagTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGc                                                                >  2:5212734/1‑150 (MQ=255)
                                                                                          tAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAgg                                                             <  1:2941753/150‑1 (MQ=255)
                                                                                          tAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAgg                                                             >  2:2458517/1‑150 (MQ=255)
                                                                                           aaGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGc                                                            <  2:4734821/150‑1 (MQ=255)
                                                                                             ggAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGt                                                          <  2:5452734/150‑1 (MQ=255)
                                                                                             ggAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGt                                                          <  1:3308543/150‑1 (MQ=255)
                                                                                             ggAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGt                                                          <  1:2775012/150‑1 (MQ=255)
                                                                                              gagaGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGGCTGGTAAGCGGGTGCTTTACAGGGATTTTTAATTGGTTATTCGCTACGTAGCAGCATAGCCTGTTTTTTAATGCACACCAGCAGAATGGTGATTGGCGCAGGCGTg                                                         >  1:2860948/1‑150 (MQ=255)
                                                                                              gagaGTTTGTTGAGCACGGTGATAACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTg                                                         >  2:3460625/1‑150 (MQ=255)
                                                                                                  gTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGc                                                     <  1:3985644/150‑1 (MQ=255)
                                                                                                   tttGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGcc                                                    <  1:4389645/150‑1 (MQ=255)
                                                                                                    ttgttgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCa                                                   <  1:1464036/150‑1 (MQ=255)
                                                                                                     tgttgAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAg                                                  >  2:4918199/1‑150 (MQ=255)
                                                                                                       ttgAGCACGGTGATTACCGGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGaa                                                >  2:2814262/1‑150 (MQ=255)
                                                                                                           gCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTCCTGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAGGGCCGCAGGCGTGACGCCAGAAATAc                                            <  1:5050038/150‑1 (MQ=255)
                                                                                                           gCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAATTTTAATTCGTTATTCGCTAAGACGCAGCATACCCAGTTTTTTCAGCCACACCACCACAAAGGAGATGGCCGCACGCGTGACGCCACAAAAAc                                            >  1:623678/1‑150 (MQ=255)
                                                                                                             aCGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATAcgc                                          <  2:87581/150‑1 (MQ=255)
                                                                                                                     tACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTg                                  <  1:3666020/150‑1 (MQ=255)
                                                                                                                       ccTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGc                                <  2:1922375/150‑1 (MQ=255)
                                                                                                                          gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTGATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGAGTGACGCCAGAAATACGCGAAGCTTGGCCGa                             >  2:1997688/1‑150 (MQ=255)
                                                                                                                              ttGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATaga                         <  2:3315521/150‑1 (MQ=255)
                                                                                                                               tGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGTAATACGCGAAGCTTGGCCGATagag                        <  2:5778713/150‑1 (MQ=255)
                                                                                                                               tGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATagag                        <  2:1616179/150‑1 (MQ=255)
                                                                                                                                gATGCGTTGCCTGGTTCGCGGGTGCTTACCAGGCATTTTTACTGCGTTATGCGCTACGACGCAGCATACCCTGCTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAgg                       <  1:5202242/150‑1 (MQ=255)
                                                                                                                                        gCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtg               <  2:135746/150‑1 (MQ=255)
                                                                                                                                          cgggTAAGCGGGTGCTTACCAGGCATTTTTCATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg             <  2:3103921/148‑1 (MQ=255)
                                                                                                                                          cTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg             <  1:858298/150‑1 (MQ=255)
                                                                                                                                          cTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg             >  1:5912823/1‑150 (MQ=255)
                                                                                                                                          cTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgtg             >  1:4090537/1‑150 (MQ=255)
                                                                                                                                           tGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGa            <  1:4512244/150‑1 (MQ=255)
                                                                                                                                           tGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGAGAGAGGCTGGTTTGTGa            >  1:3151325/1‑150 (MQ=255)
                                                                                                                                            ggTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGAt           >  1:2255114/1‑150 (MQ=255)
                                                                                                                                            ggTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGAt           >  1:2156002/1‑150 (MQ=255)
                                                                                                                                            ggTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGAt           >  1:5513236/1‑150 (MQ=255)
                                                                                                                                            ggTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGAt           >  1:198010/1‑150 (MQ=255)
                                                                                                                                             gTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATc          >  1:5725499/1‑150 (MQ=255)
                                                                                                                                             gTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATc          <  2:2375192/150‑1 (MQ=255)
                                                                                                                                              tAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATCGAGGCTGGTTTGTGAGCg         >  2:151139/1‑150 (MQ=255)
                                                                                                                                              tAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCg         >  1:1173727/1‑150 (MQ=255)
                                                                                                                                              tAAGCGGGTGCTTACCAGGCAATTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCg         <  2:2989380/150‑1 (MQ=255)
                                                                                                                                                aGCGTGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATTGAGATGTCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTTAtattt       <  2:2076295/150‑5 (MQ=255)
                                                                                                                                                 gCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTa      >  1:2659361/1‑150 (MQ=255)
                                                                                                                                                   gggTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAg    >  2:1106740/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTT  >  NC_000913/3923570‑3923868

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