Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,645 A→T 7.1% 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,6450AT 7.1% 1019.5 / 44.6 333L25Q (CTG→CAG) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base A (141/168);  new base T (11/13);  total (152/181)
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

Only 100 of 348 total aligned reads displayed.

GGTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCC  >  CP009273/4150496‑4150794
                                                                                                                                                     |                                                                                                                                                     
ggTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCa                                                                                                                                                       >  1:2055313/1‑150 (MQ=255)
 gTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAg                                                                                                                                                      <  2:688470/150‑1 (MQ=255)
 gTCGCTGTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAg                                                                                                                                                      <  2:4031061/150‑1 (MQ=255)
       gTAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCATGCTTAACCAGGCCCAt                                                                                                                                                >  2:4539864/1‑150 (MQ=255)
        tAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCAtt                                                                                                                                               <  2:447484/150‑1 (MQ=255)
        tAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCAtt                                                                                                                                               <  1:1717623/150‑1 (MQ=255)
        tAAAGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCAtt                                                                                                                                               >  2:1596245/1‑150 (MQ=255)
           aGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCa                                                                                                                                            >  2:718440/1‑150 (MQ=255)
           aGTGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCa                                                                                                                                            >  2:4228935/1‑150 (MQ=255)
             tGGGTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAgc                                                                                                                                          >  1:812305/1‑150 (MQ=255)
                gTTTTGATTGAATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGcc                                                                                                                                       <  1:955471/150‑1 (MQ=255)
                     gattgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCTATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgc                                                                                                                                  >  1:2317588/1‑150 (MQ=255)
                       ttgaATTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTcgccg                                                                                                                                >  2:3812004/1‑150 (MQ=255)
                          aaTTCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGccgccg                                                                                                                             <  1:274433/150‑1 (MQ=255)
                            ttCTATAATGCGGCGGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCggg                                                                                                                           <  1:1010073/150‑1 (MQ=255)
                             tCTATAATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCgggg                                                                                                                          <  2:5258657/150‑1 (MQ=255)
                                 taATGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCgg                                                                                                                      >  2:2362130/1‑150 (MQ=255)
                                  aaTGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGa                                                                                                                     >  2:3712842/1‑150 (MQ=255)
                                   aTGCGGCTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGaa                                                                                                                    <  1:1267370/150‑1 (MQ=255)
                                         cTGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTAtt                                                                                                              <  2:1336081/150‑1 (MQ=255)
                                          tGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTa                                                                                                             >  1:518737/1‑150 (MQ=255)
                                          tGACGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTa                                                                                                             <  1:4482169/150‑1 (MQ=255)
                                            aCGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACt                                                                                                           <  2:196570/150‑1 (MQ=255)
                                             cGGCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTa                                                                                                          >  1:5271885/1‑150 (MQ=255)
                                              ggCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTAt                                                                                                         <  2:1508952/150‑1 (MQ=255)
                                              ggCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTAt                                                                                                         >  2:989368/1‑150 (MQ=255)
                                               gCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATg                                                                                                        <  1:490408/150‑1 (MQ=255)
                                               gCGTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATg                                                                                                        <  1:3873436/150‑1 (MQ=255)
                                                 gTGACGGAGAGCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGc                                                                                                      >  2:3154959/1‑150 (MQ=255)
                                                         gagCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgtaat                                                                                              <  2:5009583/150‑1 (MQ=255)
                                                         gagCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgtaat                                                                                              <  1:4065603/150‑1 (MQ=255)
                                                         gagCTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgtaat                                                                                              >  2:4294327/1‑150 (MQ=255)
                                                            cTTTCGACGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCatcgtaatcgt                                                                                           >  2:3606663/1‑150 (MQ=255)
                                                              ttCGACGGGAGGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATCACTGCGACGCGCGCACCTGGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAg                                                                                         <  2:2297763/150‑1 (MQ=255)
                                                                  aCGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAt                                                                                     <  2:3124990/150‑1 (MQ=255)
                                                                  aCGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGAGAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAt                                                                                     <  2:43909/150‑1 (MQ=255)
                                                                   cGGGATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAtg                                                                                    <  1:1195980/150‑1 (MQ=255)
                                                                     ggATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAtgtg                                                                                  <  2:2604599/150‑1 (MQ=255)
                                                                     ggATGGTGCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAAtgtg                                                                                  <  1:5001729/150‑1 (MQ=255)
                                                                            gCCCCAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGt                                                                           <  1:2280094/150‑1 (MQ=255)
                                                                              cccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAg                                                                         >  2:3235016/1‑150 (MQ=255)
                                                                               ccAGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAgg                                                                        <  2:2793384/150‑1 (MQ=255)
                                                                                 aGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                      >  2:3968752/1‑150 (MQ=255)
                                                                                 aGTGGGTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGc                                                                      >  1:1159158/1‑150 (MQ=255)
                                                                                   tGGGTGCAACCGCCGCCAACATTTTGATATCGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCtt                                                                    <  2:3407035/150‑1 (MQ=255)
                                                                                     ggTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAc                                                                  >  1:1842354/1‑150 (MQ=255)
                                                                                     ggTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAc                                                                  <  1:3240420/150‑1 (MQ=255)
                                                                                     ggTGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTAc                                                                  <  1:1764986/150‑1 (MQ=255)
                                                                                       tGCAACCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCt                                                                <  1:5044404/150‑1 (MQ=255)
                                                                                          aaCCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGtt                                                             >  1:354997/1‑150 (MQ=255)
                                                                                          aaCCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGtt                                                             >  2:2498043/1‑150 (MQ=255)
                                                                                          aaCCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGtt                                                             >  2:111612/1‑150 (MQ=255)
                                                                                           aCCGCCGCCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTc                                                            <  1:3142616/150‑1 (MQ=255)
                                                                                             cgccgcCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCtt                                                          >  1:2867579/1‑150 (MQ=255)
                                                                                             cgccgcCAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCtt                                                          >  1:628249/1‑150 (MQ=255)
                                                                                              gccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTta                                                         <  2:4172571/150‑1 (MQ=255)
                                                                                              gccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTta                                                         >  2:414887/1‑150 (MQ=255)
                                                                                               ccgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAAGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTtat                                                        >  2:2759580/1‑150 (MQ=255)
                                                                                                cgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTtata                                                       <  2:556113/150‑1 (MQ=255)
                                                                                                cgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTtata                                                       <  1:4593267/150‑1 (MQ=255)
                                                                                                cgccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTtata                                                       >  2:353576/1‑150 (MQ=255)
                                                                                                 gccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATAc                                                      >  1:1092833/1‑150 (MQ=255)
                                                                                                 gccAACATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATAc                                                      >  1:4546590/1‑150 (MQ=255)
                                                                                                    aaCATTTTGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATa                                                   <  1:3842133/150‑1 (MQ=255)
                                                                                                        ttttGATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAg                                               <  2:3202475/150‑1 (MQ=255)
                                                                                                            gATAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAAc                                           >  2:1689279/1‑150 (MQ=255)
                                                                                                             aTAACGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACa                                          >  1:1605363/1‑150 (MQ=255)
                                                                                                               aaCGCTCGATAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGTGCTTACCTGTTCTTATACATAAAAGCAACAGa                                        <  1:3736409/150‑1 (MQ=255)
                                                                                                                  gCTCGATAACTGCGACGCGCGCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATg                                     >  2:4433648/1‑150 (MQ=255)
                                                                                                                       aTAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAAc                                >  1:2455103/1‑150 (MQ=255)
                                                                                                                        tAACTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACa                               >  1:2514155/1‑150 (MQ=255)
                                                                                                                          aCTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACAtt                             <  2:2161430/150‑1 (MQ=255)
                                                                                                                          aCTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACAtt                             <  1:2820291/150‑1 (MQ=255)
                                                                                                                           cTGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACAttt                            <  2:2407549/150‑1 (MQ=255)
                                                                                                                            tGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGACTGGTAACAtttt                           >  1:3317553/1‑150 (MQ=255)
                                                                                                                            tGCGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACAtttt                           <  1:2843635/150‑1 (MQ=255)
                                                                                                                              cGACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTAt                         <  1:2381344/150‑1 (MQ=255)
                                                                                                                               gACGCGCGCACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATc                        >  1:2023201/1‑150 (MQ=255)
                                                                                                                                 cgcgcgCACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATcgc                      >  1:2641675/1‑150 (MQ=255)
                                                                                                                                    gcgcACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCggg                   <  1:709171/150‑1 (MQ=255)
                                                                                                                                    gcgcACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCggg                   >  2:901470/1‑150 (MQ=255)
                                                                                                                                       cACCTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTaa                <  2:2964114/150‑1 (MQ=255)
                                                                                                                                       cACCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTaa                <  1:756811/150‑1 (MQ=255)
                                                                                                                                        aCCTTGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAg               <  2:785424/150‑1 (MQ=255)
                                                                                                                                         ccTTGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGc              >  1:3639092/1‑150 (MQ=255)
                                                                                                                                           ttGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCa            >  2:2268599/1‑150 (MQ=255)
                                                                                                                                            tGCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCaa           <  1:4155/150‑1 (MQ=255)
                                                                                                                                            tGCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCaa           >  2:1022930/1‑150 (MQ=255)
                                                                                                                                             gCTTAACCTGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAAt          <  1:1540536/150‑1 (MQ=255)
                                                                                                                                             gCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAAt          <  1:1359016/150‑1 (MQ=255)
                                                                                                                                             gCTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGAAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTATATCGCGGGTAAGCCAAt          >  2:4727318/1‑150 (MQ=255)
                                                                                                                                              cTTAACCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAAtt         <  1:4228935/150‑1 (MQ=255)
                                                                                                                                                  aCCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATc     >  1:1329585/1‑150 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  >  2:4873155/1‑150 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  >  1:4255856/1‑150 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  <  2:1841789/150‑1 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  <  1:4629232/150‑1 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  <  1:484756/150‑1 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  >  2:2869165/1‑150 (MQ=255)
                                                                                                                                                     aGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATcccc  <  2:1073872/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: