Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,442,594 G→T 6.1% L124L (CTG→CTT)  PP_3060 → tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,442,5940GT 6.1% 618.4 / 20.3 180L124L (CTG→CTT) PP_3060tail sheath protein
Reads supporting (aligned to +/- strand):  ref base G (79/90);  new base T (4/7);  total (83/97)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

GGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTCCC  >  NC_002947/3442445‑3442743
                                                                                                                                                     |                                                                                                                                                     
ggCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTg                                                                                                                                                       >  2:1789658/1‑150 (MQ=255)
     ggCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCa                                                                                                                                                  >  1:3270333/1‑150 (MQ=255)
       cTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACg                                                                                                                                                <  2:3258183/150‑1 (MQ=255)
        tGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACgg                                                                                                                                               <  1:736184/150‑1 (MQ=255)
        tGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACgg                                                                                                                                               <  1:468392/150‑1 (MQ=255)
        tGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACgg                                                                                                                                               >  1:2479478/1‑150 (MQ=255)
         gACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGc                                                                                                                                              <  1:343299/150‑1 (MQ=255)
            ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                           <  1:127191/150‑1 (MQ=255)
            ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                           <  2:2975884/150‑1 (MQ=255)
             gCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCTTg                                                                                                                                          <  2:629530/150‑1 (MQ=255)
               aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                        >  1:4413784/1‑150 (MQ=255)
               aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                        <  2:1981215/150‑1 (MQ=255)
               aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                        >  2:4287072/1‑150 (MQ=255)
               aaCGCCTTCAACGGCGTATTCCGAGGCTCGGTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                        <  2:3548369/150‑1 (MQ=255)
                    cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAc                                                                                                                                   <  1:2931295/150‑1 (MQ=255)
                      tCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACct                                                                                                                                 <  2:4512208/150‑1 (MQ=255)
                        aaCGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctac                                                                                                                               <  1:1430645/150‑1 (MQ=255)
                           ggCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCtactac                                                                                                                            <  2:410225/150‑1 (MQ=255)
                                aTTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCt                                                                                                                       <  1:1008313/150‑1 (MQ=255)
                                 ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                      <  1:801025/150‑1 (MQ=255)
                                 ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                      <  1:3546598/150‑1 (MQ=255)
                                 ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                      <  1:3126810/150‑1 (MQ=255)
                                 ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                      <  1:2767545/150‑1 (MQ=255)
                                 ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACAAGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGcgc                                                                                                                      >  2:3513403/1‑148 (MQ=255)
                                   ccGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGa                                                                                                                    <  2:4181332/150‑1 (MQ=255)
                                    cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                   >  2:4169035/1‑150 (MQ=255)
                                    cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                   <  2:238098/150‑1 (MQ=255)
                                    cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                   >  1:2022338/1‑150 (MQ=255)
                                    cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                   <  1:1205210/150‑1 (MQ=255)
                                    cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                   <  1:3374467/150‑1 (MQ=255)
                                    cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                   <  1:528410/150‑1 (MQ=255)
                                      aGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGt                                                                                                                 >  1:2119588/1‑150 (MQ=255)
                                          tCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggc                                                                                                             >  1:4353405/1‑150 (MQ=255)
                                                aaGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAg                                                                                                       <  2:3963049/150‑1 (MQ=255)
                                                    tCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCa                                                                                                   >  1:3446014/1‑150 (MQ=255)
                                                       aGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCgg                                                                                                                                               >  1:2739366/1‑103 (MQ=255)
                                                        gAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTa                                                                                               >  1:1966674/1‑150 (MQ=255)
                                                           gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                            <  2:3018997/150‑1 (MQ=255)
                                                           gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                            <  2:4485473/150‑1 (MQ=255)
                                                           gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                            <  1:1107977/150‑1 (MQ=255)
                                                               gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAATCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                                        <  1:2425759/150‑1 (MQ=255)
                                                               gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                                        <  1:1407320/150‑1 (MQ=255)
                                                               gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCATGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                                        <  2:2675259/150‑1 (MQ=255)
                                                                    ccTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGAccc                                                                                   <  1:1750785/150‑1 (MQ=255)
                                                                      ttGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCgg                                                                                 >  1:3991701/1‑150 (MQ=255)
                                                                         cGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCg                                                                              >  1:1946024/1‑150 (MQ=255)
                                                                            tcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                           >  2:1920631/1‑150 (MQ=255)
                                                                            tcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                           <  2:2479478/150‑1 (MQ=255)
                                                                            tcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                           >  1:3441057/1‑150 (MQ=255)
                                                                                 gtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAg                                                                      >  2:527109/1‑150 (MQ=255)
                                                                                    aCCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgt                                                                   <  1:784998/150‑1 (MQ=255)
                                                                                    aCCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgt                                                                   >  2:829856/1‑150 (MQ=255)
                                                                                    aCCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgt                                                                   >  1:2338804/1‑150 (MQ=255)
                                                                                      ccTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGc                                                                                                                                                   <  2:2780437/68‑1 (MQ=255)
                                                                                      ccTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgtgt                                                                 <  2:1282389/150‑1 (MQ=255)
                                                                                            cgGCACCCAGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAATCTAGAAGTCGGAGGCGAGGTCATATACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCaa                                                           <  1:922038/150‑1 (MQ=255)
                                                                                             ggCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAAc                                                          >  1:2056475/1‑150 (MQ=255)
                                                                                                aCCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGc                                                       >  2:1690335/1‑150 (MQ=255)
                                                                                                 cccTGGAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCa                                                      >  2:1441367/1‑150 (MQ=255)
                                                                                                  ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATTAAGCGTGTCATCAACGGCAc                                                     <  1:2792115/150‑1 (MQ=255)
                                                                                                  ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                                     <  2:3377015/150‑1 (MQ=255)
                                                                                                     gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACTGa                                                  >  1:150735/1‑150 (MQ=255)
                                                                                                     gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACTGa                                                  >  1:3701962/1‑150 (MQ=255)
                                                                                                     gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGa                                                  <  2:2123474/150‑1 (MQ=255)
                                                                                                             tGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTgg                                          >  1:3810943/1‑150 (MQ=255)
                                                                                                             tGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTgg                                          >  1:3194182/1‑150 (MQ=255)
                                                                                                              ggACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGc                                         <  2:2876627/150‑1 (MQ=255)
                                                                                                               gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                        >  1:1493128/1‑150 (MQ=255)
                                                                                                               gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                        <  1:1772209/150‑1 (MQ=255)
                                                                                                                aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                       <  1:3977054/150‑1 (MQ=255)
                                                                                                                aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                       <  2:3775825/150‑1 (MQ=255)
                                                                                                                aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                       <  1:2689569/150‑1 (MQ=255)
                                                                                                                    gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                                   >  1:1533461/1‑150 (MQ=255)
                                                                                                                    gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                                   >  2:1733599/1‑150 (MQ=255)
                                                                                                                      gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgc                                 <  1:2936607/150‑1 (MQ=255)
                                                                                                                          cTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgcgcgc                             <  1:3824312/150‑1 (MQ=255)
                                                                                                                           tGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGcc                            <  2:1236700/150‑1 (MQ=255)
                                                                                                                             gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGa                          >  2:4289671/1‑150 (MQ=255)
                                                                                                                              aaGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCCCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGAc                         <  1:1789658/150‑1 (MQ=255)
                                                                                                                                 gCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTc                      <  2:1716837/150‑1 (MQ=255)
                                                                                                                                   cGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCgg                    >  2:2966000/1‑150 (MQ=255)
                                                                                                                                    ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                   <  1:642988/150‑1 (MQ=255)
                                                                                                                                    ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                   <  1:722130/150‑1 (MQ=255)
                                                                                                                                    ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                   <  2:1127949/150‑1 (MQ=255)
                                                                                                                                       gACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg                >  1:2183234/1‑150 (MQ=255)
                                                                                                                                        aCAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCCGTCAGCGCGCCGACCTCGGCCTga               >  2:4392108/1‑149 (MQ=255)
                                                                                                                                        aCAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt               >  2:4392244/1‑150 (MQ=255)
                                                                                                                                         cAATGCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTa              >  1:3911868/1‑150 (MQ=255)
                                                                                                                                          aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa             <  1:3042061/150‑1 (MQ=255)
                                                                                                                                          aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa             >  2:2360596/1‑150 (MQ=255)
                                                                                                                                            tgCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc           >  2:1953662/2‑150 (MQ=255)
                                                                                                                                             gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccc          >  1:1604712/1‑150 (MQ=255)
                                                                                                                                               cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc        <  2:4380065/150‑1 (MQ=255)
                                                                                                                                               cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc        >  1:3590329/1‑150 (MQ=255)
                                                                                                                                                  gCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCttt     <  1:4169035/150‑1 (MQ=255)
                                                                                                                                                   cTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTc    >  2:369901/1‑150 (MQ=255)
                                                                                                                                                    tGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGAAACCCCCTATcc   >  1:3602163/1‑150 (MQ=255)
                                                                                                                                                     gAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTccc  <  2:4113721/150‑1 (MQ=255)
                                                                                                                                                     gAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTccc  <  2:4113667/150‑1 (MQ=255)
                                                                                                                                                     gAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTccc  >  2:2815820/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTCCC  >  NC_002947/3442445‑3442743

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