Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,591,139 C→A 100% 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,1390CA97.7% 560.7 / 3.6 188D66Y (GAT→TAT) dnaTDNA biosynthesis protein (primosomal protein I)
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (97/84);  total (100/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 188 total aligned reads displayed.

TCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTG  >  CP009273/4590994‑4591283
                                                                                                                                                 |                                                                                                                                                
tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGc                                                                                                                                                <  1:2029079/148‑1 (MQ=255)
tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGc                                                                                                                                                <  1:227687/148‑1 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                             <  1:1443111/148‑1 (MQ=255)
   ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTc                                                                                                                                             <  1:3107662/148‑1 (MQ=255)
                 ccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGc                                                                                                                               >  1:548852/1‑148 (MQ=255)
                 ccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGc                                                                                                                               <  1:1879171/148‑1 (MQ=255)
                  cATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                              <  2:2229431/148‑1 (MQ=255)
                     aGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCtttt                                                                                                                           >  1:1457076/1‑148 (MQ=255)
                     aGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCtttt                                                                                                                           <  2:2828066/148‑1 (MQ=255)
                        gccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttctt                                                                                                                        <  1:637704/148‑1 (MQ=255)
                               ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcgc                                                                                                                 >  2:779438/1‑148 (MQ=255)
                               ggCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcgc                                                                                                                 <  1:2764934/148‑1 (MQ=255)
                                    aTAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagca                                                                                                            >  2:242653/1‑148 (MQ=255)
                                       aaaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGtt                                                                                                         >  2:1942278/1‑148 (MQ=255)
                                         aTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCa                                                                                                       >  2:2750482/1‑148 (MQ=255)
                                         aTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCa                                                                                                       >  2:2958829/1‑148 (MQ=255)
                                         aTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCa                                                                                                       >  2:3044916/1‑148 (MQ=255)
                                          tCGGCATCGGGTTGCCAGTCCGGTTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAg                                                                                                      <  1:268675/148‑1 (MQ=255)
                                          tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAg                                                                                                      >  1:3429908/1‑148 (MQ=255)
                                          tCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAg                                                                                                      >  2:1317641/1‑148 (MQ=255)
                                            ggCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGcc                                                                                                    >  1:388718/1‑148 (MQ=255)
                                             gCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                   >  1:2062329/1‑148 (MQ=255)
                                                tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggc                                                                                                <  1:1831914/148‑1 (MQ=255)
                                                tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggc                                                                                                >  1:1470056/1‑148 (MQ=255)
                                                tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggc                                                                                                >  1:2652763/1‑148 (MQ=255)
                                                 cGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcg                                                                                               <  1:3330512/148‑1 (MQ=255)
                                                    gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                            >  1:2714550/1‑148 (MQ=255)
                                                    gTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgc                                                                                            >  1:3363769/1‑148 (MQ=255)
                                                     ttGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgca                                                                                           >  2:1405602/1‑148 (MQ=255)
                                                     ttGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgca                                                                                           >  2:1511891/1‑148 (MQ=255)
                                                        ccAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggc                                                                                        >  2:1111546/1‑148 (MQ=255)
                                                         cAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcg                                                                                       >  2:1611169/1‑148 (MQ=255)
                                                            tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGa                                                                                    >  1:1579768/1‑148 (MQ=255)
                                                            tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGa                                                                                    >  1:428783/1‑148 (MQ=255)
                                                            tCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCCCCAGCTTTTCTTCAAGCGCCACCAGTTCAGCCACGCGTGCACGCGTGa                                                                                    >  1:2220757/1‑148 (MQ=255)
                                                             ccGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGAc                                                                                   >  2:985267/1‑148 (MQ=255)
                                                                 aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCa                                                                               >  2:2775461/1‑148 (MQ=255)
                                                                   aCATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATa                                                                             >  1:2169667/1‑148 (MQ=255)
                                                                    cATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaa                                                                            <  1:3325207/148‑1 (MQ=255)
                                                                    cATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaa                                                                            <  1:3325175/148‑1 (MQ=255)
                                                                     aTGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaa                                                                           <  2:918419/148‑1 (MQ=255)
                                                                      tGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaaa                                                                          >  1:1152862/1‑148 (MQ=255)
                                                                       ggCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAAc                                                                         >  1:3485986/1‑148 (MQ=255)
                                                                        gCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACg                                                                        >  2:586783/1‑148 (MQ=255)
                                                                        gCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACg                                                                        >  1:2761381/1‑148 (MQ=255)
                                                                         cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGc                                                                       <  2:595292/148‑1 (MQ=255)
                                                                         cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGc                                                                       >  2:2518745/1‑148 (MQ=255)
                                                                            aTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgg                                                                    >  2:1070644/1‑148 (MQ=255)
                                                                            aTTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgg                                                                    <  2:2263617/148‑1 (MQ=255)
                                                                              ttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgggg                                                                  >  2:2289164/1‑148 (MQ=255)
                                                                               tCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGc                                                                 <  2:3475289/148‑1 (MQ=255)
                                                                               tCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGc                                                                 <  2:3474988/148‑1 (MQ=255)
                                                                                ccccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCTACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCa                                                                >  1:3533988/1‑148 (MQ=255)
                                                                                     tGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGtt                                                           <  1:526751/148‑1 (MQ=255)
                                                                                     tGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGtt                                                           >  2:3242242/1‑148 (MQ=255)
                                                                                        gTACAGCAACGGGAGCGGCTTGCGGTTCCGGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCAAGCGCCAGCAGTTCAGCCAGGAGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGc                                                        >  1:393643/1‑148 (MQ=255)
                                                                                         tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCa                                                       >  2:472531/1‑148 (MQ=255)
                                                                                           cAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCaaa                                                     <  1:1853637/148‑1 (MQ=255)
                                                                                           cAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCaaa                                                     <  1:2433267/148‑1 (MQ=255)
                                                                                            aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                    <  2:1632939/148‑1 (MQ=255)
                                                                                              cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                  >  2:2636513/1‑148 (MQ=255)
                                                                                                 cGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGc                                               >  2:213402/1‑148 (MQ=255)
                                                                                                        ggCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcaca                                        <  2:1907331/148‑1 (MQ=255)
                                                                                                            tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                    <  1:2570681/148‑1 (MQ=255)
                                                                                                            tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                    >  2:1464940/1‑148 (MQ=255)
                                                                                                            tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGc                                    >  2:1465623/1‑148 (MQ=255)
                                                                                                             gCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGcc                                   >  1:293711/1‑148 (MQ=255)
                                                                                                              cGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCt                                  >  1:1453829/1‑148 (MQ=255)
                                                                                                              cGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCt                                  >  1:2373322/1‑148 (MQ=255)
                                                                                                               ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACTGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCtt                                 <  2:495528/148‑1 (MQ=255)
                                                                                                                 ttCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCa                               >  2:302029/1‑148 (MQ=255)
                                                                                                                  tCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAg                              >  1:528005/1‑148 (MQ=255)
                                                                                                                  tCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAg                              <  1:3036154/148‑1 (MQ=255)
                                                                                                                   ccTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGc                             <  1:693565/148‑1 (MQ=255)
                                                                                                                      gATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                          >  2:2572356/1‑148 (MQ=255)
                                                                                                                      gATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                          >  2:2572452/1‑148 (MQ=255)
                                                                                                                       aTAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtttt                         >  1:67293/1‑148 (MQ=255)
                                                                                                                          gagTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGcc                      >  2:1161930/1‑148 (MQ=255)
                                                                                                                           agTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCa                     >  1:3555590/1‑148 (MQ=255)
                                                                                                                           agTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCa                     >  2:393160/1‑148 (MQ=255)
                                                                                                                            gTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAg                    <  1:1779272/148‑1 (MQ=255)
                                                                                                                             ttGATCGTCCAGAGCTACATAGCTTCCCGGACGCGCCAGATTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGa                   <  1:2711001/148‑1 (MQ=255)
                                                                                                                             ttGATCGTCCAGAGCGACATTGCGTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTGAGCCAGGCGTGCAGGCGTGACGGCATATAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGa                   <  2:79542/148‑1 (MQ=255)
                                                                                                                             ttGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGa                   <  1:1840993/148‑1 (MQ=255)
                                                                                                                              tGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGaa                  >  2:2503718/1‑148 (MQ=255)
                                                                                                                              tGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGaa                  <  1:994972/148‑1 (MQ=255)
                                                                                                                              tGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGaa                  >  1:850334/1‑148 (MQ=255)
                                                                                                                              tGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGaa                  >  1:850332/1‑148 (MQ=255)
                                                                                                                                 tCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACgg               >  1:3333805/1‑148 (MQ=255)
                                                                                                                                       agagCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGt         <  1:96902/148‑1 (MQ=255)
                                                                                                                                       agagCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGt         <  2:921446/148‑1 (MQ=255)
                                                                                                                                        gagCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg        >  2:2761028/1‑148 (MQ=255)
                                                                                                                                        gagCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg        <  2:2589445/148‑1 (MQ=255)
                                                                                                                                          gCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGAt      >  1:371103/1‑148 (MQ=255)
                                                                                                                                             aCATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt   <  1:2767955/148‑1 (MQ=255)
                                                                                                                                             aCATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt   >  2:2144369/1‑148 (MQ=255)
                                                                                                                                             aCATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt   >  1:2420207/1‑148 (MQ=255)
                                                                                                                                              cATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAGACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg  <  2:864670/148‑1 (MQ=255)
                                                                                                                                              cATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg  <  2:3142316/148‑1 (MQ=255)
                                                                                                                                              cATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg  >  1:3219211/1‑148 (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: ‑