Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,591,139 C→A 89.2% D66Y (GAT→TAT)  dnaT ← DNA biosynthesis protein (primosomal protein I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,591,1390CA89.2% 684.4 / 80.8 257D66Y (GAT→TAT) dnaTDNA biosynthesis protein (primosomal protein I)
Reads supporting (aligned to +/- strand):  ref base C (8/20);  new base A (114/115);  total (122/135)
Fisher's exact test for biased strand distribution p-value = 4.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

CTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGT  >  CP009273/4590993‑4591284
                                                                                                                                                  |                                                                                                                                                 
ctctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCt                                                                                                                                                <  1:486073/150‑1 (MQ=255)
 tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCtt                                                                                                                                               >  2:683664/1‑150 (MQ=255)
  ctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                              >  1:3203026/1‑150 (MQ=255)
   tctTAGCTCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTcc                                                                                                                                             >  1:2779048/1‑150 (MQ=255)
    ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTccc                                                                                                                                            <  1:394123/150‑1 (MQ=255)
    ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTccc                                                                                                                                            <  2:3143924/150‑1 (MQ=255)
     ttAGCGCCACGCCCCATAGCGCCGCCAGGCGGATATATTCGGCATCGGGTTGCCAGGCCGGATACATGGCGAATGTCCCCAGGGGTCCAGCATCGGGAGCGGGTTGCGGTTCCTTATAGAGTTGTTCGTCCAGAGCGACATAGCTTCCCg                                                                                                                                           <  2:2037782/150‑1 (MQ=255)
        gcgcCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGAc                                                                                                                                        <  2:1060424/150‑1 (MQ=255)
         cgcCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGAcg                                                                                                                                       >  1:702301/1‑150 (MQ=255)
         cgcCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGAcg                                                                                                                                       <  1:2857198/150‑1 (MQ=255)
              cGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCa                                                                                                                                  <  2:446535/150‑1 (MQ=255)
                  ccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCtt                                                                                                                              <  2:3681035/150‑1 (MQ=255)
                   cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCttt                                                                                                                             <  1:4073287/150‑1 (MQ=255)
                   cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCttt                                                                                                                             <  2:843563/150‑1 (MQ=255)
                   cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCttt                                                                                                                             <  2:3577324/150‑1 (MQ=255)
                    aTAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCtttt                                                                                                                            <  1:3913451/150‑1 (MQ=255)
                    aTAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCtttt                                                                                                                            <  1:1210673/150‑1 (MQ=255)
                     tAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttc                                                                                                                           >  2:4414950/1‑150 (MQ=255)
                           cgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAg                                                                                                                     >  2:2747282/1‑150 (MQ=255)
                            gccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgc                                                                                                                    <  1:2008687/150‑1 (MQ=255)
                             ccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcg                                                                                                                   >  2:1866061/1‑150 (MQ=255)
                              cAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcgc                                                                                                                  >  1:4582015/1‑150 (MQ=255)
                               aGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGcc                                                                                                                 <  1:72996/150‑1 (MQ=255)
                                ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCca                                                                                                                >  2:3989572/1‑150 (MQ=255)
                                  cGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagc                                                                                                              <  1:147848/150‑1 (MQ=255)
                                   ggATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagca                                                                                                             >  2:4622317/1‑150 (MQ=255)
                                    gATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagcag                                                                                                            <  1:3395447/150‑1 (MQ=255)
                                      tAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGtt                                                                                                          <  2:4079056/150‑1 (MQ=255)
                                      tAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGtt                                                                                                          <  1:3122300/150‑1 (MQ=255)
                                         aaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAg                                                                                                       >  2:1392185/1‑150 (MQ=255)
                                         aaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAg                                                                                                       >  1:1120244/1‑150 (MQ=255)
                                          aTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGc                                                                                                      >  1:908679/1‑150 (MQ=255)
                                           tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGcc                                                                                                     >  2:3551632/1‑150 (MQ=255)
                                            cGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                    >  1:1534321/1‑150 (MQ=255)
                                            cGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                    >  2:3258660/1‑150 (MQ=255)
                                            cGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                    >  1:2333456/1‑150 (MQ=255)
                                             ggCATCGGGTTGCCAGTCCGGATTCATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcag                                                                                                   <  2:1795015/150‑1 (MQ=255)
                                                  cGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtg                                                                                              <  1:3067797/150‑1 (MQ=255)
                                                  cGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtg                                                                                              >  1:981846/1‑150 (MQ=255)
                                                  cGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtg                                                                                              <  1:3477381/150‑1 (MQ=255)
                                                   gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                             <  2:1183861/150‑1 (MQ=255)
                                                   gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                             <  1:1935870/150‑1 (MQ=255)
                                                     gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcag                                                                                           <  2:3883151/150‑1 (MQ=255)
                                                          cAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcgtg                                                                                      <  1:3998974/150‑1 (MQ=255)
                                                            gTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGAc                                                                                    >  1:2680248/1‑150 (MQ=255)
                                                            gTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGAc                                                                                    >  2:305859/1‑150 (MQ=255)
                                                             tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACg                                                                                   <  1:2609638/150‑1 (MQ=255)
                                                             tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACg                                                                                   >  1:2577659/1‑150 (MQ=255)
                                                              ccGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACgg                                                                                  >  1:3170736/1‑150 (MQ=255)
                                                              ccGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACgg                                                                                  <  1:1635739/150‑1 (MQ=255)
                                                               cGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGc                                                                                 >  2:4181297/1‑150 (MQ=255)
                                                                ggATACATGGCGAATTTCCCCATGGGTACANCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCa                                                                                >  1:4275047/1‑150 (MQ=255)
                                                                ggATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCa                                                                                >  2:1562505/1‑150 (MQ=255)
                                                                  agaCATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATa                                                                              <  2:4087665/148‑1 (MQ=255)
                                                                  aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATa                                                                              >  1:2684737/1‑150 (MQ=255)
                                                                     cATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaaa                                                                           <  1:2457911/150‑1 (MQ=255)
                                                                     cATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaaa                                                                           >  2:795483/1‑150 (MQ=255)
                                                                            aaTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCggg                                                                    >  2:746324/1‑150 (MQ=255)
                                                                            aaTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCggg                                                                    >  2:2498930/1‑150 (MQ=255)
                                                                             aTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgggg                                                                   <  2:1227446/150‑1 (MQ=255)
                                                                             aTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgggg                                                                   <  1:1741230/150‑1 (MQ=255)
                                                                               ttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCa                                                                 <  1:482598/150‑1 (MQ=255)
                                                                                tCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCAt                                                                >  2:2179783/1‑150 (MQ=255)
                                                                                  cccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTg                                                              <  1:653357/150‑1 (MQ=255)
                                                                                     aTGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTa                                                           >  1:1714924/1‑150 (MQ=255)
                                                                                     aTGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTa                                                           <  2:2584433/150‑1 (MQ=255)
                                                                                           aCAGCAACGGGAGCTGCTTGCGGTTCCTGATAGAGTTGATCGTGCATAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGAGCCAGCATTTCTTACAGGAGTGATTGAGTGACGGCATAACAAGCATGGGCATTGTTAGCAAAt                                                     <  2:2678161/150‑1 (MQ=255)
                                                                                           aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                     <  1:678876/150‑1 (MQ=255)
                                                                                           aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                     <  1:678977/150‑1 (MQ=255)
                                                                                           aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                     <  1:3523468/150‑1 (MQ=255)
                                                                                            cAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATa                                                    >  2:2886510/1‑150 (MQ=255)
                                                                                             aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                   <  2:2004461/150‑1 (MQ=255)
                                                                                             aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                   >  2:2352813/1‑150 (MQ=255)
                                                                                             aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                   <  2:4318848/150‑1 (MQ=255)
                                                                                             aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                   >  2:1675656/1‑150 (MQ=255)
                                                                                              gCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACg                                                  >  1:4559925/1‑150 (MQ=255)
                                                                                               cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACgg                                                 >  1:3808743/1‑150 (MQ=255)
                                                                                               cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACgg                                                 >  1:785662/1‑150 (MQ=255)
                                                                                                 aCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCa                                               <  1:3277188/150‑1 (MQ=255)
                                                                                                  cGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCaa                                              <  1:1245656/150‑1 (MQ=255)
                                                                                                   gggAGCGGCTTGCGGTTCCTGATAGAGTTTATCGTCCAGAGCGACATAGCTTCCCTGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAAc                                             <  2:668201/150‑1 (MQ=255)
                                                                                                    ggAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAAcc                                            <  2:213870/150‑1 (MQ=255)
                                                                                                      aGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcac                                          <  2:3165192/150‑1 (MQ=255)
                                                                                                           cTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                     >  2:2562300/1‑150 (MQ=255)
                                                                                                            ttGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGcc                                    >  2:4633999/1‑150 (MQ=255)
                                                                                                               cGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                 >  2:2300211/1‑150 (MQ=255)
                                                                                                                ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCa                                >  2:1124342/1‑150 (MQ=255)
                                                                                                                 gTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAg                               <  2:2092870/150‑1 (MQ=255)
                                                                                                                      tGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtttt                          >  2:1298765/1‑150 (MQ=255)
                                                                                                                          agagTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCa                      >  2:3913412/1‑150 (MQ=255)
                                                                                                                              ttGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAAc                  <  1:4098667/150‑1 (MQ=255)
                                                                                                                              ttGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTACAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAAc                  <  1:3082789/150‑1 (MQ=255)
                                                                                                                                    gTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg            >  1:106678/1‑150 (MQ=255)
                                                                                                                                    gTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg            <  2:63311/150‑1 (MQ=255)
                                                                                                                                    gTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg            >  1:2359137/1‑150 (MQ=255)
                                                                                                                                    gTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg            >  1:3675094/1‑150 (MQ=255)
                                                                                                                                      ccAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGt          >  2:3474424/1‑150 (MQ=255)
                                                                                                                                       cAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg         <  1:4486863/150‑1 (MQ=255)
                                                                                                                                       cAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg         <  1:4486776/150‑1 (MQ=255)
                                                                                                                                            cGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt    >  2:271152/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGT  >  CP009273/4590993‑4591284

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: