Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,965,940 T→G 10.2% intergenic (‑128/+136) fliE ← / ← atoC flagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,965,9400TG10.2% 825.0 / 6.6 313intergenic (‑128/+136)fliE/atoCflagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC
Reads supporting (aligned to +/- strand):  ref base T (87/172);  new base G (38/2);  total (139/174)
Fisher's exact test for biased strand distribution p-value = 3.11e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

GAATTCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGG  >  NC_002947/4965792‑4966088
                                                                                                                                                    |                                                                                                                                                    
gAATTCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCACCAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCtt                                                                                                                                                     <  1:6965709/150‑1 (MQ=255)
  aTTCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGc                                                                                                                                                   >  2:265612/1‑150 (MQ=255)
  aTTCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGc                                                                                                                                                   >  2:4674346/1‑150 (MQ=255)
   ttCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCg                                                                                                                                                  >  1:199207/1‑150 (MQ=255)
   ttCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCg                                                                                                                                                  >  1:3437540/1‑150 (MQ=255)
      aaCACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCATCTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTcg                                                                                                                                               >  2:5892302/1‑150 (MQ=255)
      aaCACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTcg                                                                                                                                               >  2:3710381/1‑150 (MQ=255)
      aaCACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCCTGCgaag                                                                                                                                               >  1:4514838/1‑147 (MQ=255)
           cTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTcgc                                                                                                                                              >  1:7458961/1‑146 (MQ=255)
            ttGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGCGCGATAg                                                                                                                                         >  1:953641/1‑150 (MQ=255)
            ttGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                         >  2:7766151/1‑150 (MQ=255)
                cTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                     <  1:2095998/150‑1 (MQ=255)
                 tCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGc                                                                                                                                                   <  2:8805308/135‑1 (MQ=255)
                   aTGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCa                                                                                                                                  >  2:5812068/1‑150 (MQ=255)
                      aCCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCagag                                                                                                                               <  1:8611670/150‑1 (MQ=255)
                         cTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCCATATGGCCGCAGAGCgg                                                                                                                            >  1:6675914/1‑150 (MQ=255)
                          ttCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGc                                                                                                                           >  2:3148432/1‑150 (MQ=255)
                               tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCgg                                                                                                                      >  2:6888677/1‑150 (MQ=255)
                               tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGGCGCCGAGCGGcccccg                                                                                                                      >  1:1884234/1‑148 (MQ=255)
                               tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGAGAGGGCCGCCGAGCGGCCCCgg                                                                                                                      >  1:569906/1‑150 (MQ=255)
                                ctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTCTGCGGccctatc                                                                                                                     <  2:3437540/150‑5 (MQ=255)
                                ctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccctatc                                                                                                                     <  1:3610643/150‑5 (MQ=255)
                                 tcCTGCGGCCGCATTTTTTTGACAGTCAGGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccctatcg                                                                                                                    <  1:300269/150‑6 (MQ=255)
                                    tGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGGCGCGATAGGGCCGCACAGCGGCCCCCGCGAta                                                                                                                 >  2:4013636/1‑149 (MQ=255)
                                     gCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCt                                                                                                                <  1:1105731/150‑1 (MQ=255)
                                        gCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGCGCGATTGGGCCCCAGAGCGGCCCCCCCGATCGTTg                                                                                                             >  1:1072639/1‑150 (MQ=255)
                                         ccGCATTTTTTTGACCGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgc                                                                                                            <  2:1624000/150‑1 (MQ=255)
                                         ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCCGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgc                                                                                                            >  2:1066388/1‑150 (MQ=255)
                                          cGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGCGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGCTCTTTgcg                                                                                                           >  2:6704566/1‑150 (MQ=255)
                                           gCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcgg                                                                                                          <  1:507520/150‑1 (MQ=255)
                                           gCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcgg                                                                                                          >  1:7162191/1‑150 (MQ=255)
                                             aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                        <  2:8696520/150‑1 (MQ=255)
                                             aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                        >  1:8027610/1‑150 (MQ=255)
                                               ttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCGTTTCGGCGGa                                                                                                      >  2:4559280/1‑150 (MQ=255)
                                                tttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTGGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTTCGGCGGAg                                                                                                     >  2:3513168/1‑150 (MQ=255)
                                                   ttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTg                                                                                                  <  1:911478/150‑1 (MQ=255)
                                                      aCAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCCGCGCGGCCCCGGCGCGCGTTGCGGCGGCGCGGAAg                                                                                               >  1:4590000/1‑150 (MQ=255)
                                                            aTGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgg                                                                                         <  2:3601455/150‑1 (MQ=255)
                                                             tGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTggg                                                                                        <  1:336900/150‑1 (MQ=255)
                                                               cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAt                                                                                                                  <  1:1699834/122‑1 (MQ=255)
                                                               cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGc                                                                                      <  1:4177273/150‑1 (MQ=255)
                                                               cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGc                                                                                      <  1:2853191/150‑1 (MQ=255)
                                                                 gcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCCTGGGGCGCCGGAGGGCCCCGAAGCGGCCCCCGCGATGGGTTGGGCGGGGGAGAAGACCTGGGGGCg                                                                                    >  1:266297/1‑150 (MQ=255)
                                                                           tacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGccc                                                                          <  1:4102374/150‑1 (MQ=255)
                                                                            acctaAAGCAACAAAGGGGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGcccc                                                                         <  1:2398817/150‑1 (MQ=255)
                                                                              ctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGGGGCGCGGTAGGGCCGCCCAGCCGCCCCCCCCCGCTTTTGGGCGGGGCGGGAGACCAGGGGGGGCGGCGCGccccca                                                                       >  1:3175090/1‑148 (MQ=255)
                                                                              ctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGCGCGCTAGGGGCGCCCGGCGGCCCCCGCGGTCGTTGGGGCGGAGGTGGAGAACAGGGGGCGCCTCGCGcccccg                                                                       >  1:8364538/1‑148 (MQ=255)
                                                                                aaaGCAACACAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTc                                                                     <  1:5280230/150‑1 (MQ=255)
                                                                                aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTc                                                                     <  2:797757/150‑1 (MQ=255)
                                                                                aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTc                                                                     <  2:7868489/150‑1 (MQ=255)
                                                                                 aaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTcg                                                                    <  2:616799/150‑1 (MQ=255)
                                                                                 aaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATAGGGGCGCAGAGCGGCCCCCGCGAACGTTGCGGCGGGGCTGAAGACCAGGGGCCGCGTCGcccccccaacg                                                                    >  1:1341467/1‑145 (MQ=255)
                                                                                   gCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGATAGGGGCGCAGCGCGGCCCCCGCGGGCTTTGCGGCGGGGCGGTAGACCAGGGGCCGCGTCGCGCCCCaatagcc                                                                  >  2:3942423/1‑144 (MQ=255)
                                                                                    caacaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGa                                                                 <  1:6842016/150‑1 (MQ=255)
                                                                                    caacaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGa                                                                 <  2:418899/150‑1 (MQ=255)
                                                                                    caacaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCGATCTTTGCGGCGGGGCTGAAGAACTGGGGCCCCGGCGCGCCCCAGGcgcgc                                                                 >  2:2537544/1‑149 (MQ=38)
                                                                                        aaaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCt                                                                                                                                                      >  1:4957974/1‑61 (MQ=255)
                                                                                        aaaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCt                                                                                                                                                      <  2:4957974/61‑1 (MQ=255)
                                                                                         caGGTGCCAGCCTTCGCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTGGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                            <  2:7000261/149‑1 (MQ=255)
                                                                                         aaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                            <  2:2311583/150‑1 (MQ=255)
                                                                                         aaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCGATCTTTGCGGCGGGGCTGAAGACCAGGGGCCGCTTCGCGCCCCCATcgcgccccca                                                            >  2:8106427/1‑144 (MQ=38)
                                                                                          aGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAg                                                           <  2:1224969/150‑1 (MQ=255)
                                                                                          aGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGCTAGGGCCCCAGAGCGGCCCCCGCGCGCTGTGCGGCGGGGCTGGAGACCCGGGGCCGCGGCGCGCCCCCATCCCGACACAAg                                                           >  1:6029176/1‑150 (MQ=255)
                                                                                           ggTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccc                                                                                                                         >  2:2103383/1‑87 (MQ=38)
                                                                                             tGCCAGCCTTCGCAACTCCTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGcc                                                        <  2:1511623/150‑1 (MQ=255)
                                                                                                cAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGtt                                                     <  2:5525484/150‑1 (MQ=255)
                                                                                                cAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCACAGCGGCCCCGGCGAACGTTGCGGCGGAGCGGGAGACCAGGGGCCGCGTCGCGCCCCAATCGCGACACAAGGCCggt                                                     >  2:6613477/1‑148 (MQ=255)
                                                                                                 aGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGCTAGGGCCGCAGAGCGGCCCCCGCGATCTTTGCGGCGGAGCGGGAGACCAGGGGCCGCGTCGCGCCCCAGTCGCGACACAAGGCCGTTc                                                    >  2:473202/1‑150 (MQ=39)
                                                                                                 aGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTc                                                    <  2:199207/150‑1 (MQ=255)
                                                                                                  gCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTcc                                                   <  2:3315857/150‑1 (MQ=255)
                                                                                                    cTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                 <  2:2268445/150‑1 (MQ=255)
                                                                                                    cTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCCGCCCCCGCGATCTTTTCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                 >  1:5862490/1‑150 (MQ=255)
                                                                                                      tCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCcaca                                               <  1:1147640/150‑1 (MQ=255)
                                                                                                        gCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                             <  2:8339297/150‑1 (MQ=255)
                                                                                                          aaCTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATTCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCGAACTTTTCGGCGGAGCTGAAGAACAGGGGCCGCTTCGCGCCCCATTAGCGACACAAGGCCGTTACCACAGAgg                                           >  2:3283074/1‑150 (MQ=255)
                                                                                                            cTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTc                                         <  2:3017971/150‑1 (MQ=255)
                                                                                                               aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                      <  2:7162191/150‑1 (MQ=255)
                                                                                                               aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                      <  2:2839047/150‑1 (MQ=255)
                                                                                                                 ttCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCgcgc                                    <  1:4300984/150‑1 (MQ=255)
                                                                                                                   gacAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCa                                  <  2:7553674/148‑1 (MQ=255)
                                                                                                                   cccAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCa                                  <  2:426186/150‑1 (MQ=255)
                                                                                                                   cccAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCa                                  <  1:3710381/150‑1 (MQ=255)
                                                                                                                    ccAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAGAGCCGCCCCCGCGAGCGTTGCGGCGGCGCTGGAGACCGGGGGCCGCGTCGCGCCCCCGTCGCGCCACAACGCCGGTCCCACAGAGGGCCGCGcccg                                 >  2:1884195/1‑148 (MQ=255)
                                                                                                                      agagATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGtt                               <  2:1171657/150‑1 (MQ=255)
                                                                                                                         gaTCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCa                            <  1:3284930/150‑1 (MQ=255)
                                                                                                                           tCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATTGGGCCGCAGTGCGGCCCCGGCGTTCTTTGCGGCGGTGCTGAAGACCTGGGGCCGCGTCGCGCCCCCGTCGCGACACAAGGCCGTTCCCACCGTGGGCCGCGCCCGTTGCAga                          >  2:562032/1‑150 (MQ=37)
                                                                                                                           tCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCCGAGCGGCCCCCGCGCGCGTTGCGGCGGCGCGGAAGACCCGGGGCCGCGTCGCGCCCCCCGCGCCACCCCACGCCCGTCCCACACAGGTCCGCGCCCGGTGccgc                          >  2:5161297/1‑147 (MQ=255)
                                                                                                                            cacacaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAgat                         <  2:1738115/150‑1 (MQ=255)
                                                                                                                               acaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCGGCGCGCCCCAGTCGCCACACAAGGCCGTTCCCACACAGGTCCGCCCCAGTTTCAcatgat                      >  2:7086146/1‑150 (MQ=37)
                                                                                                                                     ccTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGCGCTGAAGACCAGGGGCCGCGTCGCGCCCCAATCGCGACACCAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCa                >  2:8618059/1‑150 (MQ=38)
                                                                                                                                         cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCa            <  1:100368/150‑1 (MQ=255)
                                                                                                                                            gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGcc         >  2:5408811/1‑150 (MQ=255)
                                                                                                                                            gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGcc         >  1:550927/1‑150 (MQ=255)
                                                                                                                                            gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAACGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGcc         >  1:6982434/1‑150 (MQ=255)
                                                                                                                                                ggCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAAc     <  1:4674346/150‑1 (MQ=255)
                                                                                                                                                 gCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACa    >  1:5051355/1‑150 (MQ=255)
                                                                                                                                                 gCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACa    <  1:788878/150‑1 (MQ=255)
                                                                                                                                                 gCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACa    <  1:6483505/150‑1 (MQ=255)
                                                                                                                                                  ccTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCCGCCGACCAg   >  1:6260266/1‑150 (MQ=255)
                                                                                                                                                   cTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAgg  >  2:4378233/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAATTCAACACCTTGGCTCATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGG  >  NC_002947/4965792‑4966088

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: