Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,816,637 T→A 6.7% I21N (ATT→AAT)  PP_2472 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,816,6370TA 6.7% 214.2 / 3.0 72I21N (ATT→AAT) PP_2472MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (29/37);  new base A (3/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTGC  >  NC_002947/2816489‑2816784
                                                                                                                                                    |                                                                                                                                                   
ccAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAtt                                                                                                                                                    <  1:969982/150‑1 (MQ=255)
 cAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTg                                                                                                                                                   <  1:1208442/150‑1 (MQ=255)
       aGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATGGTGAAg                                                                                                                                             <  2:1379027/150‑1 (MQ=255)
        ggggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCCACTTCACCATTGGTGAAGt                                                                                                                                            <  2:504154/150‑1 (MQ=255)
           gAAGAGACCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACCAGGAACTTCCGCCCATACCGGGCAAACGCCACTTCACCAATGGTGAAGtagc                                                                                                                                         <  1:1279858/150‑4 (MQ=255)
              gagaTCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGa                                                                                                                                      <  2:1449631/150‑1 (MQ=255)
                  tCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTAGGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtg                                                                                                                                  >  2:1058830/1‑150 (MQ=255)
                   cccGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgt                                                                                                                                 <  1:1743155/150‑1 (MQ=255)
                     cGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGtt                                                                                                                                   <  1:928860/146‑1 (MQ=255)
                     cGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGtt                                                                                                                                   >  2:928860/1‑146 (MQ=255)
                         cACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCg                                                                                                                           <  1:337250/150‑1 (MQ=255)
                           ccGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTa                                                                                                                         >  2:1713698/1‑150 (MQ=255)
                             gCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTaaa                                                                                                                       <  1:955035/150‑1 (MQ=255)
                             gCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTaaa                                                                                                                       <  1:170905/150‑1 (MQ=255)
                               aCGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGc                                                                                                                     <  2:389717/150‑1 (MQ=255)
                                   cgcgTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCcaca                                                                                                                 <  1:1631623/150‑1 (MQ=255)
                                    gcgTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCacac                                                                                                                >  2:1464169/1‑150 (MQ=255)
                                       tcgtcACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGtg                                                                                                             >  1:1495980/1‑150 (MQ=255)
                                        cgtcACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgtgtg                                                                                                                              >  2:910597/1‑132 (MQ=255)
                                        cgtcACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgtgtg                                                                                                                              <  1:910597/132‑1 (MQ=255)
                                              cTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGt                                                                                                      >  2:954139/1‑150 (MQ=255)
                                                   gtccgaggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATGGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                 <  2:461327/144‑1 (MQ=255)
                                                   gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                 <  1:954139/150‑1 (MQ=255)
                                                    tCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCACTCGTGAAGTGAGCGAGTTATGTGCCGTAAAGCCACACGTGCTGCGGTATTggg                                                                                                >  2:660464/1‑150 (MQ=255)
                                                        gggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCa                                                                                            <  1:1713698/150‑1 (MQ=255)
                                                         ggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTCAAGCCACACGTGCTGCGGTATTGGGAGCAg                                                                                           >  2:683661/1‑150 (MQ=255)
                                                         ggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAg                                                                                           >  1:1611863/1‑150 (MQ=255)
                                                         ggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAg                                                                                           <  2:273801/150‑1 (MQ=255)
                                                           cAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGa                                                                                         <  2:754435/150‑1 (MQ=255)
                                                           cAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGa                                                                                         <  2:1341066/150‑1 (MQ=255)
                                                               aGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTATAGCCACACGTGCTGCCCTTCTTTGAGCAGGAGTTc                                                                                     <  1:543506/150‑1 (MQ=255)
                                                                   gAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGAcc                                                                                 >  2:337161/1‑150 (MQ=255)
                                                                       gCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCt                                                                             >  2:1728735/1‑150 (MQ=255)
                                                                         ccGGGTTGAGGCCTATGGTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                 <  2:1117781/128‑1 (MQ=255)
                                                                         ccGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                 >  1:1117781/1‑128 (MQ=255)
                                                                             gTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAAccc                                                                       <  1:1232509/150‑1 (MQ=255)
                                                                              ttGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCt                                                                      <  1:547156/150‑1 (MQ=255)
                                                                                  ggCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGa                                                                  >  1:1518555/1‑150 (MQ=255)
                                                                                      tATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAgcg                                                              >  1:573099/1‑150 (MQ=255)
                                                                                          cTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGcggcgg                                                          <  2:1518555/150‑1 (MQ=255)
                                                                                          cTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGcggcgg                                                          <  1:1414814/150‑1 (MQ=255)
                                                                                               aCCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGaa                                                     >  1:684308/1‑150 (MQ=255)
                                                                                                   aGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcgg                                                 >  1:1370096/1‑150 (MQ=255)
                                                                                                      cATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGtt                                                                                      <  1:1298183/110‑1 (MQ=255)
                                                                                                      cATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGtt                                                                                      >  2:1298183/1‑110 (MQ=255)
                                                                                                      cATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggcgg                                              <  1:1421067/150‑1 (MQ=255)
                                                                                                           cgacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTa                                         <  1:1728735/150‑1 (MQ=255)
                                                                                                           cgacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTa                                         <  2:1520528/150‑1 (MQ=255)
                                                                                                            gacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGTAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTAc                                        >  1:961624/1‑150 (MQ=255)
                                                                                                            gacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGaa                                                     <  2:542980/137‑1 (MQ=255)
                                                                                                            gacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGaa                                                     >  1:542980/1‑137 (MQ=255)
                                                                                                            gacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTAc                                        <  1:298027/150‑1 (MQ=255)
                                                                                                                     ccccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAg                               <  2:792956/150‑1 (MQ=255)
                                                                                                                     ccccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAg                               >  1:1153970/1‑150 (MQ=255)
                                                                                                                        cccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAt                            <  2:1294258/150‑1 (MQ=255)
                                                                                                                        cccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAt                            >  2:1064274/1‑150 (MQ=255)
                                                                                                                        cccATACCGGGCAAACGCTACTTCACCATTGGGGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAt                            >  2:1071673/1‑150 (MQ=255)
                                                                                                                         ccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGtt                                                                                                                                   >  2:928934/1‑46 (MQ=255)
                                                                                                                         ccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGtt                                                                                                                                   <  1:928934/46‑1 (MQ=255)
                                                                                                                               cGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtga                     <  1:667832/150‑1 (MQ=255)
                                                                                                                                 ggCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtgatg                   >  2:1426365/1‑150 (MQ=255)
                                                                                                                                       cGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGc             <  1:1426365/150‑1 (MQ=255)
                                                                                                                                         tcaCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCa           >  1:829320/3‑150 (MQ=255)
                                                                                                                                          caCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAg          >  1:1489167/2‑150 (MQ=255)
                                                                                                                                           aCTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGa         <  1:1033976/150‑1 (MQ=255)
                                                                                                                                           aCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGa         >  2:70633/1‑150 (MQ=255)
                                                                                                                                             ttCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGAtt       >  1:1306814/1‑150 (MQ=255)
                                                                                                                                             ttCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGAtt       >  1:1074910/1‑150 (MQ=255)
                                                                                                                                               cACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCg     >  1:1858872/1‑150 (MQ=255)
                                                                                                                                               cACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCg     >  1:1858705/1‑150 (MQ=255)
                                                                                                                                                aCCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGt    >  1:1127441/1‑150 (MQ=255)
                                                                                                                                                 ccATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTg   <  2:1611863/150‑1 (MQ=255)
                                                                                                                                                  cATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTgc  <  2:295613/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTGC  >  NC_002947/2816489‑2816784

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