Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,969,020 A→C 5.9% I368I (ATA→ATC)  PP_3498 → transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,969,0200AC 5.9% 176.1 / 3.9 68I368I (ATA→ATC) PP_3498transposase
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base C (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAC  >  NC_002947/3968874‑3969154
                                                                                                                                                  |                                                                                                                                      
cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATgg                                                                                                                                     >  1:159249/1‑150 (MQ=32)
cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATgg                                                                                                                                     <  1:242893/150‑1 (MQ=32)
   atgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTTCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGcc                                                                                                                                  >  2:800820/1‑150 (MQ=17)
         cAGGGCGCGCATTCCTTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                                                                                                            >  2:1179147/1‑150 (MQ=17)
             gcgcgcATTCCTTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                                                                                                        <  1:1179147/150‑1 (MQ=17)
             gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGGGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                                                                                                        >  2:923585/1‑150 (MQ=17)
               gcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCt                                                                                                                                <  2:234350/140‑1 (MQ=33)
               gcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCt                                                                                                                                >  1:234350/1‑140 (MQ=33)
               gcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGc                                                                                                                      <  2:326586/150‑1 (MQ=32)
                 gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTatat                                                                                                                                       <  1:1150629/131‑1 (MQ=31)
                 gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTatat                                                                                                                                       >  2:1150629/1‑131 (MQ=31)
                    ttCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAGGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCa                                                                                                                 >  1:596667/1‑150 (MQ=17)
                     tCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAg                                                                                                                <  1:193114/150‑1 (MQ=32)
                       cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGt                                                                                                              <  2:159249/150‑1 (MQ=32)
                        gTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGtt                                                                                                             <  1:148690/150‑1 (MQ=32)
                         tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGttt                                                                                                            >  2:1186403/1‑150 (MQ=32)
                                  aGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCa                                                                                                   >  1:507481/1‑150 (MQ=32)
                                         gCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAgg                                                                                                               >  2:847216/1‑131 (MQ=33)
                                         gCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAgg                                                                                                               <  1:847216/131‑1 (MQ=33)
                                                tCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACAAAAAGaaaaa                                                                                     <  2:871467/150‑1 (MQ=12)
                                                   gATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                                  >  2:1205734/1‑150 (MQ=32)
                                                   gATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                                  >  2:1207954/1‑150 (MQ=32)
                                                   attagccagcctgatcgccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatca‑TCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                                                                                                       >  2:313868‑M2/68‑128 (MQ=255)
                                                   attagccagcctgatcgccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatca‑TCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                                                                                                       <  1:313868‑M2/61‑1 (MQ=255)
                                                    aTCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCgg                                                                                 >  2:1147697/1‑150 (MQ=32)
                                                     tCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                <  1:1048733/150‑1 (MQ=32)
                                                           aGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcggc                                                                          >  1:1020363/1‑150 (MQ=32)
                                                               gatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATg                                                                      >  1:431677/1‑150 (MQ=32)
                                                                 tcgccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAg                                                                    >  1:718776‑M2/54‑150 (MQ=255)
                                                                 tCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAg                                                                    >  1:701509/1‑150 (MQ=32)
                                                                    ccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGa                                                                     >  1:1006771‑M2/51‑146 (MQ=255)
                                                                    ccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGa                                                                     <  2:1006771‑M2/96‑1 (MQ=255)
                                                                    ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATg                                                                                                                                      >  1:501005/1‑81 (MQ=33)
                                                                    ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATg                                                                                                                                      <  2:501005/81‑1 (MQ=33)
                                                                     cgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCctc                                                                <  1:363708/150‑1 (MQ=32)
                                                                        ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                             <  1:1205734/150‑1 (MQ=32)
                                                                         tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                                           <  2:447118/135‑1 (MQ=32)
                                                                         tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                                           >  1:447118/1‑135 (MQ=33)
                                                                         tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGAc                                                            <  1:1174353/150‑1 (MQ=32)
                                                                         tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGAc                                                            >  1:160815/1‑150 (MQ=32)
                                                                         tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATATGCGGATTATATGTGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGAc                                                            >  1:979703/1‑150 (MQ=17)
                                                                         tccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                                                                                            >  2:839241‑M2/46‑118 (MQ=255)
                                                                         tccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                                                                                            <  1:839241‑M2/73‑1 (MQ=255)
                                                                                cTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCgg                                                     <  2:160815/150‑1 (MQ=32)
                                                                                  ccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                             <  2:399435‑M2/104‑1 (MQ=255)
                                                                                  ccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                             >  1:399435‑M2/37‑140 (MQ=255)
                                                                                         ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCgg                                                     <  1:1278813/141‑1 (MQ=32)
                                                                                         ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCgg                                                     >  2:1278813/1‑141 (MQ=32)
                                                                                         ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAAACCGAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCa                                            <  1:977711/150‑1 (MQ=12)
                                                                                                     cAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAg                                <  1:27380/150‑1 (MQ=32)
                                                                                                      caagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAgg                               <  2:718776‑M2/134‑1 (MQ=255)
                                                                                                         gaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGgcgc                            <  2:701509‑M2/137‑1 (MQ=255)
                                                                                                         gaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGgcgc                            <  1:923585‑M2/137‑1 (MQ=255)
                                                                                                              ggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCaa                       <  2:507481‑M2/142‑1 (MQ=255)
                                                                                                              ggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTGTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCaa                       <  2:1020363‑M2/142‑1 (MQ=255)
                                                                                                               ccgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTAACAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAAc                      <  1:848436‑M2/143‑1 (MQ=255)
                                                                                                                cgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACa                                                                                                                          <  2:166052‑M2/43‑1 (MQ=255)
                                                                                                                cgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACa                                                                                                                          >  1:166052‑M2/7‑49 (MQ=255)
                                                                                                                  atcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAAcccc                   <  1:800820‑M2/146‑1 (MQ=255)
                                                                                                                      ggctggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                <  2:215011/85‑1 (MQ=35)
                                                                                                                      ggctggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                >  1:215011/1‑85 (MQ=35)
                                                                                                                       gctggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTa              >  1:170179/1‑150 (MQ=32)
                                                                                                                           gctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAgg                               >  2:128271/1‑129 (MQ=33)
                                                                                                                           gctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAgg                               <  1:128271/129‑1 (MQ=33)
                                                                                                                           gctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGa          >  1:221316/1‑150 (MQ=32)
                                                                                                                             tAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAAc        >  1:662971/1‑150 (MQ=32)
                                                                                                                                tCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGTGGCTATGAGCCTCTGACCACCCGGGAAAGGCCACCGGGGCGTTAGGCGCTCCAACCCCTGCTAGTGAACAtt     >  1:1049155/1‑150 (MQ=11)
                                                                                                                                  gAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTATACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGa   <  2:662971/150‑1 (MQ=17)
                                                                                                                                         gCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAAcc                                                       >  2:1060260/1‑91 (MQ=33)
                                                                                                                                         gCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAAcc                                                       <  1:1060260/91‑1 (MQ=34)
                                                                                                                                                 tatGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCCGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc  <  2:312444/136‑1 (MQ=14)
                                                                                                                                                 tatGGGCCCTGTTACAGAGCGATCAGGCGTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc  >  1:312444/1‑136 (MQ=14)
                                                                                                                                                  atGGGACCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCgggg                                      >  1:945654/1‑99 (MQ=14)
                                                                                                                                                  atGGGACCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCgggg                                      <  2:945654/99‑1 (MQ=14)
                                                                                                                                                  |                                                                                                                                      
CTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAC  >  NC_002947/3968874‑3969154

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