Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,957,744 C→G 5.8% G334A (GGC→GCC)  fliG ← flagellar motor switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,957,7440CG 5.8% 214.8 / 2.7 69G334A (GGC→GCC) fliGflagellar motor switch protein
Reads supporting (aligned to +/- strand):  ref base C (38/27);  new base G (2/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TCCCTTGTAGGAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGTTCCGAAGCACGTTTG  >  NC_002947/4957599‑4957887
                                                                                                                                                 |                                                                                                                                               
tCCCTTGTAGGAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCtt                                                                                                                                             >  1:697169/1‑150 (MQ=255)
   cTTGTAGGAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCc                                                                                                                                          >  2:18481‑M2/1‑149 (MQ=255)
      gTAGGAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCcgcc                                                                                                                                       >  1:358968‑M2/1‑146 (MQ=255)
        aGGAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCcgccga                                                                                                                                     >  2:825329‑M2/1‑144 (MQ=255)
          gAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGc                                                                                                                                   >  1:903770/1‑150 (MQ=255)
          gAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGc                                                                                                                                   >  1:529835/1‑150 (MQ=255)
              gggTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGa                                                                                                                               >  1:896226/1‑150 (MQ=255)
                            aTGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACGTTAGAGCATTGCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCg                                                                                                                 >  2:1040668/1‑150 (MQ=255)
                                   ggCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATgc                                                                                                          >  1:632603/1‑150 (MQ=255)
                                         gACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCggg                                                                                                    <  1:371304/150‑1 (MQ=255)
                                         gACAACCTCGATTGCCAGGACCGACGCCTTCGAGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGAACACGTGCCGCCGAGCACGATCTAGCCGGCCTGGGCCATgcgttggg                                                                                                    >  2:459755/1‑145 (MQ=255)
                                             aCCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGGCAAGGGcgccgaggaaatgatctaa                                                                                                                        <  1:15494‑M2/126‑20 (MQ=255)
                                             aCCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGGCAAGGGcgccgaggaaatgatctaa                                                                                                                        >  2:15494‑M2/1‑107 (MQ=255)
                                               cTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATgg                                                                                              >  2:756672/1‑150 (MQ=255)
                                               cTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATgg                                                                                              >  1:965206/1‑150 (MQ=255)
                                                tCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGt                                                                                             <  2:632603/150‑1 (MQ=255)
                                                tCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGt                                                                                             >  1:1084809/1‑150 (MQ=255)
                                                  gATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGa                                                                                           <  1:685072/150‑1 (MQ=255)
                                                       ccAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGAtt                                                                                      >  2:182178/1‑150 (MQ=255)
                                                        cAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGAttt                                                                                     >  2:358086/1‑150 (MQ=255)
                                                             ccGACGCATTCGCGGGTGAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGAtttctttc                                                                                <  1:1145255/150‑1 (MQ=255)
                                                              cGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGa                                                                                        <  1:605348/141‑1 (MQ=255)
                                                              cGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGa                                                                                        >  2:605348/1‑141 (MQ=255)
                                                                     ttCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGt                                                                        >  2:678806/1‑150 (MQ=255)
                                                                     ttCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATAGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGt                                                                        >  2:275012/1‑150 (MQ=255)
                                                                      tCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGtt                                                                       <  2:530664/150‑1 (MQ=255)
                                                                         cgGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTcc                                                                    <  2:1084809/150‑1 (MQ=255)
                                                                           ggTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGACTAGCTTAGATCATTTCCTCGGCTCCCTTGCCGCCGAGCAGGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCAc                                                                  >  1:812857/1‑150 (MQ=255)
                                                                           ggTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCAc                                                                  <  1:127178/150‑1 (MQ=255)
                                                                                aCCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCa                                                                   >  2:465025/1‑144 (MQ=255)
                                                                                aCCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCa                                                                   <  1:465025/144‑1 (MQ=255)
                                                                                 cccGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCt                                                            <  1:600343/150‑1 (MQ=255)
                                                                                 cccGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCt                                                            >  2:868087/1‑150 (MQ=255)
                                                                                  ccGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGGCAAGGGcgccgaggaaatgatctaag                                                                                                                       <  1:840386‑M2/90‑21 (MQ=255)
                                                                                  ccGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGGCAAGGGcgccgaggaaatgatctaag                                                                                                                       >  2:840386‑M2/1‑70 (MQ=255)
                                                                                    gCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGAc                                                         <  2:903770/150‑1 (MQ=255)
                                                                                         tACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCAGGCGCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCAGGCCGCTGACACGTa                                                    >  1:386010/1‑150 (MQ=255)
                                                                                         tACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGAtttct                                                                                   <  2:473283/119‑1 (MQ=255)
                                                                                         tACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGAtttct                                                                                   >  1:473283/1‑119 (MQ=255)
                                                                                            aaaaGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATg                                                                                               >  1:155887/1‑104 (MQ=255)
                                                                                            aaaaGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATg                                                                                               <  2:155887/104‑1 (MQ=255)
                                                                                              aaGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAAt                                                                                                <  2:23589/101‑1 (MQ=255)
                                                                                              aaGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAAt                                                                                                >  1:23589/1‑101 (MQ=255)
                                                                                                    gcgTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTg                                         <  1:1057037/150‑1 (MQ=255)
                                                                                                    gcgTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAAACGTTAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTg                                         <  1:30117/150‑1 (MQ=255)
                                                                                                           cTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCa                                  >  2:166199/1‑150 (MQ=255)
                                                                                                             aTGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGTCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGa                                >  2:662014/1‑150 (MQ=255)
                                                                                                               gggATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATc                              <  2:318243/150‑1 (MQ=255)
                                                                                                               gggATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATc                              <  1:358086/150‑1 (MQ=255)
                                                                                                                       gAATACCTTAGATCATTTCCTCGGCGCACTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTc                           <  2:1144045/145‑1 (MQ=255)
                                                                                                                       gAATACCTTAGATCATTTCCTCGGCGCACTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTc                           >  1:1144045/1‑145 (MQ=255)
                                                                                                                          tACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACgcagca                   <  2:980907/150‑1 (MQ=255)
                                                                                                                            ccTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAAt                                                                                                <  1:534274/71‑1 (MQ=255)
                                                                                                                            ccTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAAt                                                                                                >  2:534274/1‑71 (MQ=255)
                                                                                                                             cTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCAc                         <  2:1027905/141‑1 (MQ=255)
                                                                                                                             cTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCAc                         >  1:1027905/1‑141 (MQ=255)
                                                                                                                             cTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGtt                >  1:638095/1‑150 (MQ=255)
                                                                                                                             cTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGtt                >  1:332131/1‑150 (MQ=255)
                                                                                                                             cTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGtt                >  2:403617/1‑150 (MQ=255)
                                                                                                                                  aTCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTgcgc                                                                           >  1:436524/1‑86 (MQ=255)
                                                                                                                                  aTCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTgcgc                                                                           <  2:436524/86‑1 (MQ=255)
                                                                                                                                     aTTTCCTCGGCGCCCTTGTCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGTTCCGAAGCa        <  1:662014/150‑1 (MQ=255)
                                                                                                                                      tttCCTCGGCGCCCTTGCcgccgagcaaatgatctaaggtattcgcgatcccatagccgacgccgccttttgtaggagcgggttt                                                                        <  1:900640‑M2/85‑68 (MQ=255)
                                                                                                                                      tttCCTCGGCGCCCTTGCcgccgagcaaatgatctaaggtattcgcgatcccatagccgacgccgccttttgtaggagcgggttt                                                                        >  2:900640‑M2/1‑18 (MQ=255)
                                                                                                                                         ccTCGGCGCCCTTGCcgccgaggaaatgatctaaggtattcgcgatcccatagccgacgccgccttttgtaggagcgggtttacccgcgaaggcgtcggtcctggcaatcgaggttgtcgctgccgacgcattcgcgggtgaacccgctc    <  2:675983‑M2/150‑136 (MQ=255)
                                                                                                                                           tCGGCGCCCTTGCcgccgaggaaatgatctaaggtattcgcgatcccatagccgacgccgccttttgta                                                                                   <  1:1064837‑M2/69‑57 (MQ=255)
                                                                                                                                           tCGGCGCCCTTGCcgccgaggaaatgatctaaggtattcgcgatcccatagccgacgccgccttttgta                                                                                   >  2:1064837‑M2/1‑13 (MQ=255)
                                                                                                                                           tCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGTTCCGAAGCACGTTTg  >  1:1082753/1‑150 (MQ=255)
                                                                                                                                           tCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGTTCCGAAGCACGTTTg  >  1:1088051/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCCCTTGTAGGAGCGGGTTCACCCGCGAATGCGTCGGCAGCGACAACCTCGATTGCCAGGACCGACGCCTTCGCGGGTAAACCCGCTCCTACAAAAGGCGGCGTCGGCTATGGGATCGCGAATACCTTAGATCATTTCCTCGGCGCCCTTGCCGCCGAGCACGATCTCGCCGGCCTCGGCCATGCGCCGGGCAATGGTGAGGATTTCTTTCTGCGCGGTTTCCACGTCGCTGACCCGTACCGGCCCTTTGGCCTCCAGATCGTCACGCAGCAGTTCCGAAGCACGTTTG  >  NC_002947/4957599‑4957887

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