Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,274,881 C→G 5.0% A199G (GCC→GGC)  eamA → O‑acetylserine/cysteine export transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,274,8810CG 5.0% 381.7 / 8.0 119A199G (GCC→GGC) eamAO‑acetylserine/cysteine export transporter
Reads supporting (aligned to +/- strand):  ref base C (65/48);  new base G (3/3);  total (68/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 120 total aligned reads displayed.

CGGTGCGATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGTGAT  >  NC_002947/3274732‑3275028
                                                                                                                                                      |                                                                                                                                                   
cGGTGCGATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGCAACGTGCCGATTGCATTGGCCATGc                                                                                                                                                     <  2:1808014/150‑1 (MQ=255)
      gATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGc                                                                                                                                               >  1:2188934/1‑150 (MQ=255)
       aTGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGcc                                                                                                                                              >  2:920644/1‑150 (MQ=255)
       aTGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGcc                                                                                                                                              <  1:490398/150‑1 (MQ=255)
          tgGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAat                                                                                                                                           >  2:891292/1‑150 (MQ=255)
           ggggCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAata                                                                                                                                          >  1:2851151/1‑150 (MQ=255)
             ggCGTGCGCAACAGCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATc                                                                                                                                        <  2:390642/150‑1 (MQ=255)
              gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGTATTGGCCAAGGCATGGCCAatgc                                                                                                                                         <  1:714234/148‑3 (MQ=255)
              gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGTATTGGCCAAGGCATGGCCAatgc                                                                                                                                         >  2:714234/1‑146 (MQ=255)
              gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCAAGGCATGGCCAatgcgc                                                                                                                                       <  1:2856567/150‑5 (MQ=255)
              gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCAAGGCATGGCCAat                                                                                                                                           <  1:2370240/146‑1 (MQ=255)
              gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCAAGGCATGGCCAat                                                                                                                                           >  2:2370240/1‑146 (MQ=255)
              gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCATGGCACTGTCGTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCa                                                                                                                                       >  2:482198/1‑150 (MQ=255)
                cTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGGTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGt                                                                                                                                     <  2:2011147/150‑1 (MQ=255)
                    gCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGa                                                                                                                                 >  2:3040789/1‑150 (MQ=255)
                    gCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGAAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGa                                                                                                                                 >  1:1634097/1‑150 (MQ=255)
                      aaCATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAg                                                                                                                                      >  1:1663526/1‑143 (MQ=255)
                      aaCATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAg                                                                                                                                      <  2:1663526/143‑1 (MQ=255)
                         ttcatTATCCACCAGTTCGGGTCATGCGACCTTGGCGCGTTGTTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCg                                                                                                                            <  2:2188934/149‑1 (MQ=255)
                           catTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGt                                                                                                                          <  1:891292/150‑1 (MQ=255)
                           catTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGt                                                                                                                          <  2:2781303/150‑1 (MQ=255)
                            atTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTGGGAGTTCGGTg                                                                                                                         <  2:907783/150‑1 (MQ=255)
                             ttACCCGACGGTTCGGCTCGGTCGACCTNGTCGCGTTNGTGATCTGGNGTGGGCTGATAGCGCAACTTCGATTCCTGGCACTCCAGTAGGTGTTGGTAGGCCCCGAGCGCATTGGCCAGGCCTTGGCCAATATAAGTTTGGGTTCGGGGc                                                                                                                        >  1:993965/1‑150 (MQ=255)
                             ttACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGc                                                                                                                        <  2:2851151/150‑1 (MQ=255)
                                cccGACGGTTGGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTgg                                                                                                                     <  2:1323937/150‑1 (MQ=255)
                                 ccGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGc                                                                                                                    <  2:1794599/150‑1 (MQ=255)
                                  cGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGcc                                                                                                                   >  1:946431/1‑150 (MQ=255)
                                   gACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCCCCGCCCATGGGGCGTTAAGGCCCCGAGCGCCTTGGCCCTGCGATGGCGAATCTCAGTGGGATGTCGGTGCTGGccc                                                                                                                  >  1:2833353/1‑150 (MQ=255)
                                   gACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGccc                                                                                                                  >  2:615816/1‑150 (MQ=255)
                                       gTTCGCGTCGGTCGAGCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATCCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCTTGGc                                                                                                              <  2:1514505/150‑1 (MQ=255)
                                              tCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTGCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTg                                                                                                       >  2:2259085/1‑150 (MQ=255)
                                                  tCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCTACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGa                                                                                                                                 >  2:1548739/1‑120 (MQ=255)
                                                      ccTGGTCGCGTTGGTGATCTGGGGGGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTtgt                                                                                               >  1:1886369/1‑150 (MQ=255)
                                                       cTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTtgtg                                                                                              <  1:919437/150‑1 (MQ=255)
                                                         ggtatccttGCTGATCTGGTGTTTGCTTATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGTCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGc                                                                                            <  2:1709673/143‑1 (MQ=255)
                                                         ggTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGCGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGc                                                                                            >  2:2145544/1‑150 (MQ=255)
                                                          gTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTgg                                                                                                                                                <  1:476705/97‑1 (MQ=255)
                                                          gTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTgg                                                                                                                                                >  2:476705/1‑97 (MQ=255)
                                                           tCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGcca                                                                                          <  2:1634097/150‑1 (MQ=255)
                                                                    tGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCg                                                                                 <  1:769185/150‑1 (MQ=255)
                                                                     gATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGtt                                                                                                                                    >  2:2193571/1‑98 (MQ=255)
                                                                       tCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCaa                                                                                                                                            <  1:738362/88‑1 (MQ=255)
                                                                       tCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCaa                                                                                                                                            >  2:738362/1‑88 (MQ=255)
                                                                       tCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGt                                                                                                                                     <  1:1987160/95‑1 (MQ=255)
                                                                            ggTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGc                                                                         <  2:1263197/150‑1 (MQ=255)
                                                                                  cTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTgg                                                                   <  1:1563039/150‑1 (MQ=255)
                                                                                  cTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTgg                                                                   >  2:867382/1‑150 (MQ=255)
                                                                                   tGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTAttt                                                                                                        >  2:2924639/1‑112 (MQ=255)
                                                                                   tGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTAttt                                                                                                        <  1:2924639/112‑1 (MQ=255)
                                                                                     aTACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGTCCATGCCTTGGCCAATATCAGATGCCAAGGGGCGTTGCCGCTGTGGTAGTTGGTCTTTGTGGCCACCATGCTCGGCTACAACCTGCGgcgc                                                                >  1:1897294/1‑150 (MQ=255)
                                                                                     aTACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATg                                                                                     <  2:117330/129‑1 (MQ=255)
                                                                                     aTACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATg                                                                                     >  1:117330/1‑129 (MQ=255)
                                                                                     aTACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGc                                                                <  2:1799443/150‑1 (MQ=255)
                                                                                         cGCCACTTCCATTCCTGTCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGt                                                            >  1:239438/1‑150 (MQ=255)
                                                                                         cGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGt                                                            >  2:1850123/1‑150 (MQ=255)
                                                                                          gCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGtt                                                           <  1:2080901/150‑1 (MQ=255)
                                                                                           ccACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTATAGCCTGTGGAGCAAGTTg                                                          >  1:2371290/1‑150 (MQ=255)
                                                                                           ccACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTg                                                          >  2:512254/1‑150 (MQ=255)
                                                                                           ccACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTAGTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTg                                                          >  2:2059733/1‑150 (MQ=255)
                                                                                            cACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGc                                                         >  2:2117115/1‑150 (MQ=255)
                                                                                             aCTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCATGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCt                                                        >  2:1422790/1‑150 (MQ=255)
                                                                                              cTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCg                                                                                                                            <  2:2675823/81‑1 (MQ=255)
                                                                                              cTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCg                                                                                                                            >  1:2675823/1‑81 (MQ=255)
                                                                                                tCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCATGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTc                                                     <  1:1422790/150‑1 (MQ=255)
                                                                                                       cTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCtg                                                                      <  2:1805288/126‑1 (MQ=255)
                                                                                                       cTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCtg                                                                      >  1:1805288/1‑126 (MQ=255)
                                                                                                         ggCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCAccc                                            <  1:615816/150‑1 (MQ=255)
                                                                                                         ggCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTAGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTCCGCTGGCCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCAccc                                            >  1:2490463/1‑150 (MQ=255)
                                                                                                           cACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCgg                                          >  1:1906174/1‑150 (MQ=255)
                                                                                                            aCTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCggc                                         >  2:2671338/1‑150 (MQ=255)
                                                                                                              tGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGc                                                                                                     >  1:811783/1‑88 (MQ=255)
                                                                                                              tGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGc                                                                                                     <  2:811783/88‑1 (MQ=255)
                                                                                                              tGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCCCCCggcag                                       >  2:174695/1‑150 (MQ=255)
                                                                                                                  cTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGa                                                                  <  1:2746686/119‑1 (MQ=255)
                                                                                                                  cTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGc                                                                >  2:2013710/1‑121 (MQ=255)
                                                                                                                  cTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGc                                                                <  1:2013710/121‑1 (MQ=255)
                                                                                                                    tggtggCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCaaa                                  >  1:660631/1‑150 (MQ=255)
                                                                                                                           tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAata                                                                                                                                          <  2:354808/39‑1 (MQ=255)
                                                                                                                           tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAata                                                                                                                                          >  1:354808/1‑39 (MQ=255)
                                                                                                                           tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTtgt                                                                                               >  2:2647950/1‑82 (MQ=255)
                                                                                                                           tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCAcc                                             >  2:2008915/1‑132 (MQ=255)
                                                                                                                           tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGgcgc                           >  2:2323800/1‑150 (MQ=255)
                                                                                                                           tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGgcgc                           >  1:3035759/1‑150 (MQ=255)
                                                                                                                             gAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCt                                                                                   >  2:566999/1‑92 (MQ=255)
                                                                                                                              aaGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGt                        <  1:785568/150‑1 (MQ=255)
                                                                                                                                   cccGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTTTGGCCACCATGCGCGGCTACAGCCGGTGGAGCAAGTTGCTGTGTCGTCACCCGTCTGGCAAAGGGGCGCCGttatag                   >  2:1411992/1‑146 (MQ=255)
                                                                                                                                        tgataTTGGCCAAGGCATGGCCAATGTCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCgctgct              >  2:1680487/5‑150 (MQ=255)
                                                                                                                                        gcgcATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCt                                                                       <  1:588109/93‑1 (MQ=255)
                                                                                                                                        gcgcATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCt                                                                       >  2:588109/1‑93 (MQ=255)
                                                                                                                                        gcgcATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCgctgct              >  1:1664683/1‑150 (MQ=255)
                                                                                                                                         cgcATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGc                                                                                                     >  2:2901155/1‑62 (MQ=255)
                                                                                                                                           cATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCtgt                                                                     >  1:2083520/1‑92 (MQ=255)
                                                                                                                                           cATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCtgt                                                                     <  2:2083520/92‑1 (MQ=255)
                                                                                                                                            aTTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTa          >  1:1437602/1‑150 (MQ=255)
                                                                                                                                            aTTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTa          >  1:1184455/1‑150 (MQ=255)
                                                                                                                                             ttGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTAc         <  1:210873/150‑1 (MQ=255)
                                                                                                                                             ttGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTAc         >  1:1251328/1‑150 (MQ=255)
                                                                                                                                                   aTGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGTGa   >  1:2236582/1‑150 (MQ=255)
                                                                                                                                                    tGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGTAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGTGAt  <  1:1728444/150‑1 (MQ=255)
                                                                                                                                                    tGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGTGAt  <  2:1251328/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGGTGCGATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGTGAT  >  NC_002947/3274732‑3275028

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