Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,364 G→A 5.3% intergenic (+81/+45) groL → / ← pykA chaperonin GroEL/pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,551,3640GA 5.3% 344.9 / 6.6 115intergenic (+81/+45)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (52/54);  new base A (2/4);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 123 total aligned reads displayed.

CCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCACCCAGT  >  NC_002947/1551215‑1551503
                                                                                                                                                     |                                                                                                                                               
ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCg                                                                                                                                                 <  2:736085/150‑1 (MQ=255)
  ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCCTGATGTAAGCAAGAGAACCTGCCCTGCACCATCAAAACACACCGGTTCGCGAGAGCCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                                                               >  1:522621/1‑150 (MQ=255)
  ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                                                               <  1:2207887/150‑1 (MQ=255)
  ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                                                               >  1:2208376/1‑150 (MQ=255)
   gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCCGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgcag                                                                                                                                              <  2:1254319/150‑1 (MQ=255)
      gACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGcac                                                                                                                                           <  2:1454351/150‑1 (MQ=255)
         aaGCCAGCTGCTGACATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATAATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACaac                                                                                                                                        >  2:1485562/1‑149 (MQ=255)
          aGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAgg                                                                                                                                       <  2:2299991/150‑1 (MQ=255)
            ccAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCt                                                                                                                                     >  1:2177743/1‑150 (MQ=255)
                   ctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGTCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTg                                                                                                                          <  1:368659/150‑1 (MQ=255)
                    tgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGa                                                                                                                         <  1:1920563/150‑1 (MQ=255)
                    tgGCATGCCTGATATGGCTGGCATGTGTGGCATCGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGCTCGCGAGAACCGGGGCTTTTTTATTGCGGTCGGCGCCACCTCGCAGCACAAGCGT‑‑‑GCGCttta                                                                                                                         <  2:272664/150‑3 (MQ=255)
                         tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCCGTGCGGTGGCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCtt                                                                                                                    <  2:1230381/150‑1 (MQ=255)
                         tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCtt                                                                                                                    <  2:900705/150‑1 (MQ=255)
                           tggtggagcaaggtgactgggtgatcctgaccaagggtgacagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCACAGCACAAGGCT‑‑‑GCGCtt                                                                                                                           <  2:1322878‑M1/22‑1 (MQ=255)
                           tggtggagcaaggtgactgggtgatcctgaccaagggtgacagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCACAGCACAAGGCT‑‑‑GCGCtt                                                                                                                           >  1:1322878‑M1/120‑141 (MQ=255)
                           ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                                                               >  2:2166071/1‑125 (MQ=255)
                           ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                                                               <  1:2166071/125‑1 (MQ=255)
                             tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGc                                                                                                                >  1:1218654/1‑150 (MQ=255)
                                   caaggtgactgggtgatcctgaccaagggtgacagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCACAGCACAAGGCT‑‑‑GCGCTTGAGCCTtgt                                                                                                                  >  2:2099233‑M1/112‑142 (MQ=255)
                                   caaggtgactgggtgatcctgaccaagggtgacagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCACAGCACAAGGCT‑‑‑GCGCTTGAGCCTtgt                                                                                                                  <  1:2099233‑M1/31‑1 (MQ=255)
                                        gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCtg                                                                                                              >  1:1838615/1‑141 (MQ=255)
                                        gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCtg                                                                                                              <  2:1838615/141‑1 (MQ=255)
                                         catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                                                    <  2:1218654/150‑1 (MQ=255)
                                         catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCCACATCAGGGGCTTTTTTATTGCTGGCGTGGCCTCTTCTCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                                                    <  1:870884/150‑1 (MQ=255)
                                           tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                                                  <  2:1151001/150‑1 (MQ=255)
                                              ggtgatcctgaccaagggtgacagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCGCAGCACaa                                                                                                                                         <  1:1124183‑M1/12‑1 (MQ=255)
                                              ggtgatcctgaccaagggtgacagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCGCAGCACaa                                                                                                                                         >  2:1124183‑M1/101‑112 (MQ=255)
                                                 gcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAc                                                                                            >  1:1386427/1‑150 (MQ=255)
                                                   atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCg                                                                                          >  1:686059/1‑150 (MQ=255)
                                                                 gTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAcc                                                                            >  2:2272817/1‑150 (MQ=255)
                                                                  acagctaccacaccatcggtggcaccaatggcatgaagatcctgcacgtcggtgatccgctggtcggttgatccactcgtTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAg                                                                                                     <  2:2229423‑M1/44‑1 (MQ=255)
                                                                     gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACg                                                                        >  2:1891032/1‑150 (MQ=255)
                                                                     gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACg                                                                        >  1:1711976/1‑150 (MQ=255)
                                                                       cagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTg                                                                      >  2:1509195/1‑150 (MQ=255)
                                                                        agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCtg                                                                                                              >  1:107949/1‑109 (MQ=255)
                                                                        agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCtg                                                                                                              <  2:107949/109‑1 (MQ=255)
                                                                                cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                                                    <  1:1396222/111‑1 (MQ=255)
                                                                                cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCtt                                                             >  1:2188090/1‑150 (MQ=255)
                                                                                   ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCa                                                                    >  1:640057/1‑140 (MQ=255)
                                                                                   ccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCa                                                                    <  2:640057/140‑1 (MQ=255)
                                                                                    cTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATg                                                         <  2:1386427/150‑1 (MQ=255)
                                                                                      gCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGCTAAACGTACCAACGGATCACCGACGTGCAGGATCTTCATGcc                                                       >  2:1133109/1‑150 (MQ=255)
                                                                                          cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGaaa                                                                                                         <  2:2270892/96‑1 (MQ=255)
                                                                                          cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTg                                                   <  1:1891032/150‑1 (MQ=255)
                                                                                          cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCCTGAGCCTTGTGCTGTGaaa                                                                                                         >  1:2270892/1‑96 (MQ=255)
                                                                                           aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTGTATTGCTGTCGGGACCTCGTAGCAGCACAAGGCT‑‑‑GGGCTTGAGCCTTGTGCTGTGGAAAGATTGGGTCAACCGACGAGCGAATAACAGACATGCAGGAGCTTCATGCCATTgg                                                  >  2:648878/1‑150 (MQ=255)
                                                                                            tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGGTCTTCATGCGATTgtt                                                 >  1:965971/1‑148 (MQ=255)
                                                                                            tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                                                 <  1:1860969/150‑1 (MQ=255)
                                                                                            tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                                                 >  1:737402/1‑150 (MQ=255)
                                                                                            tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                                                 >  1:2307248/1‑150 (MQ=255)
                                                                                              aaaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                                                    <  1:831854/97‑1 (MQ=255)
                                                                                              aaaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                                                    >  2:831854/1‑97 (MQ=255)
                                                                                              aaaaaGCCCCGGTTCGCGAGAACCCGCTACGTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc                                               <  1:1170693/150‑1 (MQ=255)
                                                                                               aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGc                                                        >  2:1569629/1‑140 (MQ=255)
                                                                                               aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGc                                                        <  1:1569629/140‑1 (MQ=255)
                                                                                                aaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                                                  <  1:339717/97‑1 (MQ=255)
                                                                                                aaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                                                  >  2:339717/1‑97 (MQ=255)
                                                                                                 aaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc                                            <  2:1711976/150‑1 (MQ=255)
                                                                                                   gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTTAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCACGCCATTGGTGCCACCg                                          <  2:1169167/150‑1 (MQ=255)
                                                                                                   gCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCg                                          >  2:1760535/1‑150 (MQ=255)
                                                                                                    ccccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCtgt                                                                                                             >  1:62525/1‑82 (MQ=255)
                                                                                                    ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCTGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa                                         >  1:1650041/1‑150 (MQ=255)
                                                                                                    ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                                                          >  1:2291596/1‑133 (MQ=255)
                                                                                                     cccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGc                                                                                                                       <  2:1953924/71‑1 (MQ=38)
                                                                                                     cccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGc                                                                                                                       >  1:1953924/1‑71 (MQ=255)
                                                                                                      ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCCGCCACCGTAC‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg                                       >  2:1374825/1‑150 (MQ=255)
                                                                                                      ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg                                       >  1:2094644/1‑150 (MQ=255)
                                                                                                      ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg                                       >  2:1328730/1‑150 (MQ=255)
                                                                                                       cGGTTCGCGAGAACCGGGGCTTTTTTCTGTGAGGCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATgg                                      >  1:1814101/1‑150 (MQ=255)
                                                                                                       cGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATgg                                      <  2:1650041/150‑1 (MQ=255)
                                                                                                        ggTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt                                     <  1:1509195/150‑1 (MQ=255)
                                                                                                         gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGaaa                                                                                                         <  1:518205/81‑1 (MQ=255)
                                                                                                         gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGaaa                                                                                                         >  2:518205/1‑81 (MQ=255)
                                                                                                         gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtg                                    <  2:522621/150‑1 (MQ=255)
                                                                                                                  aaCCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTg                           <  2:1396892/150‑1 (MQ=255)
                                                                                                                     cGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGCGTGGATCAACCGACCCGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGCGGTAGCTGTcc                        >  1:1540809/1‑149 (MQ=255)
                                                                                                                        ggCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCAccc                     >  2:833347/1‑150 (MQ=255)
                                                                                                                        ggCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCAccc                     <  2:686059/150‑1 (MQ=255)
                                                                                                                         gCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCt                    >  2:2039235/1‑150 (MQ=255)
                                                                                                                           ttttttATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCCTg                  <  1:1760535/150‑1 (MQ=255)
                                                                                                                               ttATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCa              <  2:2188090/150‑1 (MQ=255)
                                                                                                                                tATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAg             >  2:1500946/1‑150 (MQ=255)
                                                                                                                                     cTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCa        <  1:1500946/150‑1 (MQ=255)
                                                                                                                                     cTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCa        <  2:2307248/150‑1 (MQ=255)
                                                                                                                                        atccactcgtTTCACAGCACAAGGCTCAAGCGC‑‑‑AGCCTTGTGCTGCGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCAccc     <  1:1462484‑M1/140‑1 (MQ=255)
                                                                                                                                         tccactcgtTTCACAGCACAAGGCT‑‑‑CAAGCGCAGCCTTGTGCTGCGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc                                            <  1:1275747‑M1/101‑1 (MQ=255)
                                                                                                                                         cGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCACCCa    <  2:2094644/150‑1 (MQ=255)
                                                                                                                                         cGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCCTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGCGTGGTAGCTGTCACCCTTGGTCAGGATCACCCa    <  1:1485562/150‑1 (MQ=255)
                                                                                                                                          ggggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCt                    >  2:2298787/1‑133 (MQ=255)
                                                                                                                                          ggggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCt                    <  1:2298787/133‑1 (MQ=255)
                                                                                                                                          ggggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGATCAGGATCACCCAg   >  2:1887741/1‑150 (MQ=255)
                                                                                                                                           gggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt                                     <  2:577476/115‑1 (MQ=255)
                                                                                                                                           gggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGtgtg                                  >  1:779489/1‑118 (MQ=255)
                                                                                                                                           gggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGtgtg                                  <  2:779489/118‑1 (MQ=255)
                                                                                                                                           gggCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCACCCAGt  >  2:209008/1‑150 (MQ=255)
                                                                                                                                                tctTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                                                         >  2:1127241/1‑74 (MQ=255)
                                                                                                                                                tctTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                                                         <  1:1127241/74‑1 (MQ=255)
                                                                                                                                                  TTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAgg            <  2:574002‑M1/133‑1 (MQ=255)
                                                                                                                                                  TTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAgg            >  1:574002‑M1/1‑133 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCT‑‑‑GCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGGTAGCTGTCACCCTTGGTCAGGATCACCCAGT  >  NC_002947/1551215‑1551503

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