Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 568,464 G→T 5.5% E30* (GAA→TAA)  ninE → DLP12 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273568,4640GT 5.5% 1228.7 / 31.1 316E30* (GAA→TAA) ninEDLP12 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (126/172);  new base T (10/8);  total (136/180)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

TGCCACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCC  >  CP009273/568315‑568613
                                                                                                                                                     |                                                                                                                                                     
  ccACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGaa                                                                                                                                                     >  2:32339/1‑150 (MQ=255)
   cACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGaaa                                                                                                                                                    <  1:3943570/150‑1 (MQ=255)
     ccTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATa                                                                                                                                                  <  1:4124173/150‑1 (MQ=255)
        gACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCa                                                                                                                                               >  2:2751063/1‑150 (MQ=255)
        gACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCa                                                                                                                                               >  1:1707789/1‑150 (MQ=255)
           tATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAaca                                                                                                                                            >  1:1962642/1‑150 (MQ=255)
            aTTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAacac                                                                                                                                           >  1:3415607/1‑150 (MQ=255)
            aTTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAacac                                                                                                                                           >  1:3264484/1‑150 (MQ=255)
            aTTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAacac                                                                                                                                           >  1:184692/1‑150 (MQ=255)
              tCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTg                                                                                                                                         >  2:1835900/1‑150 (MQ=255)
               ccACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGc                                                                                                                                        >  1:2193696/1‑150 (MQ=255)
               ccACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGc                                                                                                                                        >  1:4734362/1‑150 (MQ=255)
                 acaGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTc                                                                                                                                      >  1:4548222/1‑150 (MQ=255)
                   agagCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCgg                                                                                                                                    >  1:2895546/1‑150 (MQ=255)
                   agagCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCgg                                                                                                                                    <  1:1231528/150‑1 (MQ=255)
                       cTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGatat                                                                                                                                >  1:4311355/1‑150 (MQ=255)
                          cGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATAcc                                                                                                                             <  1:2546351/150‑1 (MQ=255)
                            aGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGc                                                                                                                           <  2:3211094/150‑1 (MQ=255)
                               tCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAg                                                                                                                        <  1:1476158/150‑1 (MQ=255)
                               tCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAg                                                                                                                        <  1:667937/150‑1 (MQ=255)
                                ccAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGc                                                                                                                       >  2:2247856/1‑150 (MQ=255)
                                 cAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCt                                                                                                                      >  2:3122785/1‑150 (MQ=255)
                                 cAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCt                                                                                                                      >  2:2926627/1‑150 (MQ=255)
                                  aGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCtt                                                                                                                     <  1:3364248/150‑1 (MQ=255)
                                     tgctgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGt                                                                                                                  <  1:2405812/150‑1 (MQ=255)
                                        tgcGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGa                                                                                                               >  1:2477113/1‑150 (MQ=255)
                                           gcgACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaa                                                                                                            <  2:1250067/150‑1 (MQ=255)
                                              aCAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaaaaa                                                                                                         >  1:2042977/1‑150 (MQ=255)
                                               cAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATaaaaaa                                                                                                        >  2:50004/1‑150 (MQ=255)
                                                aaGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAAt                                                                                                       >  2:2855516/1‑150 (MQ=255)
                                                 aGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAACTg                                                                                                      >  2:4628576/1‑150 (MQ=255)
                                                 aGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATg                                                                                                      <  1:4475969/150‑1 (MQ=255)
                                                  gCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATgt                                                                                                     <  1:3498386/150‑2 (MQ=255)
                                                   cAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGt                                                                                                    >  1:3701619/1‑150 (MQ=255)
                                                    aggaggCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGtt                                                                                                   >  1:1310072/1‑150 (MQ=255)
                                                              aTGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgca                                                                                         >  2:3742947/1‑150 (MQ=255)
                                                               tGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcag                                                                                        <  2:13074/150‑1 (MQ=255)
                                                                 gCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGgcagca                                                                                      >  1:2864041/1‑150 (MQ=255)
                                                                   tACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAag                                                                                    >  1:4233996/1‑150 (MQ=255)
                                                                   tACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAag                                                                                    <  1:2968457/150‑1 (MQ=255)
                                                                    acacCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAagg                                                                                   >  2:1407367/1‑150 (MQ=255)
                                                                           ttATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTc                                                                            >  1:4531672/1‑150 (MQ=255)
                                                                             aTTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCAt                                                                          >  2:2322107/1‑150 (MQ=255)
                                                                               tCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATgg                                                                        <  2:4488458/150‑1 (MQ=255)
                                                                                 gtgtCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCt                                                                      <  1:375149/150‑1 (MQ=255)
                                                                                      aTGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGAttt                                                                 <  2:2556423/150‑1 (MQ=255)
                                                                                       tGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTg                                                                >  1:1126281/1‑150 (MQ=255)
                                                                                         aaCGGACACATCTACAGAGTACCAACTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTgag                                                              <  2:1851664/150‑1 (MQ=255)
                                                                                         aaCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTgag                                                              <  2:2887428/150‑1 (MQ=255)
                                                                                         aaCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTgag                                                              <  2:682655/150‑1 (MQ=255)
                                                                                         aaCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTgag                                                              <  2:1339367/150‑1 (MQ=255)
                                                                                         aaCGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTgag                                                              >  1:4599811/1‑150 (MQ=255)
                                                                                             gACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGaaaa                                                          <  2:3186458/150‑1 (MQ=255)
                                                                                                acaTCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAagca                                                       <  2:3497587/150‑1 (MQ=255)
                                                                                                        agagTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTgg                                               <  2:1529193/150‑1 (MQ=255)
                                                                                                        agagTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTgg                                               <  2:4577715/150‑1 (MQ=255)
                                                                                                             aCCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCggg                                          <  2:3792220/150‑1 (MQ=255)
                                                                                                             aCCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCggg                                          <  2:2210653/150‑1 (MQ=255)
                                                                                                             aCCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCggg                                          >  1:3879532/1‑150 (MQ=255)
                                                                                                                 aaTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATg                                      <  2:566234/150‑1 (MQ=255)
                                                                                                                  aTCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGc                                     >  2:4673500/1‑150 (MQ=255)
                                                                                                                     gtcgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAg                                  <  2:3269537/150‑1 (MQ=255)
                                                                                                                     gtcgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAg                                  <  2:4273384/150‑1 (MQ=255)
                                                                                                                     gtcgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAg                                  <  1:3570347/150‑1 (MQ=255)
                                                                                                                      tcgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAgg                                 <  2:3605618/150‑1 (MQ=255)
                                                                                                                      tcgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAgg                                 <  1:320167/150‑1 (MQ=255)
                                                                                                                       cgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGt                                >  1:1517715/1‑150 (MQ=255)
                                                                                                                       cgtAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGt                                >  2:805634/1‑150 (MQ=255)
                                                                                                                        gtAAGCGTAAACCTGAGCTTAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCTGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATTCCAGGTa                               <  2:2978557/150‑1 (MQ=255)
                                                                                                                         tAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTaa                              <  2:4363390/150‑1 (MQ=255)
                                                                                                                         tAAGCGTAAACCTGAGCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTaa                              >  1:2115254/1‑150 (MQ=255)
                                                                                                                          aaGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAg                             <  2:4735004/150‑1 (MQ=255)
                                                                                                                           aGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGa                            <  2:1210876/150‑1 (MQ=255)
                                                                                                                            gCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGaa                           <  1:1738231/150‑1 (MQ=255)
                                                                                                                             cGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAAt                          <  2:2797275/150‑1 (MQ=255)
                                                                                                                              gTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATc                         >  2:1365395/1‑150 (MQ=255)
                                                                                                                              gTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATc                         <  2:1143126/150‑1 (MQ=255)
                                                                                                                               tAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATcc                        <  1:503016/150‑1 (MQ=255)
                                                                                                                                 aaCCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCt                      >  2:82113/1‑150 (MQ=255)
                                                                                                                                  aCCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTTGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTg                     <  2:1980803/150‑1 (MQ=255)
                                                                                                                                    cTGAGCTGAAGCCATACGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGc                   <  2:611760/150‑1 (MQ=255)
                                                                                                                                     tGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCg                  <  1:979440/150‑1 (MQ=255)
                                                                                                                                     tGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCg                  >  1:1617308/1‑150 (MQ=255)
                                                                                                                                       aGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTa                <  1:2692759/150‑1 (MQ=255)
                                                                                                                                       aGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTa                <  1:103650/150‑1 (MQ=255)
                                                                                                                                       aGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTa                <  2:4572756/150‑1 (MQ=255)
                                                                                                                                        gCTGAAGCCATCCTAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTAt               <  2:3889055/150‑1 (MQ=255)
                                                                                                                                        gCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTAt               <  1:3122866/150‑1 (MQ=255)
                                                                                                                                          tGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGt             >  2:2533228/1‑150 (MQ=255)
                                                                                                                                             tgCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTact          >  2:2999922/2‑148 (MQ=255)
                                                                                                                                               ccATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATgg        >  2:1699620/1‑150 (MQ=255)
                                                                                                                                                cATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGc       <  2:4548222/150‑1 (MQ=255)
                                                                                                                                                cATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGc       <  2:2479804/150‑1 (MQ=255)
                                                                                                                                                cATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGc       <  1:4085649/150‑1 (MQ=255)
                                                                                                                                                cATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGc       >  2:2117207/1‑150 (MQ=255)
                                                                                                                                                cATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGc       >  2:4204094/1‑150 (MQ=255)
                                                                                                                                                 aTCCGAAAGACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCa      <  2:4704577/150‑1 (MQ=255)
                                                                                                                                                  tCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCaa     >  1:4429802/1‑150 (MQ=255)
                                                                                                                                                   ccGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCTTGGCTGATTTGAGAAAAGCATCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAAc    <  1:3058117/150‑1 (MQ=255)
                                                                                                                                                     gAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAAccc  >  1:2083329/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCACCTGACTATTCCACAGAGCTGCGAGTTCCAGCTGCTGCGCGACAAGCAGGAGGCGTAATGGCTACACCGCTTATTCGTGTCATGAACGGACACATCTACAGAGTACCAAATCGTCGTAAGCGTAAACCTGAGCTGAAGCCATCCGAAATACCAACACTGCTCGGATATACCGCCAGCTTGGTTGATAAAAAATGGTTGCGACTGGCAGCAAGGAGGAGTCATGGCTGATTTGAGAAAAGCAGCGCGTGGTCGGGAATGCCAGGTAAGAATCCCTGGCGTATGTAATGGCAACCC  >  CP009273/568315‑568613

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: