Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,907 T→C 6.1% I227I (ATT→ATC yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,9070TC 6.1% 459.8 / 12.1 146I227I (ATT→ATCyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (59/78);  new base C (6/3);  total (65/81)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 160 total aligned reads displayed.

TCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCT  >  NC_002947/1486758‑1487051
                                                                                                                                                     |                                                                                                                                                
tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                                                  >  2:1193960/1‑150 (MQ=255)
tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                                                  >  1:2727218/1‑150 (MQ=255)
tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                                                  >  2:2483813/1‑150 (MQ=255)
 cGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtg                                                                                                                                                 <  2:1313655/150‑1 (MQ=255)
  ggCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgt                                                                                                                                                <  1:1077059/150‑1 (MQ=255)
  ggCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgt                                                                                                                                                >  2:2715506/1‑150 (MQ=255)
  ggCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgt                                                                                                                                                >  1:2012626/1‑150 (MQ=255)
   gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg                                                                                                                                               <  1:3403146/150‑1 (MQ=255)
   gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg                                                                                                                                               <  2:1935590/150‑1 (MQ=255)
   gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg                                                                                                                                               <  2:2226378/150‑1 (MQ=255)
   gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg                                                                                                                                               <  2:2851776/150‑1 (MQ=255)
     aCTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGa                                                                                                                                             <  1:3080033/150‑1 (MQ=255)
      ctcgACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGACGGCATCACAATCATCTGCGTAACGCATGAAATGGGTTTTGCGCGCAAGATCGCCGATCGCGTCATCTTTATGGATCAGGGCAAGATCGTCGAg                                                                                                                                            <  1:5387771/147‑1 (MQ=255)
       tGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGc                                                                                                                                           >  1:2496736/1‑150 (MQ=255)
       tGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGc                                                                                                                                           <  2:395364/150‑1 (MQ=255)
         gatccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAg                                                                                                                                         >  2:1480009/4‑150 (MQ=255)
         gatccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTCGAACAg                                                                                                                                         >  2:3786894/4‑150 (MQ=255)
         gatccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTCGAACAg                                                                                                                                         >  2:1046939/4‑150 (MQ=255)
         gACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGACGGCATCACAATCATCTGCGTAACGCATGAAATGGGTTTTGCGCGCAAGATCGCCGATCGCGTCATCTTTATGGATCAGGGCAAGATCGTCGAggaa                                                                                                                                         >  2:1759407/1‑147 (MQ=255)
          atccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAga                                                                                                                                        <  2:3319847/148‑2 (MQ=255)
          aCCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAgg                                                                                                                                        <  2:1456016/150‑1 (MQ=255)
           tccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAgaa                                                                                                                                       >  1:1125438/2‑148 (MQ=255)
            ccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGcc                                                                                                                                      >  2:2917055/1‑150 (MQ=255)
              ggAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAgaattc                                                                                                                                    >  2:5154691/1‑145 (MQ=255)
              ggAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCcgc                                                                                                                                    >  1:1745607/1‑150 (MQ=255)
              ggAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCcgc                                                                                                                                    <  2:1322586/150‑1 (MQ=255)
                 aaTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGcc                                                                                                                                 <  2:4720706/150‑1 (MQ=255)
                 aaTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGcc                                                                                                                                 >  1:161881/1‑150 (MQ=255)
                     gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             >  1:978644/1‑150 (MQ=255)
                     gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             >  2:4385062/1‑150 (MQ=255)
                     gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             >  1:3720333/1‑150 (MQ=255)
                      tCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATg                                                                                                                            >  1:3793011/1‑150 (MQ=255)
                        aaGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGAc                                                                                                                          <  1:2828245/150‑1 (MQ=255)
                          ggAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActt                                                                                                                        >  2:1556243/1‑150 (MQ=255)
                             aGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActtctt                                                                                                                     <  2:3574062/150‑1 (MQ=255)
                             aGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActtctt                                                                                                                     >  2:2732454/1‑150 (MQ=255)
                              gTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                                                    <  1:4184614/150‑1 (MQ=255)
                              gTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                                                    <  1:3715159/150‑1 (MQ=255)
                               tGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCAGGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCg                                                                                                                   <  2:2363109/150‑1 (MQ=255)
                                gCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGa                                                                                                                  >  1:2832973/1‑150 (MQ=255)
                                  tCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGAcc                                                                                                                <  2:1873028/150‑1 (MQ=255)
                                    gATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCgg                                                                                                              >  2:5399226/1‑150 (MQ=255)
                                             gtggGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGt                                                                                                     >  1:2984623/1‑150 (MQ=255)
                                                      gCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGAcc                                                                                            <  1:3870417/150‑1 (MQ=255)
                                                        cGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCac                                                                                          >  1:406512/1‑149 (MQ=255)
                                                        cGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCaa                                                                                          >  1:3579926/1‑150 (MQ=255)
                                                           ggATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAgtt                                                                                       <  1:308318/150‑1 (MQ=255)
                                                            gATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAgttg                                                                                      <  2:4445459/150‑1 (MQ=255)
                                                             aTGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGTCCGGCCGCGCAGTGACCGGACCAAgttgt                                                                                     <  2:2068065/150‑1 (MQ=255)
                                                             aTGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAgttgt                                                                                     >  2:3147446/1‑150 (MQ=255)
                                                             aTGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAgttgt                                                                                     <  1:741987/150‑1 (MQ=255)
                                                                 cATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACTAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCt                                                                                 <  2:1578987/150‑1 (MQ=255)
                                                                 cATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCt                                                                                 <  1:2917050/150‑1 (MQ=255)
                                                                   tGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCa                                                                               <  1:5133311/150‑1 (MQ=255)
                                                                     aCCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGc                                                                             <  1:4821745/150‑1 (MQ=255)
                                                                      ccATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGcc                                                                            <  1:1556243/150‑1 (MQ=255)
                                                                      ccATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGcc                                                                            >  2:5208982/1‑150 (MQ=255)
                                                                         tGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGa                                                                         <  1:5147808/150‑1 (MQ=255)
                                                                            tctGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATcc                                                                      >  2:200100/1‑150 (MQ=255)
                                                                             ctGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCt                                                                     <  1:5208982/150‑1 (MQ=255)
                                                                               gCGTGACCCACGAATTGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGc                                                                   <  2:3824374/150‑1 (MQ=255)
                                                                                cGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCg                                                                                                               >  1:4457171/1‑105 (MQ=255)
                                                                                cGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCg                                                                                                               <  2:4457171/105‑1 (MQ=255)
                                                                                 gTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAt                                                                 <  1:4603006/150‑1 (MQ=255)
                                                                                    aCCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGa                                                              <  1:1821061/150‑1 (MQ=255)
                                                                                      ccACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATg                                                            >  2:1117834/1‑150 (MQ=255)
                                                                                      ccACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATg                                                            >  2:2128933/1‑150 (MQ=255)
                                                                                        aCGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGcc                                                                                                            <  2:1707591/100‑1 (MQ=255)
                                                                                        aCGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGcc                                                                                                            >  1:1707591/1‑100 (MQ=255)
                                                                                            aaTGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGttt                                                      <  1:1583733/150‑1 (MQ=255)
                                                                                             aTGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTg                                                     <  1:2583821/150‑1 (MQ=255)
                                                                                                 gCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTg                                                     <  2:228383/146‑1 (MQ=255)
                                                                                                 gCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTg                                                     >  1:228383/1‑146 (MQ=255)
                                                                                                     tGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTgcgc                                             <  2:2041973/150‑1 (MQ=255)
                                                                                                          atttcgtggCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAG                                        >  2:2990542‑M1/5‑150 (MQ=255)
                                                                                                            aCCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCC                                      <  1:3877604‑M1/150‑1 (MQ=255)
                                                                                                            aCCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCC                                      <  1:4843714‑M1/150‑1 (MQ=255)
                                                                                                                tGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTC                                  <  2:362780‑M1/150‑1 (MQ=255)
                                                                                                                tGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTC                                  >  1:3273389‑M1/1‑150 (MQ=255)
                                                                                                                tGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCCGCGCAGCCCaaa                                  <  1:4919641‑M1/150‑4 (MQ=255)
                                                                                                                tGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGAGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTC                                  >  1:3273370‑M1/1‑150 (MQ=255)
                                                                                                                  gCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGC                                >  2:153857‑M1/1‑150 (MQ=255)
                                                                                                                     aaCCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGG                             >  2:4606603‑M1/1‑150 (MQ=255)
                                                                                                                        cGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTACTTCGACCGGCCGCGCAGTGACCGGACCAAGATGTTCCACAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCcc                          >  1:282192‑M1/1‑148 (MQ=255)
                                                                                                                            ttATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCC                      <  2:695598‑M1/150‑1 (MQ=255)
                                                                                                                              aTCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTggg                                                   >  2:363781/1‑119 (MQ=255)
                                                                                                                              aTCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTggg                                                   <  1:363781/119‑1 (MQ=255)
                                                                                                                                cTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTC                  <  1:1760502‑M1/150‑1 (MQ=255)
                                                                                                                                 ttCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCC                 <  2:4223023‑M1/150‑1 (MQ=255)
                                                                                                                                 ttCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCC                 <  1:610802‑M1/150‑1 (MQ=255)
                                                                                                                                  tCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCC                >  2:2073529‑M1/1‑150 (MQ=255)
                                                                                                                                   cATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCGGCCCTTCGCGGGCACGCCCGCTCCCA               <  2:3793011‑M1/150‑1 (MQ=255)
                                                                                                                                      ggACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG            <  1:815162‑M1/150‑1 (MQ=255)
                                                                                                                                       gACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAGTGATGCGGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGG           >  2:4936633‑M1/1‑150 (MQ=255)
                                                                                                                                          aaGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt        >  2:1232995/1‑150 (MQ=255)
                                                                                                                                          aaGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt        <  2:2400644/150‑1 (MQ=255)
                                                                                                                                          aaGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGAGGCGGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCGCCCCCAGGttt        >  2:535258/1‑150 (MQ=255)
                                                                                                                                           aGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTa       <  1:2073529/150‑1 (MQ=255)
                                                                                                                                               cGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAgcgc   >  1:5043362/1‑150 (MQ=255)
                                                                                                                                                gagaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCt  >  2:1313309/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCT  >  NC_002947/1486758‑1487051

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: