Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,355,859 G→A 7.0% G76G (GGC→GGT pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,355,8590GA 7.0% 964.3 / 41.3 317G76G (GGC→GGTpnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (170/123);  new base A (16/6);  total (186/131)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 321 total aligned reads displayed.

GGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCG  >  NC_002947/5355711‑5356006
                                                                                                                                                    |                                                                                                                                                   
ggtcggtcTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGc                                                                                                                                                    <  2:315165/150‑1 (MQ=255)
ggtcggtcTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGc                                                                                                                                                    <  1:2179456/150‑1 (MQ=255)
 gtcggtcTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGcc                                                                                                                                                   <  1:4762829/150‑1 (MQ=255)
     gtcTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAGGAAGCCACCt                                                                                                                                               >  1:1789147/1‑150 (MQ=255)
     gtcTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCt                                                                                                                                               >  1:5394515/1‑150 (MQ=255)
      tcttctTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTg                                                                                                                                              >  2:4511862/1‑150 (MQ=255)
        ttcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTggg                                                                                                                                            >  2:9069193/1‑150 (MQ=255)
           ttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATc                                                                                                                                         >  2:5000060/1‑150 (MQ=255)
             gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                                                       >  2:173041/1‑150 (MQ=255)
             gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                                                       >  1:332867/1‑150 (MQ=255)
               tggtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGATGCCACCTGGGATCTTGc                                                                                                                                     >  1:1340739/1‑150 (MQ=255)
               tggtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGc                                                                                                                                     >  1:201319/1‑150 (MQ=255)
               tggtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGc                                                                                                                                     <  1:2159879/150‑1 (MQ=255)
               tggtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGc                                                                                                                                     >  2:7798110/1‑150 (MQ=255)
                 gtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCa                                                                                                                                   <  1:1043690/150‑1 (MQ=255)
                  tggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAg                                                                                                                                  >  2:2580778/1‑150 (MQ=255)
                  tggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAg                                                                                                                                  >  2:8385192/1‑150 (MQ=255)
                   ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                                                                                                 <  1:8972674/150‑1 (MQ=255)
                    gAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcg                                                                                                                                <  1:9101525/150‑1 (MQ=255)
                         cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                                                                                                           <  2:3733265/150‑1 (MQ=255)
                         cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                                                                                                           >  2:419389/1‑150 (MQ=255)
                          aCGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAg                                                                                                                          >  2:8180622/1‑150 (MQ=255)
                          aCGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAg                                                                                                                          >  2:6892953/1‑150 (MQ=255)
                                  gacgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCCCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTc                                                                                                                  >  1:7489667/1‑150 (MQ=255)
                                            gtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGttg                                                                                                        <  2:4846608/147‑3 (MQ=255)
                                            aCTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGa                                                                                                        <  1:6974547/150‑1 (MQ=255)
                                            aCTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGa                                                                                                        <  1:6974754/150‑1 (MQ=255)
                                              ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                      >  2:1993037/1‑150 (MQ=255)
                                              ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                      >  1:6339921/1‑150 (MQ=255)
                                              ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                      <  1:8817973/150‑1 (MQ=255)
                                              ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                      >  1:3862112/1‑150 (MQ=255)
                                              ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTGTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                                                      >  2:7203946/1‑150 (MQ=255)
                                                 gTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGa                                                                                                   >  2:5321805/1‑150 (MQ=255)
                                                      tGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCg                                                                                              >  2:8468410/1‑150 (MQ=255)
                                                       gAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                                                                                             >  2:7585097/1‑150 (MQ=255)
                                                       gAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                                                                                             <  1:8704128/150‑1 (MQ=255)
                                                        aaGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCggg                                                                                            <  1:3133790/150‑1 (MQ=255)
                                                          gCCTTCCGGTAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaa                                                                                          <  2:5776999/150‑1 (MQ=255)
                                                          gCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaa                                                                                          <  2:612220/150‑1 (MQ=255)
                                                            cTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaaga                                                                                        >  1:268599/1‑150 (MQ=255)
                                                              tCCGCGAACAGCGGACGGATCGGACCGACGATCAGGCGCGACCACAGCGTCTCAATCTCCCAACCACGGCCTTCACGCTTGAAGAACCCACCTGGGATCTTGCCAGCGGCGTAGGTCATTTCCTGCTAGTGAACCGACAGCCCGAAGaaa                                                                                      <  2:5635174/150‑1 (MQ=255)
                                                                 gggAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAccc                                                                                   <  2:3789633/150‑1 (MQ=255)
                                                                 gggAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAccc                                                                                   <  2:3744368/150‑1 (MQ=255)
                                                                  ggAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCt                                                                                  >  2:6613888/1‑150 (MQ=255)
                                                                   gAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCtt                                                                                 <  1:7917106/150‑1 (MQ=255)
                                                                   gAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCtt                                                                                 >  2:6729345/1‑150 (MQ=255)
                                                                   gAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtg                                                                                 >  1:6340884/1‑149 (MQ=255)
                                                                    aaaaGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtga                                                                                <  2:1321170/147‑3 (MQ=255)
                                                                     aCAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGc                                                                               <  2:2740444/150‑1 (MQ=255)
                                                                     aCAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTACACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCATTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGc                                                                               <  2:298800/150‑1 (MQ=255)
                                                                      aaGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtgacc                                                                              >  2:7831803/2‑146 (MQ=255)
                                                                       aGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCt                                                                             >  1:512652/1‑150 (MQ=255)
                                                                        gCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTg                                                                            >  2:6137451/1‑150 (MQ=255)
                                                                         cggacggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTgg                                                                           >  1:583349/1‑150 (MQ=255)
                                                                               gaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGc                                                                     >  2:6332240/1‑150 (MQ=255)
                                                                                 tCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                                                    <  1:6507495/133‑1 (MQ=255)
                                                                                 tCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGCTCAGCCt                                                                   >  1:6139694/1‑150 (MQ=255)
                                                                                    gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                                                                <  1:8249627/150‑1 (MQ=255)
                                                                                     aCGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGttt                                                               <  2:296810/150‑1 (MQ=255)
                                                                                      cGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTg                                                              >  1:7347758/1‑150 (MQ=255)
                                                                                       ggTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTgg                                                             >  2:8827598/1‑150 (MQ=255)
                                                                                       ggTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTgg                                                             >  2:2520882/1‑150 (MQ=255)
                                                                                          cGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGgcgc                                                          <  1:3764917/150‑1 (MQ=255)
                                                                                             tCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTa                                                       <  1:7983753/150‑1 (MQ=255)
                                                                                             tCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTa                                                       >  1:3894511/1‑150 (MQ=255)
                                                                                             tCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTAGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTa                                                       >  1:4476215/1‑150 (MQ=255)
                                                                                                ggCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACca                                                    <  1:3487493/150‑1 (MQ=255)
                                                                                                 gcgcGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcac                                                   >  2:5193928/1‑150 (MQ=255)
                                                                                                  cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                                                  <  2:5139482/150‑1 (MQ=255)
                                                                                                  cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                                                  <  1:1162368/150‑1 (MQ=255)
                                                                                                         tCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCa                                           >  2:1200053/1‑150 (MQ=255)
                                                                                                               tctcTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACgg                                     >  2:7118954/1‑150 (MQ=255)
                                                                                                               tctcTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACgg                                     >  2:9132544/1‑150 (MQ=255)
                                                                                                               tctcTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACgg                                     >  1:4992296/1‑150 (MQ=255)
                                                                                                                ctctTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGt                                    <  2:2675911/150‑1 (MQ=255)
                                                                                                                   tttCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAc                                 >  2:8088584/1‑150 (MQ=255)
                                                                                                                      ctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGATACGCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtg                              >  2:2246122/1‑149 (MQ=255)
                                                                                                                      ctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtc                              <  1:6861865/150‑1 (MQ=255)
                                                                                                                      ctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtc                              <  1:6007823/150‑1 (MQ=255)
                                                                                                                       tcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtcg                             >  1:1595441/1‑150 (MQ=255)
                                                                                                                       tcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTAGCGCCTACCACAGTCACCAGCACGGTGAcgtcg                             >  1:6432841/1‑150 (MQ=255)
                                                                                                                         ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                           >  2:2434471/1‑150 (MQ=255)
                                                                                                                         ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                           >  1:8379973/1‑150 (MQ=255)
                                                                                                                         ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                           >  1:3575900/1‑150 (MQ=255)
                                                                                                                         ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                           >  1:5706675/1‑150 (MQ=255)
                                                                                                                         ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt                           <  2:7347758/150‑1 (MQ=255)
                                                                                                                          gAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTg                          >  1:783177/1‑150 (MQ=255)
                                                                                                                              gACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGa                      >  2:7425504/1‑150 (MQ=255)
                                                                                                                                   ccTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTa                 >  1:1933316/1‑150 (MQ=255)
                                                                                                                                    cTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGcc                                                                              <  2:7275076/88‑1 (MQ=255)
                                                                                                                                    cTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTaa                <  1:6058311/150‑1 (MQ=255)
                                                                                                                                    cTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTaa                >  2:6504465/1‑150 (MQ=255)
                                                                                                                                     ttCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAAc               >  2:5626253/1‑150 (MQ=255)
                                                                                                                                       cACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCa             >  2:5410850/1‑150 (MQ=255)
                                                                                                                                          gCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCa          >  1:1738041/1‑150 (MQ=255)
                                                                                                                                           cTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCAc         <  1:6605763/150‑1 (MQ=255)
                                                                                                                                              gaagaaGCCAACTGGGATCTTGCCAGCGGCGTAGGACTTTTCCTGGTAGTTAACCGACAGCGGGAAGAAACCCTAGCCTGGATCAGCCAGGATGGCGCCTACCACAGACACCAGCACGGTGAAGACCTAGACGACGGGAACCAGCACGgc      >  2:6035750/1‑150 (MQ=255)
                                                                                                                                                agaagCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGgcgc    >  2:3009906/1‑150 (MQ=255)
                                                                                                                                                agaagCCACCTGGGATCTTGCCAGCGGCGTAGGTCCTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGgcgc    <  2:3629477/150‑1 (MQ=255)
                                                                                                                                                 gaagCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGcc   >  1:294896/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTTGTCGACGGTAACCAGCACGGCGCCG  >  NC_002947/5355711‑5356006

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: