Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,591,139 C→A 87.8% 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,1390CA87.8% 438.6 / 63.3 189D66Y (GAT→TAT) dnaTDNA biosynthesis protein (primosomal protein I)
Reads supporting (aligned to +/- strand):  ref base C (10/13);  new base A (69/96);  total (80/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 189 total aligned reads displayed.

TCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTG  >  CP009273/4590994‑4591283
                                                                                                                                                 |                                                                                                                                                
tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGc                                                                                                                                                <  1:1103075/148‑1 (MQ=255)
tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGc                                                                                                                                                <  1:2759041/148‑1 (MQ=255)
tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGc                                                                                                                                                >  1:2521641/1‑148 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                             <  2:231572/148‑1 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                             <  2:1469739/148‑1 (MQ=255)
     tAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTccc                                                                                                                                           >  1:1214084/1‑148 (MQ=255)
       gcgcCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCgg                                                                                                                                         >  2:1857239/1‑148 (MQ=255)
           cACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGAcgc                                                                                                                                     >  1:892232/1‑148 (MQ=255)
             cGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACgcgc                                                                                                                                   >  2:1040805/1‑148 (MQ=255)
               ccccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCa                                                                                                                                 >  1:837780/1‑148 (MQ=255)
               ccccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCa                                                                                                                                 >  2:964396/1‑148 (MQ=255)
               ccccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCa                                                                                                                                 >  2:2752338/1‑148 (MQ=255)
                 ccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGc                                                                                                                               <  1:1697647/148‑1 (MQ=255)
                 ccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGc                                                                                                                               >  2:230588/1‑148 (MQ=255)
                  cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCt                                                                                                                              <  1:147017/148‑1 (MQ=255)
                  cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCt                                                                                                                              <  1:1818481/148‑1 (MQ=255)
                  cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                              >  2:2888008/1‑148 (MQ=255)
                  cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                              >  2:2457432/1‑148 (MQ=255)
                   aTAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCtt                                                                                                                             <  1:3192492/148‑1 (MQ=255)
                   aTAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCtt                                                                                                                             <  2:502837/148‑1 (MQ=255)
                   aTAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCtt                                                                                                                             <  2:502860/148‑1 (MQ=255)
                    tAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCttt                                                                                                                            <  2:2978705/148‑1 (MQ=255)
                      gcgcCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttc                                                                                                                          >  2:2071814/1‑148 (MQ=255)
                      gcgcCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttc                                                                                                                          >  2:2633641/1‑148 (MQ=255)
                        gccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttctt                                                                                                                        <  1:799935/148‑1 (MQ=255)
                        gccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttctt                                                                                                                        <  1:799978/148‑1 (MQ=255)
                         ccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTttcttc                                                                                                                       <  1:1530921/148‑1 (MQ=255)
                         ccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttcttc                                                                                                                       >  1:2699941/1‑148 (MQ=255)
                         ccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttcttc                                                                                                                       <  1:2176590/148‑1 (MQ=255)
                         ccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttcttc                                                                                                                       <  1:1942083/148‑1 (MQ=255)
                              aGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcg                                                                                                                  >  2:677556/1‑148 (MQ=255)
                               ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATTGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcgc                                                                                                                 >  1:2217569/1‑148 (MQ=255)
                               ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcgc                                                                                                                 <  1:769869/148‑1 (MQ=255)
                               ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcgc                                                                                                                 >  1:617629/1‑148 (MQ=255)
                                gCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGcc                                                                                                                >  2:1124302/1‑148 (MQ=255)
                                gCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGcc                                                                                                                <  1:1268374/148‑1 (MQ=255)
                                    aTAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagca                                                                                                            >  2:2863567/1‑148 (MQ=255)
                                     tAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagcag                                                                                                           <  1:956182/148‑1 (MQ=255)
                                     tAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagcag                                                                                                           <  1:1202833/148‑1 (MQ=255)
                                     tAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagcag                                                                                                           >  1:1299903/1‑148 (MQ=255)
                                      aaaaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGt                                                                                                          <  1:374611/148‑1 (MQ=255)
                                              cATCGGGTTGTCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcag                                                                                                  >  2:288135/1‑148 (MQ=255)
                                               aTCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGGTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcagg                                                                                                 <  2:2783398/148‑1 (MQ=255)
                                                tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggc                                                                                                >  1:2460242/1‑148 (MQ=255)
                                                  gggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgt                                                                                              <  2:2990502/148‑1 (MQ=255)
                                                    gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCGAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                            <  1:281121/148‑1 (MQ=255)
                                                    gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCGAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                            <  1:1070724/148‑1 (MQ=255)
                                                    gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                            <  2:2784865/148‑1 (MQ=255)
                                                        ccAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggc                                                                                        <  2:1997743/148‑1 (MQ=255)
                                                         cAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcg                                                                                       <  2:2477112/148‑1 (MQ=255)
                                                          aGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcgt                                                                                      <  2:2220368/148‑1 (MQ=255)
                                                             ccGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGAc                                                                                   >  1:1231588/1‑148 (MQ=255)
                                                             ccGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGAc                                                                                   >  1:112069/1‑148 (MQ=255)
                                                                     aTGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaa                                                                           <  2:1741898/148‑1 (MQ=255)
                                                                     aTGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaa                                                                           <  2:1129049/148‑1 (MQ=255)
                                                                      tGGCGAATGTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaaa                                                                          <  1:1930066/148‑1 (MQ=255)
                                                                        gCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACg                                                                        >  1:2699575/1‑148 (MQ=255)
                                                                            aTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgg                                                                    >  1:1086644/1‑148 (MQ=255)
                                                                            aTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgg                                                                    >  1:2738641/1‑148 (MQ=255)
                                                                             tttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCggg                                                                   >  2:947721/1‑148 (MQ=255)
                                                                             tttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCggg                                                                   >  2:2135039/1‑148 (MQ=255)
                                                                                 cccATGGGTACAGCACCGGGAGCGGCTGGCGGTTCCTGATAGAGTGGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCAt                                                               <  1:371317/148‑1 (MQ=255)
                                                                                   cATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTg                                                             <  2:1231125/148‑1 (MQ=255)
                                                                                     tGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGtt                                                           >  1:3035382/1‑148 (MQ=255)
                                                                                       ggTACAGCACCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAg                                                         <  1:1829974/148‑1 (MQ=255)
                                                                                         tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCa                                                       >  2:282910/1‑148 (MQ=255)
                                                                                         tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCa                                                       >  1:2886250/1‑148 (MQ=255)
                                                                                          aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCaa                                                      <  1:494986/148‑1 (MQ=255)
                                                                                              cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                  >  2:66080/1‑148 (MQ=255)
                                                                                                      gCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACca                                          <  1:2611877/148‑1 (MQ=255)
                                                                                                       cGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcac                                         <  2:802155/148‑1 (MQ=255)
                                                                                                        ggCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcaca                                        <  2:2480514/148‑1 (MQ=255)
                                                                                                        ggCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcaca                                        <  2:2934691/148‑1 (MQ=255)
                                                                                                         gCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCacac                                       >  2:1922639/1‑148 (MQ=255)
                                                                                                          cTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACAcc                                      >  1:321955/1‑148 (MQ=255)
                                                                                                          cTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGACGGGGCATTGTTAGCAAATACGGCAACCACAcc                                      >  2:3292151/1‑148 (MQ=255)
                                                                                                            tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCCGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                    <  1:939626/148‑1 (MQ=255)
                                                                                                             gCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGcc                                   >  1:1398156/1‑148 (MQ=255)
                                                                                                                gTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCTGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                <  1:1258216/148‑1 (MQ=255)
                                                                                                                  tCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAg                              >  2:795005/1‑148 (MQ=255)
                                                                                                                   ccTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGc                             >  1:482660/1‑148 (MQ=255)
                                                                                                                     tGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtt                           >  2:632558/1‑148 (MQ=255)
                                                                                                                     tGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtt                           <  2:511719/148‑1 (MQ=255)
                                                                                                                      gATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                          >  2:2409955/1‑148 (MQ=255)
                                                                                                                        tAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTg                        <  2:2161879/148‑1 (MQ=255)
                                                                                                                        tAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTg                        >  2:470612/1‑148 (MQ=255)
                                                                                                                        tAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTg                        <  2:2162174/148‑1 (MQ=255)
                                                                                                                         agagTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGc                       <  2:1889315/148‑1 (MQ=255)
                                                                                                                             ttGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGa                   <  1:516885/148‑1 (MQ=255)
                                                                                                                                 tCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACgg               <  1:2863571/148‑1 (MQ=255)
                                                                                                                                  cGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGt              <  1:835622/148‑1 (MQ=255)
                                                                                                                                   gTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGtt             <  1:380634/148‑1 (MQ=255)
                                                                                                                                     ccAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTg           >  2:1486680/1‑148 (MQ=255)
                                                                                                                                         agCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa       >  1:2946529/1‑148 (MQ=255)
                                                                                                                                         agCGACATAGCTTCCCGGACGCGCCCGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa       <  1:667577/148‑1 (MQ=255)
                                                                                                                                         agCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa       >  1:1106806/1‑148 (MQ=255)
                                                                                                                                           cGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATc     >  1:2025936/1‑148 (MQ=255)
                                                                                                                                              cATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg  <  1:140535/148‑1 (MQ=255)
                                                                                                                                              cATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg  <  1:890843/148‑1 (MQ=255)
                                                                                                                                              cATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg  <  1:3233666/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTG  >  CP009273/4590994‑4591283

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: