Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,480,664 G→C 7.4% S166R (AGC→AGG) ‡ kefG ← protein required for KefB activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,480,6640GC 7.4% 310.6 / 13.1 94S166R (AGC→AGG) ‡kefGprotein required for KefB activity
Reads supporting (aligned to +/- strand):  ref base G (37/50);  new base C (3/4);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GCGCCAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACT  >  NC_000913/3480516‑3480813
                                                                                                                                                    |                                                                                                                                                     
gcgcCAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGc                                                                                                                                                      <  2:1315451/150‑1 (MQ=255)
 cgcCAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCt                                                                                                                                                     <  1:1419892/150‑1 (MQ=255)
 cgcCAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCt                                                                                                                                                     <  2:395286/150‑1 (MQ=255)
   ccAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCg                                                                                                                                                   >  2:1437933/1‑150 (MQ=255)
   ccAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCg                                                                                                                                                   >  2:449046/1‑150 (MQ=255)
     aaTACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGcc                                                                                                                                                 <  1:366995/150‑1 (MQ=255)
     aaTACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGcc                                                                                                                                                 <  1:1115305/150‑1 (MQ=255)
       tACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAg                                                                                                                                               <  1:147074/150‑1 (MQ=255)
       tACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGAGCCACGCCAg                                                                                                                                               <  1:641949/150‑1 (MQ=255)
       tACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGAGCCACGCCAg                                                                                                                                               <  2:302026/150‑1 (MQ=255)
         cccAGCCGCGATTCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGATTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGAGCCACGCCAgag                                                                                                                                             <  2:755381/150‑3 (MQ=255)
         cccAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCCCCGTAGGCTCTGGCGTGGCTCGCCAGCt                                                                                                                                             >  1:484851/1‑150 (MQ=255)
         cccAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCt                                                                                                                                             <  2:959737/150‑1 (MQ=255)
          ccAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTc                                                                                                                                            >  2:911783/1‑150 (MQ=255)
           cAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGACACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGAGCCACGCCAgagcc                                                                                                                                           <  1:899611/150‑5 (MQ=255)
           cAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTcc                                                                                                                                           <  1:1365562/150‑1 (MQ=255)
               cgcgATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCtgtg                                                                                                                                       <  1:699508/150‑1 (MQ=255)
                     gCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCttt                                                                                                                                 >  2:1646629/1‑150 (MQ=255)
                      ccAGCGGCACCGCAGCCACCGCCGCGATGAGAATCAGCACTCGTGCGAGTAAAAAATGGGAAGCTTCCATCAGCGGCCTCCTGGAGAGAGCGGATTTGCCAGCCAGTCAGCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTg                                                                                                                                <  1:494938/150‑1 (MQ=255)
                       cAGCGGCACCGCAGCCCCCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGc                                                                                                                               >  2:271322/1‑150 (MQ=255)
                         gCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCg                                                                                                                             >  1:146436/1‑150 (MQ=255)
                            gCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGtct                                                                                                                          >  1:440245/1‑150 (MQ=255)
                             cACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGtctc                                                                                                                         >  2:102868/1‑150 (MQ=255)
                                  cAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCa                                                                                                                    <  2:919375/150‑1 (MQ=255)
                                      cACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCa                                                                                                                                                >  1:501878/1‑118 (MQ=255)
                                      cACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCa                                                                                                                                                <  2:501878/118‑1 (MQ=255)
                                      cACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTaa                                                                                                                <  1:911783/150‑1 (MQ=255)
                                      cACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTaa                                                                                                                <  2:100948/150‑1 (MQ=255)
                                        ccgccgCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAaat                                                                                                              <  1:95040/150‑1 (MQ=255)
                                            cgcgAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAatg                                                                                                          <  2:818521/150‑1 (MQ=255)
                                             gcgAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAatga                                                                                                         <  1:1023489/150‑1 (MQ=255)
                                               gAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAatgatg                                                                                                       <  1:1094875/150‑1 (MQ=255)
                                                 agagAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATggg                                                                                                     <  1:1003227/150‑1 (MQ=255)
                                                    gaAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACt                                                                                                  >  1:1615598/1‑150 (MQ=255)
                                                      aaCAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTa                                                                                                >  2:1440320/1‑150 (MQ=255)
                                                      aaCAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTa                                                                                                >  2:233813/1‑150 (MQ=255)
                                                      aaCAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTa                                                                                                >  1:1460734/1‑150 (MQ=255)
                                                        cAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAAc                                                                                              >  2:1561538/1‑150 (MQ=255)
                                                          gCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCa                                                                                            >  1:442839/1‑150 (MQ=255)
                                                            aCTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAAt                                                                                          >  1:1415077/1‑150 (MQ=255)
                                                               ccTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCa                                                                                       <  2:1472411/150‑1 (MQ=255)
                                                                  gCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATcc                                                                                    >  1:654591/1‑150 (MQ=255)
                                                                     aGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACATCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTTGCTCGCCAGCTCCTGTGCGCTTTGCCGGCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGc                                                                                 <  1:586170/150‑1 (MQ=255)
                                                                     aGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGc                                                                                 >  1:1562229/1‑150 (MQ=255)
                                                                          aaaaTCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATg                                                                            >  1:29196/1‑150 (MQ=255)
                                                                           aaaTCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGc                                                                           >  2:11576/1‑150 (MQ=255)
                                                                           aaaTCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGc                                                                           <  1:1143233/150‑1 (MQ=255)
                                                                                   aCCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCgccgcc                                                                   <  2:1460734/150‑1 (MQ=255)
                                                                                   aCCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCgccgcc                                                                   <  2:192830/150‑1 (MQ=255)
                                                                                     cTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTTTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAg                                                                 >  1:410539/1‑150 (MQ=255)
                                                                                     cTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAAAAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAg                                                                 >  2:838030/1‑150 (MQ=255)
                                                                                      ttCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGt                                                                <  1:1143924/150‑1 (MQ=255)
                                                                                        ccATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTc                                                              <  1:205017/150‑1 (MQ=255)
                                                                                        ccATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTc                                                              >  1:564454/1‑150 (MQ=255)
                                                                                         cATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCa                                                             >  2:1426516/1‑150 (MQ=255)
                                                                                                 gCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGACAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGgcgc                                                    <  2:36217/150‑1 (MQ=255)
                                                                                                   ctcctGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCgca                                                   <  1:1640337/150‑1 (MQ=255)
                                                                                                     cctGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTTTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCgcagc                                                 <  2:410539/150‑1 (MQ=255)
                                                                                                     cctGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCgcagc                                                 <  2:927833/150‑1 (MQ=255)
                                                                                                     cctGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCgcagc                                                 <  2:1425113/150‑1 (MQ=255)
                                                                                                      ctGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCgcagca                                                <  1:1461728/150‑1 (MQ=255)
                                                                                                      ctGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCgcagca                                                <  1:449046/150‑1 (MQ=255)
                                                                                                          agaCAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCCATGCATCCGGCACATGCCCGCCGCCAGGTCAAAGGGGCGCAGCACATc                                            >  2:789620/1‑150 (MQ=255)
                                                                                                          agaCAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATc                                            <  2:442839/150‑1 (MQ=255)
                                                                                                          agaCAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATc                                            >  2:789618/1‑150 (MQ=255)
                                                                                                           gaCAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCg                                           <  1:1646629/150‑1 (MQ=255)
                                                                                                            aCAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGc                                          <  1:271322/150‑1 (MQ=255)
                                                                                                             cAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCt                                         >  2:620482/1‑150 (MQ=255)
                                                                                                               gCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCa                                       >  2:725879/1‑150 (MQ=255)
                                                                                                                 ggATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATc                                     <  1:409892/150‑1 (MQ=255)
                                                                                                                 ggATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGAAGGGACTTAACCAATGCATCCGGCACATGCCCGACGCCAGTTCAAAGGGGAGCAGCACATCGCTccta                                     >  1:710589/1‑147 (MQ=255)
                                                                                                                   aTTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCgg                                   >  2:1454449/1‑150 (MQ=255)
                                                                                                                   aTTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCgg                                   >  2:338084/1‑150 (MQ=255)
                                                                                                                    tttGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCggg                                  <  2:146436/150‑1 (MQ=255)
                                                                                                                    tttGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTAGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCggg                                  <  2:380660/150‑1 (MQ=255)
                                                                                                                      tGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTa                                <  2:29196/150‑1 (MQ=255)
                                                                                                                        ccagccagTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGc                              <  2:1274297/150‑1 (MQ=255)
                                                                                                                        ccagccagTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGc                              >  2:810634/1‑150 (MQ=255)
                                                                                                                         cagccagTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCg                             <  1:233813/150‑1 (MQ=255)
                                                                                                                         cagccagTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCg                             >  1:1224145/1‑150 (MQ=255)
                                                                                                                           gccagTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGAt                           <  1:11576/150‑1 (MQ=255)
                                                                                                                              agTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCa                        >  1:741590/1‑150 (MQ=255)
                                                                                                                               gTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAg                       <  1:1316221/150‑1 (MQ=255)
                                                                                                                                     gTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTc                 <  1:620482/150‑1 (MQ=255)
                                                                                                                                     acaggagctGGCGAGCCACGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAGCCAATGCCTCCGGCACATGCCCGCCGCCAGGTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTc                 >  1:1658121/8‑150 (MQ=255)
                                                                                                                                      gcggCTCTGGCGTAGCTCGCCAGCTCCTGTGCGCTTTGCCGGCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCa                <  1:1170445/148‑1 (MQ=255)
                                                                                                                                        ggagctGGCGAGCCACGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATa              >  1:425632/5‑150 (MQ=255)
                                                                                                                                        ggagctGGCGAGCCACGCCAGAGCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCCCATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATa              >  1:1298493/5‑150 (MQ=255)
                                                                                                                                        ggCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATa              >  1:906770/1‑150 (MQ=255)
                                                                                                                                            ctGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACgg          <  1:1440320/150‑1 (MQ=255)
                                                                                                                                              ggCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTa        >  2:1329926/1‑150 (MQ=255)
                                                                                                                                                cGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAg      <  1:789618/150‑1 (MQ=255)
                                                                                                                                                cGTGGCGCGCCAGCTCCTGTCCGCTGTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGCGGCGCAGCACAGCGCTCATCGGGTAGCGATTCAGCGCGTCATAACCGTAAg      <  1:789620/150‑1 (MQ=255)
                                                                                                                                                  tGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCGATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCCTAACGGTAAGCa    <  1:260084/150‑1 (MQ=255)
                                                                                                                                                    gCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACt  >  1:1393479/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCCAATACCCAGCCGCGATGCCAGCGGCACCGCAGCCACCGCCGCGAAGAGAAACAGCACTCCTGCGAGTAAAAAATCGGAACCTTCCATCAGCGGCCTCCTGGAGACAGCGGATTTGCCAGCCAGTCACCGTAGGCTCTGGCGTGGCTCGCCAGCTCCTGTGCGCTTTGCCGTCTCGCCCAGTAAATAATGATGGGACTTAACCAATGCATCCGGCACATGCCCGCCGCCAGTTCAAAGGGGCGCAGCACATCGCTCATCGGGTAGCGATTCAGCGCGTCATAACGGTAAGCACT  >  NC_000913/3480516‑3480813

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