Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 338,071 A→C 7.3% S167R (AGT→AGG)  PP_0280 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947338,0710AC 7.3% 264.9 / 12.1 96S167R (AGT→AGG) PP_0280amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (40/49);  new base C (3/4);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

AGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTGCA  >  NC_002947/337922‑338219
                                                                                                                                                     |                                                                                                                                                    
agcagcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCa                                                                                                                                                      <  2:1758135/150‑1 (MQ=255)
 gcagcTTCAAGCCTTGCACCAGAATCAACGAGATCACCACGGAGAATACGCCTTCGGCGAAGAAGAGTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCAc                                                                                                                                                     >  2:25765/1‑150 (MQ=255)
   agcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTg                                                                                                                                                   <  1:981519/150‑1 (MQ=255)
      ttGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGcc                                                                                                                                                >  2:1747695/1‑150 (MQ=255)
       tGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGcca                                                                                                                                               <  2:717081/150‑2 (MQ=255)
       tGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGcca                                                                                                                                               <  1:223108/150‑2 (MQ=255)
        gAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                              >  2:2344746/1‑150 (MQ=255)
        gAAGCCTTGCACCAGCATGAAAGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                              <  1:2070394/150‑1 (MQ=255)
             cTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            <  1:1121414/147‑1 (MQ=255)
             cTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            >  2:1121414/1‑147 (MQ=255)
             cTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCa                                                                                                                                         <  1:1875467/150‑1 (MQ=255)
              ttGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCAAGGTGCATTGGGAGGTTCGATGCCTGTCAGTTGCAGCAGGGCAATGGGGCTGGCGAGGGCACTGGCCTTGAGCAt                                                                                                                                        >  1:43346/1‑150 (MQ=255)
                gCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                      <  2:962098/150‑1 (MQ=255)
                gCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                      <  2:448482/150‑1 (MQ=255)
                gCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                      <  2:2145351/150‑1 (MQ=255)
                gCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                      <  1:1747695/150‑1 (MQ=255)
                 cACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAg                                                                                                                                     <  1:1918874/150‑1 (MQ=255)
                  aCCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAgg                                                                                                                                    >  2:904788/1‑150 (MQ=255)
                    cAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGAt                                                                                                                                  >  1:676343/1‑150 (MQ=255)
                    cAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGAt                                                                                                                                  <  1:904788/150‑1 (MQ=255)
                     aGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGagc                                                                                                                                 >  2:384963/1‑148 (MQ=255)
                          gAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTc                                                                                                                            <  2:284622/150‑1 (MQ=255)
                           aaCGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCg                                                                                                                           >  1:1293778/1‑150 (MQ=255)
                            aCGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGt                                                                                                                          >  2:589830/1‑150 (MQ=255)
                            aCGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGt                                                                                                                          <  1:814214/150‑1 (MQ=255)
                              gagaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTg                                                                                                                        <  2:1029070/150‑1 (MQ=255)
                               agaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGc                                                                                                                       <  1:589830/150‑1 (MQ=255)
                                 aTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCt                                                                                                                      <  2:1091973/149‑1 (MQ=255)
                                 aTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCt                                                                                                                      >  1:1091973/1‑149 (MQ=255)
                                   cACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCAATGATGGGACGGACCATGCCGGTCAGTTCGAGCATGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTa                                                                                                                   >  1:1305931/1‑150 (MQ=255)
                                    aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            <  2:2295762/124‑1 (MQ=255)
                                    aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            >  1:2295762/1‑124 (MQ=255)
                                            gAACACGCCTGCGGCGAAGATGATTACACTAGGCTGGTAGCTGTGGTTTTTGATGTTACAGGCCATGCCGTTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCATTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAg                                                                                                          <  2:1104902/150‑1 (MQ=255)
                                             aaCACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAgg                                                                                                         <  2:282349/150‑1 (MQ=255)
                                             aaCACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGCGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGCGCATCAGGATCACTTCGTTGCTGTACGCCGGCAgg                                                                                                         <  2:1331090/150‑1 (MQ=255)
                                                  gCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGAt                                                                                                    <  1:1756937/150‑1 (MQ=255)
                                                    cTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGTCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgc                                                                                                  >  1:627324/1‑150 (MQ=255)
                                                    cTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgc                                                                                                  >  2:2031274/1‑150 (MQ=255)
                                                      caggcgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGTCAGGCCGATgcgc                                                                                                <  1:1855180/148‑1 (MQ=255)
                                                         gcgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGAGCGCa                                                                                             <  1:2031274/150‑1 (MQ=255)
                                                          cgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAg                                                                                            <  1:2023911/150‑1 (MQ=255)
                                                            aagaagATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcg                                                                                          >  1:922574/1‑150 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcgc                                                                                         <  2:939951/150‑1 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcgc                                                                                         <  1:2344746/150‑1 (MQ=255)
                                                              gaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGcgcg                                                                                        >  2:808381/1‑150 (MQ=255)
                                                              gaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGcgcg                                                                                        >  1:171658/1‑150 (MQ=255)
                                                                agaTTTCCACGGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGGCGGTCATTTTCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTg                                                                                      <  2:65745/150‑1 (MQ=255)
                                                                 gaTTTCCATCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTgg                                                                                     <  2:669605/150‑1 (MQ=255)
                                                                 gaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTgg                                                                                     <  1:999347/150‑1 (MQ=255)
                                                                   tttCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgca                                                                                   <  1:808381/150‑1 (MQ=255)
                                                                   tttCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgca                                                                                   <  1:1976405/150‑1 (MQ=255)
                                                                        aCCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCatg                                                                              >  1:1350310/1‑150 (MQ=255)
                                                                          cGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATACCTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCatgat                                                                            >  1:1330833/1‑150 (MQ=255)
                                                                           ggCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCatgatg                                                                           <  1:1136733/150‑1 (MQ=255)
                                                                              aggtaggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTCAACTTATGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAg                                                                        <  1:2090174/150‑1 (MQ=255)
                                                                                     tGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            <  2:297505/75‑1 (MQ=255)
                                                                                     tGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            >  1:297505/1‑75 (MQ=255)
                                                                                     tGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTg                                                                                                                        <  1:1140309/95‑1 (MQ=255)
                                                                                     tGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTg                                                                                                                        >  2:1140309/1‑95 (MQ=255)
                                                                                      gCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGt                                                                <  1:234001/150‑1 (MQ=255)
                                                                                        gggCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGttt                                                              >  2:2322898/1‑150 (MQ=255)
                                                                                        gggCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACGCGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGttt                                                              >  1:1879068/1‑150 (MQ=255)
                                                                                              tgatgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAc                                                        >  2:192028/1‑150 (MQ=255)
                                                                                              tgatgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAc                                                        <  1:384963/150‑1 (MQ=255)
                                                                                                   gTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGCAGAACAGCGTTTTGCCACGCGAc                                                   >  2:1396424/1‑150 (MQ=255)
                                                                                                       gggCCACGCCGGTCAGTTCCAGAAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGc                                               >  2:719564/1‑150 (MQ=255)
                                                                                                          ccATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCa                                            >  2:456411/1‑150 (MQ=255)
                                                                                                          ccATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCa                                            >  1:2161240/1‑150 (MQ=255)
                                                                                                                cGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGccc                                      >  1:208559/1‑150 (MQ=255)
                                                                                                                 ggTCAGTTCCAGCAGGGTCACGGTACTAGGCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCg                                     >  1:666215/1‑150 (MQ=255)
                                                                                                                  gTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGt                                    >  1:1707480/1‑150 (MQ=255)
                                                                                                                   tCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTg                                   >  2:2050474/1‑150 (MQ=255)
                                                                                                                   tCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGAGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTg                                   <  1:209072/150‑1 (MQ=255)
                                                                                                                      gTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCg                                >  2:2176069/1‑150 (MQ=255)
                                                                                                                       ttCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGc                               >  2:1906858/1‑150 (MQ=255)
                                                                                                                          cagcagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTAGTCGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTc                            <  1:25765/150‑1 (MQ=255)
                                                                                                                             cagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAAt                         <  1:2176069/150‑1 (MQ=255)
                                                                                                                             cagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAAt                         <  2:537865/150‑1 (MQ=255)
                                                                                                                                 gTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCa                     <  2:171658/150‑1 (MQ=255)
                                                                                                                                  tCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCAc                    <  2:1293778/150‑1 (MQ=255)
                                                                                                                                    aCGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCt                  >  1:606508/1‑150 (MQ=255)
                                                                                                                                        tcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAc                                                        <  2:2046885/104‑1 (MQ=255)
                                                                                                                                        tcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAc                                                        >  1:2046885/5‑108 (MQ=255)
                                                                                                                                        tcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAg              >  2:591271/5‑150 (MQ=255)
                                                                                                                                        tACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAg              <  2:922574/150‑1 (MQ=255)
                                                                                                                                         aCTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAgc             >  2:543866/1‑149 (MQ=255)
                                                                                                                                           aggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGa           >  2:1422653/2‑150 (MQ=255)
                                                                                                                                            ggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAt          <  1:591271/150‑1 (MQ=255)
                                                                                                                                            ggCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAt          <  1:1400053/150‑1 (MQ=255)
                                                                                                                                            ggCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAt          <  2:606508/150‑1 (MQ=255)
                                                                                                                                             gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGt                                                                <  2:1223013/95‑1 (MQ=255)
                                                                                                                                             gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGt                                                                >  1:1223013/1‑95 (MQ=255)
                                                                                                                                             gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATc         >  2:1757336/1‑150 (MQ=255)
                                                                                                                                              ccAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCg        <  1:192028/150‑1 (MQ=255)
                                                                                                                                               cAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGCTGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGc       <  2:1707480/150‑1 (MQ=255)
                                                                                                                                               cAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGGTGCAGCATGATGTAGAACCGCGTGTTGCCACTCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACACTTAGTGATCgt       >  1:166523/1‑149 (MQ=255)
                                                                                                                                                    cACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTGCa  <  1:2050474/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTGCA  >  NC_002947/337922‑338219

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