Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,645 A→T 13.8% L25Q (CTG→CAG)  sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,150,6450AT13.8% 793.2 / 96.7 283L25Q (CTG→CAG) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base A (100/144);  new base T (21/18);  total (121/162)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 285 total aligned reads displayed.

GGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCC  >  CP009273/4150496‑4150794
                                                                                                                                                     |                                                                                                                                                     
ggTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCa                                                                                                                                                       <  1:1823058/150‑1 (MQ=255)
 gTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAg                                                                                                                                                      >  2:4627782/1‑150 (MQ=255)
 gTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAg                                                                                                                                                      <  1:2110923/150‑1 (MQ=255)
  tCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAgg                                                                                                                                                     <  1:1091179/150‑1 (MQ=255)
   cGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGc                                                                                                                                                    >  1:4591187/1‑150 (MQ=255)
    gCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTTACCAGGcc                                                                                                                                                   <  2:2841317/150‑1 (MQ=255)
        tAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCAtt                                                                                                                                               <  1:3792860/150‑1 (MQ=255)
        tAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCAtt                                                                                                                                               >  2:1861684/1‑150 (MQ=255)
           aGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCa                                                                                                                                            >  2:4745134/1‑150 (MQ=255)
            gTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAg                                                                                                                                           >  2:2505927/1‑150 (MQ=255)
             tGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAgc                                                                                                                                          <  1:3344668/150‑1 (MQ=255)
             tGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAgc                                                                                                                                          >  1:2177807/1‑150 (MQ=255)
                     gattgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTcgc                                                                                                                                  >  2:4752727/1‑150 (MQ=255)
                     gattgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgc                                                                                                                                  >  1:3363188/1‑150 (MQ=255)
                      attgaATTCTATACTGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgcc                                                                                                                                 <  2:4691240/150‑1 (MQ=255)
                      attgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgcc                                                                                                                                 >  1:2498146/1‑150 (MQ=255)
                      attgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgcc                                                                                                                                 <  1:914138/150‑1 (MQ=255)
                        tgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgccgc                                                                                                                               >  2:2380312/1‑150 (MQ=255)
                        tgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgccgc                                                                                                                               <  1:670487/150‑1 (MQ=255)
                              cTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGc                                                                                                                         <  2:428496/150‑1 (MQ=255)
                                ataATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCg                                                                                                                       >  2:2746656/1‑150 (MQ=255)
                                    tGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAAc                                                                                                                   >  2:1346613/1‑150 (MQ=255)
                                      cGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAAcca                                                                                                                 <  2:5046748/150‑2 (MQ=255)
                                          tGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTa                                                                                                             >  2:3509577/1‑150 (MQ=255)
                                           taCGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                            >  2:5118030/2‑150 (MQ=255)
                                           gACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                            >  2:298620/1‑150 (MQ=255)
                                           gACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                            >  1:2313953/1‑150 (MQ=255)
                                           gACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTAc                                                                                                            <  1:1467105/150‑1 (MQ=255)
                                            aCGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACt                                                                                                           <  1:1967435/150‑1 (MQ=255)
                                             cGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                          <  1:524045/150‑1 (MQ=255)
                                              ggCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTAt                                                                                                         <  1:1216097/150‑1 (MQ=255)
                                                  tGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCACCGGGGCCGGAACCTATTACTATGGCa                                                                                                     >  1:3554685/1‑150 (MQ=255)
                                                      ggAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgt                                                                                                 >  2:4879507/1‑150 (MQ=255)
                                                      ggAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgt                                                                                                 >  2:2461462/1‑150 (MQ=255)
                                                        agagCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgtaa                                                                                               <  2:1476875/150‑1 (MQ=255)
                                                         gagCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgtaat                                                                                              <  1:2971925/150‑1 (MQ=255)
                                                                cGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGa                                                                                       >  1:3036450/1‑150 (MQ=255)
                                                                    gggATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAtgt                                                                                   >  1:5088169/1‑150 (MQ=255)
                                                                    gggATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAtgt                                                                                   <  1:1746495/150‑1 (MQ=255)
                                                                     ggATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAtgtg                                                                                  >  1:3283388/1‑150 (MQ=255)
                                                                        tGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCa                                                                               <  1:3324444/150‑1 (MQ=255)
                                                                          gTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAAACTGGCCCATTGCAGCGCCTTCGACGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATg                                                                             >  1:2422865/1‑150 (MQ=255)
                                                                           tGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATgg                                                                            >  1:3688550/1‑150 (MQ=255)
                                                                           tGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATgg                                                                            <  2:2177807/150‑1 (MQ=255)
                                                                            gCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGt                                                                           <  1:3454466/150‑1 (MQ=255)
                                                                             ccccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTa                                                                          <  2:3781129/150‑1 (MQ=255)
                                                                              cccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAg                                                                         <  1:1057735/150‑1 (MQ=255)
                                                                              cccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAg                                                                         <  1:3743860/150‑1 (MQ=255)
                                                                               ccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAgg                                                                        <  1:4361638/150‑1 (MQ=255)
                                                                                 aGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                      >  2:4085989/1‑150 (MQ=255)
                                                                                  gTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCt                                                                     <  2:1293845/150‑1 (MQ=255)
                                                                                   tGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCtt                                                                    >  2:367387/1‑150 (MQ=255)
                                                                                   tGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCtt                                                                    <  2:4540511/150‑1 (MQ=255)
                                                                                    gggTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTa                                                                   <  1:3127629/150‑1 (MQ=255)
                                                                                     ggTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAc                                                                  >  2:2088026/1‑150 (MQ=255)
                                                                                       tGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCt                                                                <  2:1870727/150‑1 (MQ=255)
                                                                                       tGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCt                                                                <  1:3379234/150‑1 (MQ=255)
                                                                                        gCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTg                                                               <  2:1106135/150‑1 (MQ=255)
                                                                                              gccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTta                                                         >  1:2236840/1‑150 (MQ=255)
                                                                                               ccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTtat                                                        <  2:479672/150‑1 (MQ=255)
                                                                                               ccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTtat                                                        >  1:2830782/1‑150 (MQ=255)
                                                                                                cgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCCTGGTAGGGCTTACCTGTTCTgata                                                       >  1:3862397/1‑150 (MQ=255)
                                                                                                  acAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGACGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACa                                                     >  2:4516364/2‑150 (MQ=255)
                                                                                                   cAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACAt                                                    >  2:4913498/1‑150 (MQ=255)
                                                                                                   cAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACAt                                                    >  2:3327709/1‑150 (MQ=255)
                                                                                                       aTTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATaaaa                                                >  2:4414854/1‑150 (MQ=255)
                                                                                                       aTTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATaaaa                                                <  1:1309412/150‑1 (MQ=255)
                                                                                                         tttGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGc                                              <  2:57900/150‑1 (MQ=255)
                                                                                                           tGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCaa                                            >  2:758141/1‑150 (MQ=255)
                                                                                                           tGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCaa                                            >  2:3390987/1‑150 (MQ=255)
                                                                                                              tAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAg                                         <  2:2577275/150‑1 (MQ=255)
                                                                                                              tAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAg                                         <  1:4212360/150‑1 (MQ=255)
                                                                                                              tAACGCTCGATAACTGCGACGCGCGCACCTTGCTCAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAg                                         >  2:4181826/1‑150 (MQ=255)
                                                                                                                aCGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGaa                                       <  1:4087212/150‑1 (MQ=255)
                                                                                                                 cGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAAt                                      <  1:2699293/150‑1 (MQ=255)
                                                                                                                  gCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATg                                     <  1:4145960/150‑1 (MQ=255)
                                                                                                                  gCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATg                                     <  1:1125903/150‑1 (MQ=255)
                                                                                                                    tCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGt                                   >  1:1627539/1‑150 (MQ=255)
                                                                                                                    tCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGACTGGt                                   >  1:2617674/1‑150 (MQ=255)
                                                                                                                      gATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTaa                                 >  1:4504804/1‑150 (MQ=255)
                                                                                                                      gATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTaa                                 <  1:4464987/150‑1 (MQ=255)
                                                                                                                          aCTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACAtt                             <  1:4882981/150‑1 (MQ=255)
                                                                                                                            tGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACAtttt                           <  1:1861684/150‑1 (MQ=255)
                                                                                                                                aCGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATcg                       <  2:1688665/150‑1 (MQ=255)
                                                                                                                                aCGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATcg                       <  1:4898716/150‑1 (MQ=255)
                                                                                                                                   cgcgcACCTTGCTTAACCAGGCCCATTGCAGCTCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCgg                    >  1:4961615/1‑150 (MQ=255)
                                                                                                                                    gcgcACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCggg                   <  1:4159896/150‑1 (MQ=255)
                                                                                                                                      gcACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTa                 <  2:847399/150‑1 (MQ=255)
                                                                                                                                       cACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTaa                <  2:2991949/150‑1 (MQ=255)
                                                                                                                                        aCCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAg               >  2:3935505/1‑150 (MQ=255)
                                                                                                                                        aCCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAg               <  1:2003490/150‑1 (MQ=255)
                                                                                                                                         ccTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGc              <  2:2719491/150‑1 (MQ=255)
                                                                                                                                          cTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGcc             <  1:1394163/150‑1 (MQ=255)
                                                                                                                                          cTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGcc             <  2:2321298/150‑1 (MQ=255)
                                                                                                                                           ttGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCa            <  1:984564/150‑1 (MQ=255)
                                                                                                                                           ttGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCa            >  2:1607598/1‑150 (MQ=255)
                                                                                                                                           ttGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCa            <  1:4627782/150‑1 (MQ=255)
                                                                                                                                             gCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAAt          >  2:3275494/1‑150 (MQ=255)
                                                                                                                                              cTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAAtt         <  1:2205643/150‑1 (MQ=255)
                                                                                                                                                 aaCCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGAt      <  2:2809214/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCC  >  CP009273/4150496‑4150794

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: