Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,297,107:1 +C 100% coding (538/624 nt) ccmA ← heme export ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,297,1071.C91.5% 1814.7 / 170.5 519coding (538/624 nt)ccmAheme export ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base . (23/21);  new base C (257/218);  total (280/239)
Fisher's exact test for biased strand distribution p-value = 8.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 566 total aligned reads displayed.

GGCGATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGACGCGGCGTTGTTGCCCGGCC  >  NC_000913/2296959‑2297255
                                                                                                                                                     |                                                                                                                                                    
ggCGATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATc                                                                                                                                                     <  1:1741722/150‑1 (MQ=255)
 gCGATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATcc                                                                                                                                                    >  1:10227031/1‑150 (MQ=255)
  cGATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATccc                                                                                                                                                   <  2:9355303/150‑1 (MQ=255)
  cGATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTTCTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATcccc                                                                                                                                                  <  1:11444111/150‑1 (MQ=255)
   gATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGACTCACAATcccc                                                                                                                                                  >  1:10957918/1‑150 (MQ=255)
     tttCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATcccccc                                                                                                                                                <  2:2111273/150‑1 (MQ=255)
         ggCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCTCAATcccccccctg                                                                                                                                            >  1:11928991/1‑147 (MQ=255)
         ggCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATcccccccctg                                                                                                                                            >  2:1258277/1‑147 (MQ=255)
             cTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTcc                                                                                                                                         <  2:5077042/150‑1 (MQ=255)
             cTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTcc                                                                                                                                         <  2:8431554/150‑1 (MQ=255)
              tATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCg                                                                                                                                        >  1:8723216/1‑150 (MQ=255)
              tATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCg                                                                                                                                        >  1:727789/1‑150 (MQ=255)
              tATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCg                                                                                                                                        >  2:9110562/1‑150 (MQ=255)
                tGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTa                                                                                                                                      <  2:2598632/150‑1 (MQ=255)
                tGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTa                                                                                                                                      >  2:8816560/1‑150 (MQ=255)
                      aaCGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTAtgctgc                                                                                                                                >  2:11332544/1‑150 (MQ=255)
                       aCGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTgcg                                                                                                                               <  1:487884/150‑1 (MQ=255)
                        cGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTgcgc                                                                                                                              <  2:10625789/150‑1 (MQ=255)
                         gCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGcc                                                                                                                             <  2:7774358/150‑1 (MQ=255)
                           tACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCAt                                                                                                                           >  1:187227/1‑150 (MQ=255)
                            acacGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATa                                                                                                                          <  1:5942554/150‑1 (MQ=255)
                            acacGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATa                                                                                                                          <  1:5473629/150‑1 (MQ=255)
                              acGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACg                                                                                                                        >  2:184078/1‑150 (MQ=255)
                              acGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACg                                                                                                                        >  2:6469950/1‑150 (MQ=255)
                              acGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACg                                                                                                                        >  2:7607091/1‑150 (MQ=255)
                               cGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGc                                                                                                                       <  2:10720042/150‑1 (MQ=255)
                                gCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCt                                                                                                                      <  1:4348867/150‑1 (MQ=255)
                                gCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCt                                                                                                                      >  2:543849/1‑150 (MQ=255)
                                gCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCt                                                                                                                      <  2:2834939/150‑1 (MQ=255)
                                gCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTg                                                                                                                     >  2:3360545/1‑150 (MQ=255)
                                 cAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTg                                                                                                                     <  1:11619609/150‑1 (MQ=255)
                                 cAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTg                                                                                                                     <  1:11815580/150‑1 (MQ=255)
                                  aGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTgg                                                                                                                    <  2:4033246/150‑1 (MQ=255)
                                   gCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTggg                                                                                                                   <  1:9630489/150‑1 (MQ=255)
                                    cTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGt                                                                                                                  <  1:2712706/150‑1 (MQ=255)
                                    cTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGt                                                                                                                  <  1:10057696/150‑1 (MQ=255)
                                    cTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGt                                                                                                                  <  1:7556063/150‑1 (MQ=255)
                                    cTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGt                                                                                                                  <  1:3576363/150‑1 (MQ=255)
                                     tCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTc                                                                                                                 >  2:5968267/1‑150 (MQ=255)
                                       aaGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGa                                                                                                              <  2:6747910/150‑1 (MQ=255)
                                        aGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGa                                                                                                              <  1:3225146/150‑1 (MQ=255)
                                          aCGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACg                                                                                                            >  2:198507/1‑150 (MQ=255)
                                            ggAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATc                                                                                                         <  1:5170726/150‑1 (MQ=255)
                                             gAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCg                                                                                                        >  1:8634456/1‑150 (MQ=255)
                                                aaTGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGac                                                                                                      <  1:10477466/150‑1 (MQ=255)
                                                aaTGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGac                                                                                                      >  2:11079478/1‑150 (MQ=255)
                                                    cgcCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCg                                                                                                  <  2:5260166/150‑1 (MQ=255)
                                                     gcCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATATGCTGGGTCAGACGATCGACACCGt                                                                                                 >  2:3796182/1‑150 (MQ=255)
                                                       cAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTa                                                                                               >  1:896251/1‑150 (MQ=255)
                                                       cAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCACCATACGCTGGGTCAGACGATCGACACCGTTa                                                                                               <  1:12066083/150‑1 (MQ=255)
                                                         gAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCGGACGATCGACACCGTTAAc                                                                                             >  1:5241502/1‑150 (MQ=255)
                                                         gAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAAc                                                                                             <  2:6201673/150‑1 (MQ=255)
                                                            catcatGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCa                                                                                         >  1:4476248/1‑150 (MQ=255)
                                                             atcatGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCa                                                                                         >  2:1561962/1‑150 (MQ=255)
                                                             atcatGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCa                                                                                         >  2:11917008/1‑150 (MQ=255)
                                                                 tGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATcg                                                                                     >  1:10705784/1‑150 (MQ=255)
                                                                     gCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGt                                                                                 >  2:10258200/1‑150 (MQ=255)
                                                                     gCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGt                                                                                 <  1:2397010/150‑1 (MQ=255)
                                                                      cccTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTa                                                                                >  1:264872/1‑150 (MQ=255)
                                                                      cccTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTa                                                                                >  2:5711867/1‑150 (MQ=255)
                                                                        cTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTaaa                                                                              <  1:5679125/150‑1 (MQ=255)
                                                                          cGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAgg                                                                           <  2:8086464/150‑1 (MQ=255)
                                                                             ttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCt                                                                         <  2:2205774/150‑1 (MQ=255)
                                                                              tGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCtc                                                                        <  2:9756736/150‑1 (MQ=255)
                                                                                    aGTGAAATGCGGCGAATTTTAGTTTCAGCAACGTTGATCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCTATCGCGGTAAAAGGCTCGTCGAg                                                                  >  2:7272598/1‑150 (MQ=255)
                                                                                    aGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAg                                                                  <  2:11789556/150‑1 (MQ=255)
                                                                                        aaaTGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATc                                                              <  2:5498671/150‑1 (MQ=255)
                                                                                          aTGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCa                                                            >  2:10414243/1‑150 (MQ=255)
                                                                                          aTGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCa                                                            >  2:8153638/1‑150 (MQ=255)
                                                                                               gcgAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGGGGATCCATAACg                                                       >  2:4894085/1‑150 (MQ=255)
                                                                                                     tttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCAc                                                 <  1:6332700/150‑1 (MQ=255)
                                                                                                       tACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACgg                                               >  1:3399496/1‑150 (MQ=255)
                                                                                                       tACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACgg                                               <  1:6851449/150‑1 (MQ=255)
                                                                                                          tttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCa                                           <  2:10926559/150‑1 (MQ=255)
                                                                                                          tttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCa                                           >  1:7957970/1‑150 (MQ=255)
                                                                                                          tttCAGCAACGTTGAGCGGCTGGGGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTACCGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCAGAACGTGGCACGGGTc                                            >  1:840942/1‑150 (MQ=255)
                                                                                                           ttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCa                                           <  1:6978018/150‑1 (MQ=255)
                                                                                                           ttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCa                                           <  2:684353/150‑1 (MQ=255)
                                                                                                           ttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCa                                           <  2:11607909/150‑1 (MQ=255)
                                                                                                           ttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCa                                           <  1:9917088/150‑1 (MQ=255)
                                                                                                             cAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACTATCGACACCGTTAACGTCAATCGCTGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGc                                         <  1:1364626/150‑1 (MQ=255)
                                                                                                             cAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGc                                         <  2:6648660/150‑1 (MQ=255)
                                                                                                              aGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCca                                       >  2:4393645/1‑150 (MQ=255)
                                                                                                               gCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCca                                       <  1:3659257/150‑1 (MQ=255)
                                                                                                                cAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCcac                                      >  1:1964307/1‑150 (MQ=255)
                                                                                                                cAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCcac                                      >  1:8235390/1‑150 (MQ=255)
                                                                                                                cAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCcac                                      <  1:7802126/150‑1 (MQ=255)
                                                                                                                   cGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGa                                   >  2:8671932/1‑150 (MQ=255)
                                                                                                                    gTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGAGCGACACCGGGAACGTCAATCGCGGTCCAAGGCTCGTCGAGGAGCCAGCACGTGGCACGGGTCCGCCACCGAc                                  >  2:4680996/1‑150 (MQ=255)
                                                                                                                      tGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGAcgc                                >  1:6579367/1‑150 (MQ=255)
                                                                                                                            gCTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTaaa                          <  1:9251745/150‑1 (MQ=255)
                                                                                                                             cTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGGGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAg                         <  2:6941694/150‑1 (MQ=255)
                                                                                                                             cTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAg                         <  2:12151673/150‑1 (MQ=255)
                                                                                                                             cTGGTGGGTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAg                         >  1:226041/1‑150 (MQ=255)
                                                                                                                                    gTAGTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGAcgcg                  >  1:4716623/1‑150 (MQ=255)
                                                                                                                                       gTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGACgcggcg               >  1:5832928/1‑150 (MQ=255)
                                                                                                                                       gTCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGACgcggcg               >  1:5356855/1‑150 (MQ=255)
                                                                                                                                        tCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGACGCGGCgt              >  2:981123/1‑150 (MQ=255)
                                                                                                                                        tCAGAATCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGAAGCGGCgt              >  2:1904972/1‑150 (MQ=255)
                                                                                                                                             aTCACAATCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGACGCGGCGttgttg         >  1:9560014/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGATTTCGGCGCTATGGCGAAACGCTACACGCAGCTCAAGACGGAAAATGCGCCAGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAATCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCCATAACGTGGCACGGGTCAGCCACAGACGCGCTAAAGCGACGCGGCGTTGTTGCCCGGCC  >  NC_000913/2296959‑2297255

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