Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,625,283 T→A 6.7% intergenic (‑156/‑55) rob ← / → creA right oriC‑binding transcriptional activator, AraC family/putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,625,2830TA 6.7% 748.7 / 34.1 240intergenic (‑156/‑55)rob/creAright oriC‑binding transcriptional activator, AraC family/putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (90/134);  new base A (12/4);  total (102/138)
Fisher's exact test for biased strand distribution p-value = 8.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 245 total aligned reads displayed.

TCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATT  >  CP009273/4625137‑4625429
                                                                                                                                                  |                                                                                                                                                  
 cATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGt                                                                                                                                                  <  2:102343/148‑1 (MQ=255)
 cATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGt                                                                                                                                                  <  2:2841665/148‑1 (MQ=255)
 cATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGt                                                                                                                                                  <  2:366617/148‑1 (MQ=255)
  cgccTTTCCCCAACGAGCCCCTGACATCAGGTAATAGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTCAGATTTTGTTCTTAATCACTATGtt                                                                                                                                                 <  1:207825/146‑1 (MQ=255)
  aTCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGtt                                                                                                                                                 <  1:3295979/148‑1 (MQ=255)
   tCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTa                                                                                                                                                >  2:390134/1‑148 (MQ=255)
    ccTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTac                                                                                                                                               >  2:2241585/1‑147 (MQ=255)
      tttCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTActt                                                                                                                                             <  2:2178975/148‑1 (MQ=255)
      tttCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTActt                                                                                                                                             <  2:2520958/148‑1 (MQ=255)
       ttCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTa                                                                                                                                            <  2:2918236/148‑1 (MQ=255)
         caacaaCGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTAcc                                                                                                                                          <  1:1494320/148‑1 (MQ=255)
          aacaacGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCg                                                                                                                                         <  1:2885316/148‑1 (MQ=255)
               cGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACg                                                                                                                                    >  2:2963952/1‑148 (MQ=255)
                 aGCACCTGACATCATGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGaa                                                                                                                                  <  1:2790488/148‑1 (MQ=255)
                  gCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAAc                                                                                                                                 <  2:1149178/148‑1 (MQ=255)
                   cACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACt                                                                                                                                >  1:2257927/1‑148 (MQ=255)
                     ccTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGGCGAACTaa                                                                                                                              <  2:1552561/148‑1 (MQ=255)
                     ccTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTaa                                                                                                                              <  2:1237161/148‑1 (MQ=255)
                      cTGTCATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAAt                                                                                                                             <  2:2744706/148‑1 (MQ=255)
                        gACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTg                                                                                                                           <  1:1190661/148‑1 (MQ=255)
                          cATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCt                                                                                                                         <  1:2633878/148‑1 (MQ=255)
                            tCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCg                                                                                                                       <  1:3138006/148‑1 (MQ=255)
                             cAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCgt                                                                                                                      >  2:1063954/1‑148 (MQ=255)
                                gTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCTTGTa                                                                                                                   >  2:2421506/1‑148 (MQ=255)
                                 tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                                  >  1:1927973/1‑148 (MQ=255)
                                   aTTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATa                                                                                                                >  1:2390667/1‑148 (MQ=255)
                                    ttGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAg                                                                                                               >  1:349264/1‑148 (MQ=255)
                                    ttGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAg                                                                                                               >  1:3218246/1‑148 (MQ=255)
                                       gATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATa                                                                                                            <  2:402396/148‑1 (MQ=255)
                                       gATAATAGCATTTTTTGCTGGTTTAGCATTCAGAGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAAATAATTGATCGTGTAATAGATa                                                                                                            >  2:205763/1‑148 (MQ=255)
                                        ataataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATaa                                                                                                           >  2:726033/1‑148 (MQ=255)
                                        ataataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATaa                                                                                                           >  1:928937/1‑148 (MQ=255)
                                        ataataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATaa                                                                                                           <  1:2355717/148‑1 (MQ=255)
                                          aataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTACTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATaaaa                                                                                                         <  2:1294605/148‑1 (MQ=255)
                                           ataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATaaaaa                                                                                                        <  2:1733801/148‑1 (MQ=255)
                                           ataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAACAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATaaaaa                                                                                                        <  2:607641/148‑1 (MQ=255)
                                            taGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAAt                                                                                                       >  2:1177781/1‑148 (MQ=255)
                                              gCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                                     <  1:1735649/148‑1 (MQ=255)
                                                aTTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACTTAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGAAGAACTAATTGCTCGTGTAATAGATAAAAATGTTa                                                                                                   <  1:897535/148‑1 (MQ=255)
                                                aTTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTa                                                                                                   <  2:3109357/148‑1 (MQ=255)
                                                   ttttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACa                                                                                                <  1:3041020/148‑1 (MQ=255)
                                                   ttttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACa                                                                                                <  1:390134/148‑1 (MQ=255)
                                                     ttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATTTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAat                                                                                              <  2:144565/148‑1 (MQ=255)
                                                      tGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAata                                                                                             <  1:1771674/148‑1 (MQ=255)
                                                       gCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAatat                                                                                            <  1:534779/148‑1 (MQ=255)
                                                          gTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATAAGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGaa                                                                                         <  1:1510075/148‑1 (MQ=255)
                                                           ttttAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGaaa                                                                                        <  2:2340823/148‑1 (MQ=255)
                                                             ttAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATa                                                                                      <  1:1189495/148‑1 (MQ=255)
                                                               aGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGATTGACGGAAATAAAAGTATTGAGATTTTGTACTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACa                                                                                    <  2:2206650/148‑1 (MQ=255)
                                                                 cATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAg                                                                                  <  2:2703984/148‑1 (MQ=255)
                                                                    tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAt                                                                               <  1:833988/148‑1 (MQ=255)
                                                                    tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAt                                                                               <  2:3289792/148‑1 (MQ=255)
                                                                    tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATAAGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAt                                                                               >  1:2584878/1‑148 (MQ=255)
                                                                       gtgtTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTg                                                                            >  1:1665388/1‑148 (MQ=255)
                                                                         gtTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGAt                                                                          >  1:65103/1‑148 (MQ=255)
                                                                         gtTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGAt                                                                          <  2:1973385/148‑1 (MQ=255)
                                                                           ttttttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATcc                                                                        <  2:2698448/148‑1 (MQ=255)
                                                                           ttttttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATcc                                                                        <  2:137076/148‑1 (MQ=255)
                                                                           ttttttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATcc                                                                        <  2:2698535/148‑1 (MQ=255)
                                                                              tttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGt                                                                     <  1:1329009/148‑1 (MQ=255)
                                                                              tttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGt                                                                     <  2:701515/148‑1 (MQ=255)
                                                                               ttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTc                                                                    >  1:1002576/1‑148 (MQ=255)
                                                                                tCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCt                                                                   <  1:704887/148‑1 (MQ=255)
                                                                                 cTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCtt                                                                  >  1:2749732/1‑148 (MQ=255)
                                                                                  ttAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCttt                                                                 <  2:518644/148‑1 (MQ=255)
                                                                                     gtagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAg                                                              <  1:2752244/148‑1 (MQ=255)
                                                                                      tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGc                                                             <  1:1900290/148‑1 (MQ=255)
                                                                                      tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGc                                                             >  1:3170562/1‑148 (MQ=255)
                                                                                      tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATAAGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGc                                                             >  1:419767/1‑148 (MQ=255)
                                                                                       agagTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGcc                                                            <  2:2999430/148‑1 (MQ=255)
                                                                                       agagTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGcc                                                            >  1:2292033/1‑148 (MQ=255)
                                                                                          gTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGa                                                         <  1:2952555/148‑1 (MQ=255)
                                                                                          gTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGa                                                         <  2:494666/148‑1 (MQ=255)
                                                                                                   aGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTgg                                                <  2:1100892/148‑1 (MQ=255)
                                                                                                     ggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATAAGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTgggg                                              >  1:1350904/1‑148 (MQ=255)
                                                                                                       cTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGcc                                            <  2:1755020/148‑1 (MQ=255)
                                                                                                       cTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGcc                                            <  2:1375247/148‑1 (MQ=255)
                                                                                                         tGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCAt                                          <  1:1716916/148‑1 (MQ=255)
                                                                                                               aaaTAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTc                                    <  2:2473313/148‑1 (MQ=255)
                                                                                                                aaTAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCa                                   <  2:2353372/148‑1 (MQ=255)
                                                                                                                  tAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATg                                 >  1:2032275/1‑148 (MQ=255)
                                                                                                                  tAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATg                                 >  2:2816132/1‑148 (MQ=255)
                                                                                                                  tAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATg                                 <  1:848512/148‑1 (MQ=255)
                                                                                                                   aaaaGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGc                                <  2:947354/148‑1 (MQ=255)
                                                                                                                   aaaaGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGc                                >  2:2391241/1‑148 (MQ=255)
                                                                                                                   aaaaGTATTGAGATTTTGTTCTTAATCAATAAGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGc                                >  1:3116381/1‑148 (MQ=255)
                                                                                                                      aGTATTGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTCACAATATGAAATACAAGCATTTGATCATGTCTTTAAGCCTGATAATGCTGGGGCCCTTGGCTCATGCaca                             >  2:1599131/1‑146 (MQ=255)
                                                                                                                          ttGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAagag                         >  2:2582612/1‑148 (MQ=255)
                                                                                                                           tGAGATTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAgaga                        <  1:2403240/148‑1 (MQ=255)
                                                                                                                               aTTTTGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTgg                    <  2:57892/148‑1 (MQ=255)
                                                                                                                                ttttGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGt                   >  1:192071/1‑148 (MQ=255)
                                                                                                                                 tttGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGtt                  <  2:1927718/148‑1 (MQ=255)
                                                                                                                                 tttGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGtt                  <  2:808664/148‑1 (MQ=255)
                                                                                                                                  ttGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTc                 >  1:1636566/1‑148 (MQ=255)
                                                                                                                                  ttGTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTc                 <  2:1658065/148‑1 (MQ=255)
                                                                                                                                    gTTCTTAATCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCgg               <  2:1997636/148‑1 (MQ=255)
                                                                                                                                           aTCAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGacac        >  2:2603317/1‑148 (MQ=255)
                                                                                                                                             cAATATGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCg      >  1:299105/1‑148 (MQ=255)
                                                                                                                                               atatGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTa    <  1:1830372/148‑1 (MQ=255)
                                                                                                                                                 atGTTACTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTAtt  >  2:1349585/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATT  >  CP009273/4625137‑4625429

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: