Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,923,376 A→G 5.1% V102A (GTG→GCG) ‡ rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,3760AG 5.1% 500.5 / 10.2 151V102A (GTG→GCG) ‡rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base A (78/65);  new base G (2/6);  total (80/71)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

Only 100 of 160 total aligned reads displayed.

GCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATG  >  NC_000913/3923228‑3923524
                                                                                                                                                    |                                                                                                                                                    
gCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAacgca                                                                                                                                                      >  2:177958/1‑149 (MQ=255)
      gCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                <  4:807327/149‑1 (MQ=255)
      gCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                >  1:772472/1‑149 (MQ=255)
      gCTCACCATATCGTTCAGAGAGGCAAAAGCGCGACTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                >  3:538930/1‑149 (MQ=255)
       cTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGttt                                                                                                                                               <  2:509587/149‑1 (MQ=255)
         cACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttcc                                                                                                                                             <  1:828741/149‑3 (MQ=255)
          aCCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTAc                                                                                                                                             >  4:223826/1‑148 (MQ=255)
            cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttccttc                                                                                                                                          <  3:632812/149‑6 (MQ=255)
            cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttccttc                                                                                                                                          <  3:418443/149‑6 (MQ=255)
            cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCaa                                                                                                                                          <  4:632394/149‑1 (MQ=255)
            cATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCaa                                                                                                                                          >  1:979829/1‑149 (MQ=255)
             atatCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAg                                                                                                                                         <  3:948033/149‑1 (MQ=255)
               atCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGc                                                                                                                                       >  3:985124/1‑149 (MQ=255)
                   ttCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATc                                                                                                                                   <  4:1002865/149‑1 (MQ=255)
                      agagagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCaa                                                                                                                                >  1:606020/1‑149 (MQ=255)
                      agagagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCaa                                                                                                                                >  4:451066/1‑149 (MQ=255)
                        agagGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACa                                                                                                                              >  1:1035595/1‑149 (MQ=255)
                          agGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACaga                                                                                                                            <  3:271780/149‑1 (MQ=255)
                            gCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCGTGCGtttcc                                                                                                                                             >  3:222743/1‑128 (MQ=255)
                            gCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGt                                                                                                                          >  1:440739/1‑149 (MQ=255)
                            gCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACCCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGt                                                                                                                          <  2:815812/149‑1 (MQ=255)
                              aaaaGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGa                                                                                                                        >  2:316843/1‑149 (MQ=255)
                                aaGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaaa                                                                                                                      <  4:902663/149‑1 (MQ=255)
                                aaGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaaa                                                                                                                      <  2:864753/149‑1 (MQ=255)
                                       gCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  1:177958/149‑1 (MQ=255)
                                       gCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  2:812221/149‑1 (MQ=255)
                                        cTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTc                                                                                                              >  4:46109/1‑149 (MQ=255)
                                         tAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCa                                                                                                             >  4:171573/1‑149 (MQ=255)
                                          aaTTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAg                                                                                                            <  3:1042616/149‑1 (MQ=255)
                                            ttACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGa                                                                                                          <  3:46109/149‑1 (MQ=255)
                                                gCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCa                                                                                                      >  2:930760/1‑149 (MQ=255)
                                                   aTCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGa                                                                                                   >  1:680638/1‑149 (MQ=255)
                                                   aTCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGa                                                                                                   >  3:17303/1‑149 (MQ=255)
                                                    tCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGAt                                                                                                  >  4:762166/1‑149 (MQ=255)
                                                       aaTGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCCGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATaga                                                                                               <  4:8288/149‑1 (MQ=255)
                                                           gcggcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGt                                                                                           >  2:270890/1‑149 (MQ=255)
                                                            cggcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTg                                                                                          <  1:516148/149‑1 (MQ=255)
                                                              gcTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGa                                                                                        >  2:324182/1‑149 (MQ=255)
                                                                tctGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAg                                                                                                                                         <  1:101800/98‑1 (MQ=255)
                                                                tctGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAg                                                                                                                                         >  2:101800/1‑98 (MQ=255)
                                                                   gAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGcc                                                                                   <  2:440739/149‑1 (MQ=255)
                                                                    aaGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCt                                                                                  >  3:41382/1‑149 (MQ=255)
                                                                     attccaCTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTg                                                                                 >  2:276324/6‑149 (MQ=255)
                                                                      ggAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTgg                                                                                >  1:403035/1‑149 (MQ=255)
                                                                       gAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGc                                                                               >  1:376472/1‑149 (MQ=255)
                                                                            tcttctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                          <  1:617391/149‑1 (MQ=255)
                                                                             cttctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCt                                                                         <  2:430366/149‑1 (MQ=255)
                                                                              ttctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTg                                                                        >  2:619603/1‑149 (MQ=255)
                                                                               tctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTgg                                                                       <  1:295307/149‑1 (MQ=255)
                                                                                ctACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGt                                                                      >  2:993966/1‑149 (MQ=255)
                                                                                  aCCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTcc                                                                    >  3:553601/1‑149 (MQ=255)
                                                                                  aCCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTcc                                                                    <  1:784918/149‑1 (MQ=255)
                                                                                  aCCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTcc                                                                    >  2:521262/1‑149 (MQ=255)
                                                                                   cccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCg                                                                   >  2:881470/1‑149 (MQ=255)
                                                                                   cccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCg                                                                   >  3:751136/1‑149 (MQ=255)
                                                                                    ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGa                                                                                                   <  4:683941/116‑1 (MQ=255)
                                                                                    ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCgg                                                                  >  2:594539/1‑149 (MQ=255)
                                                                                    ccTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCgg                                                                  >  3:23698/1‑149 (MQ=255)
                                                                                     cTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGt                                                                 >  1:344875/1‑149 (MQ=255)
                                                                                     cTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGt                                                                 <  2:505917/149‑1 (MQ=255)
                                                                                      tGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGtt                                                                                                                                                <  4:583171/69‑1 (MQ=255)
                                                                                      tGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCATGGTCCGGTg                                                                >  2:349170/1‑149 (MQ=255)
                                                                                         tctGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCg                                                             >  1:939148/1‑149 (MQ=255)
                                                                                              acaggTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCGGACATc                                                        <  2:368701/146‑1 (MQ=255)
                                                                                                  gTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATa                                                    >  2:380226/1‑149 (MQ=255)
                                                                                                   ttCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATaa                                                   >  2:701613/1‑149 (MQ=255)
                                                                                                    tCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCt                                                                                  >  3:811066/1‑117 (MQ=255)
                                                                                                      aaTATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAAcc                                                >  3:576212/1‑149 (MQ=255)
                                                                                                       atatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCg                                               >  3:72873/1‑149 (MQ=255)
                                                                                                        tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGt                                              >  3:457459/1‑149 (MQ=255)
                                                                                                          ttCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTc                                            <  3:371374/149‑1 (MQ=255)
                                                                                                            ctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCAc                                          <  3:421456/149‑1 (MQ=255)
                                                                                                             tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATagaga                                                                                             <  4:160858/97‑1 (MQ=255)
                                                                                                                 gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCg                                                                   <  2:499613/119‑1 (MQ=255)
                                                                                                                 gTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGc                                     >  3:670239/1‑149 (MQ=255)
                                                                                                                       gCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATAc                               <  2:1035595/149‑1 (MQ=255)
                                                                                                                          cATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGt                            <  1:881617/149‑1 (MQ=255)
                                                                                                                             gTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               >  4:100168/1‑63 (MQ=255)
                                                                                                                             gTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCtt                                                                                                               <  3:100168/63‑1 (MQ=255)
                                                                                                                              ttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATc                                                        <  3:858093/117‑1 (MQ=255)
                                                                                                                              ttGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATc                                                        >  4:858093/1‑117 (MQ=255)
                                                                                                                                  aCCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCcaccac                    >  3:73190/1‑149 (MQ=255)
                                                                                                                                  aCCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCcaccac                    >  3:292698/1‑149 (MQ=255)
                                                                                                                                    cTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACaa                  >  2:676197/1‑149 (MQ=255)
                                                                                                                                     tGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAAt                 <  4:73190/149‑1 (MQ=255)
                                                                                                                                         ttggtaaaCGCGTGCGTTTACCAAGGCTATTCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCtt                                       >  2:614687/6‑123 (MQ=255)
                                                                                                                                         gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGaa                                                                                                                       <  1:529/43‑1 (MQ=255)
                                                                                                                                         gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGc                                                               <  2:981986/99‑1 (MQ=255)
                                                                                                                                         gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGc                                                               >  1:981986/1‑99 (MQ=255)
                                                                                                                                         gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTa             >  4:376583/1‑149 (MQ=255)
                                                                                                                                         gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTa             >  4:360144/1‑149 (MQ=255)
                                                                                                                                           aGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATc           >  4:804091/1‑149 (MQ=255)
                                                                                                                                            ggAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCg          >  3:363932/1‑149 (MQ=255)
                                                                                                                                             gAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCga         <  3:168073/149‑1 (MQ=255)
                                                                                                                                                acgcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGaa      <  2:403035/149‑1 (MQ=255)
                                                                                                                                                  gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCa                                    <  4:489958/117‑1 (MQ=255)
                                                                                                                                                  gcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAata    >  3:802904/1‑149 (MQ=255)
                                                                                                                                                    acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                          >  2:46790/1‑77 (MQ=255)
                                                                                                                                                    acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTcc                                                                          <  1:46790/77‑1 (MQ=255)
                                                                                                                                                    acgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATg  >  1:384777/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACCAGCTCACCATATCGTTCAGAGAGGCAAAAGCGCGGCTAATTACGCCATCAAATGGCGGCTCTGAAGGAAACTCTTCTACCCTGCTCTGTACTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATG  >  NC_000913/3923228‑3923524

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