Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,837,177 C→T 5.6% intergenic (‑28/‑147) yicO ← / → ade putative adenine permease/cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,837,1770CT 5.6% 1054.4 / 33.7 301intergenic (‑28/‑147)yicO/adeputative adenine permease/cryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (130/154);  new base T (7/10);  total (137/164)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 318 total aligned reads displayed.

AAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAATATG  >  CP009273/3837028‑3837326
                                                                                                                                                     |                                                                                                                                                     
aaaaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATc                                                                                                                                                       <  2:1484135/150‑1 (MQ=255)
aaaaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATc                                                                                                                                                       >  2:68483/1‑150 (MQ=255)
aaaaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATc                                                                                                                                                       >  1:2259760/1‑150 (MQ=255)
aaaaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATc                                                                                                                                                       <  2:3809085/150‑1 (MQ=255)
  aaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATccc                                                                                                                                                     <  1:3834015/150‑1 (MQ=255)
  aaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATccc                                                                                                                                                     <  1:1347728/150‑1 (MQ=255)
  aaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATccc                                                                                                                                                     >  2:2221278/1‑150 (MQ=255)
   aaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCt                                                                                                                                                    >  1:1884457/1‑150 (MQ=255)
   aaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTAATCAATTTTTTGTCCATGGGAGGAGATTTAATCCCt                                                                                                                                                    <  2:2900888/150‑1 (MQ=255)
      gTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCAGTGGAGGAGATTTAATCCCTTGc                                                                                                                                                 >  2:3207020/1‑150 (MQ=255)
       tGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCt                                                                                                                                                <  2:4328258/150‑1 (MQ=255)
          tCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCtttt                                                                                                                                             <  1:3346568/150‑1 (MQ=255)
            aTCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGt                                                                                                                                           <  2:2040156/150‑1 (MQ=255)
             tCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTa                                                                                                                                          >  1:4206416/1‑150 (MQ=255)
               cccGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTaaa                                                                                                                                        >  1:4417452/1‑150 (MQ=255)
                      tcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATg                                                                                                                                 >  2:3133219/1‑150 (MQ=255)
                      tcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATg                                                                                                                                 >  1:2833544/1‑150 (MQ=255)
                      tcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATg                                                                                                                                 <  1:1557398/150‑1 (MQ=255)
                         atTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGaca                                                                                                                              <  1:3417187/150‑1 (MQ=255)
                           tCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGacaac                                                                                                                            >  1:2064913/1‑150 (MQ=255)
                              gtgCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACaacaac                                                                                                                         >  1:180048/1‑150 (MQ=255)
                                gCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACgg                                                                                                                       <  1:2965087/150‑1 (MQ=255)
                                      ggtggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTCTTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGt                                                                                                                 >  2:3790229/1‑150 (MQ=255)
                                      ggtggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGt                                                                                                                 >  1:1374938/1‑150 (MQ=255)
                                      ggtggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGt                                                                                                                 >  2:1530521/1‑150 (MQ=255)
                                         tgtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGAt                                                                                                              <  2:2478306/149‑1 (MQ=255)
                                         ggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTACAATGGAGTAGATTTAATCCCTTGCTTTTTTAAAACAAATGAAAACAACGGTTCAGTGAt                                                                                                              >  1:3663187/1‑150 (MQ=255)
                                             ccATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAAtt                                                                                                          <  2:4269207/150‑1 (MQ=255)
                                              cATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTa                                                                                                         <  1:838458/150‑1 (MQ=255)
                                               aTGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTAt                                                                                                        <  2:1177105/150‑1 (MQ=255)
                                                   tGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATcaca                                                                                                    >  2:1070921/1‑150 (MQ=255)
                                                    gAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACAt                                                                                                   <  2:368470/150‑1 (MQ=255)
                                                      ggTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACAttt                                                                                                 <  2:1476198/150‑1 (MQ=255)
                                                           tttAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTACAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAAtt                                                                                            <  1:4475317/150‑1 (MQ=255)
                                                               aataatCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGcac                                                                                        <  1:2897512/150‑1 (MQ=255)
                                                                  aatCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACAtt                                                                                     <  1:3201337/150‑1 (MQ=255)
                                                                   ttCGCGTATGCGTGACTGAGTCATTACTGACGTAATCGGTATTGTCTTTATTCATTTTTTTGTGCTTTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACTAATGACATCAACGGTTCAGTGATAATTATCACATTTCATTTGAACATTa                                                                                    <  2:1240707/149‑1 (MQ=255)
                                                                     cgcgAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAAt                                                                                  >  2:4356019/1‑150 (MQ=255)
                                                                           aGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGTGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGatat                                                                            >  2:1297887/1‑150 (MQ=255)
                                                                           aGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGatat                                                                            <  1:3164660/150‑1 (MQ=255)
                                                                            gCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTGTCACATTTCAATTGCACATTAATGGATAtt                                                                           >  1:1222347/1‑150 (MQ=255)
                                                                              gTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCt                                                                         >  1:2881800/1‑150 (MQ=255)
                                                                               tCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCtt                                                                        <  1:51298/150‑1 (MQ=255)
                                                                               tCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCtt                                                                        <  1:2355243/150‑1 (MQ=255)
                                                                               tCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGAAATTCtt                                                                        <  2:1700230/150‑1 (MQ=255)
                                                                                 ccTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTta                                                                      <  1:1687004/150‑1 (MQ=255)
                                                                                 ccTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTta                                                                      <  2:506709/150‑1 (MQ=255)
                                                                                      ttCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTaataat                                                                 <  1:1143595/150‑1 (MQ=255)
                                                                                        cATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATATTTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct                                                               <  1:1040602/150‑1 (MQ=255)
                                                                                        cATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct                                                               <  1:2757848/150‑1 (MQ=255)
                                                                                           tACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTcgc                                                            >  1:956410/1‑150 (MQ=255)
                                                                                             cTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGa                                                          <  1:36888/150‑1 (MQ=255)
                                                                                             cTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGa                                                          >  2:2323053/1‑150 (MQ=255)
                                                                                               cACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACg                                                        <  2:1115212/150‑1 (MQ=255)
                                                                                                 cGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGtt                                                      >  1:358279/1‑150 (MQ=255)
                                                                                                    aaTCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCt                                                   >  1:4487085/1‑150 (MQ=255)
                                                                                                     aTCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCtt                                                  <  1:3027749/150‑1 (MQ=255)
                                                                                                     aTCGGTATTGTCATTATTCATTTATTAGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAACTCTCGACGTTTCta                                                  >  1:1064554/1‑149 (MQ=255)
                                                                                                       cGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTa                                                >  2:4123263/1‑150 (MQ=255)
                                                                                                       cGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTa                                                <  2:4221904/150‑1 (MQ=255)
                                                                                                        ggTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTAt                                               <  2:4382655/150‑1 (MQ=255)
                                                                                                        ggTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGATTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTAt                                               >  1:3250989/1‑150 (MQ=255)
                                                                                                         gTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATg                                              >  2:2260354/1‑150 (MQ=255)
                                                                                                          tATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGa                                             >  2:1659112/1‑150 (MQ=255)
                                                                                                          tATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGa                                             >  1:3543301/1‑150 (MQ=255)
                                                                                                           aTTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGat                                            >  2:3496885/1‑150 (MQ=255)
                                                                                                            ttGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGata                                           <  1:1279717/150‑1 (MQ=255)
                                                                                                              gTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaa                                         <  1:932370/150‑1 (MQ=255)
                                                                                                              gTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaa                                         <  2:3405229/150‑1 (MQ=255)
                                                                                                              gTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaa                                         >  1:3072328/1‑150 (MQ=255)
                                                                                                               tCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaat                                        <  2:4014931/150‑1 (MQ=255)
                                                                                                                catttttcatttttTTGTCCTTTGGAGGAGATTTAATCTCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaata                                       <  2:2657649/150‑1 (MQ=255)
                                                                                                                cATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaata                                       <  2:590945/150‑1 (MQ=255)
                                                                                                                 attattCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGataaataa                                      <  2:2833544/150‑1 (MQ=255)
                                                                                                                  ttattCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGAc                                                         <  1:1132582/130‑1 (MQ=255)
                                                                                                                      tCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATc                                 >  1:3464691/1‑150 (MQ=255)
                                                                                                                       cATTTTTTTGTCCATTGGAGGAGATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCa                                <  2:1884457/150‑1 (MQ=255)
                                                                                                                             tttGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGa                          >  2:2274127/1‑150 (MQ=255)
                                                                                                                                gTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATaa                       <  2:4422270/150‑1 (MQ=255)
                                                                                                                                gTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATaa                       >  1:4543134/1‑150 (MQ=255)
                                                                                                                                  ccATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATaaaa                     <  2:4380575/150‑1 (MQ=255)
                                                                                                                                     ttGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCa                  >  1:1600610/1‑150 (MQ=255)
                                                                                                                                      tGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCaa                 <  1:1229194/150‑1 (MQ=255)
                                                                                                                                      tGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCaa                 <  1:3897559/150‑1 (MQ=255)
                                                                                                                                        taggagATTTAATTCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTCCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCaaaa               <  2:3881364/149‑1 (MQ=255)
                                                                                                                                        gaggagATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCaaaa               <  2:3570813/150‑1 (MQ=255)
                                                                                                                                         aggagATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAAt              >  1:4254900/1‑150 (MQ=255)
                                                                                                                                          ggagATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATg             <  2:764007/150‑1 (MQ=255)
                                                                                                                                          ggagATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATg             <  1:2199385/150‑1 (MQ=255)
                                                                                                                                           gagaTTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATga            >  1:2161287/1‑150 (MQ=255)
                                                                                                                                           gagaTTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATga            >  1:4071960/1‑150 (MQ=255)
                                                                                                                                             gaTTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATgaga          <  2:1615071/150‑1 (MQ=255)
                                                                                                                                               tttAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGaaa        <  2:64648/150‑1 (MQ=255)
                                                                                                                                                  aaTCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAat     >  1:2083221/1‑150 (MQ=255)
                                                                                                                                                  aaTCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAat     >  1:26196/1‑150 (MQ=255)
                                                                                                                                                  aaTCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAat     >  1:604649/1‑150 (MQ=255)
                                                                                                                                                   aTCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAata    <  2:2185778/150‑1 (MQ=255)
                                                                                                                                                    tCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCTACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAatat   <  2:3663187/150‑1 (MQ=255)
                                                                                                                                                    tCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAatat   >  2:2090365/1‑150 (MQ=255)
                                                                                                                                                     cccTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAATATg  <  1:2506736/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGCGACGTTTCTTTATGATAAATAATAATCAAATTGATAAAATCAAAATGAGAAAAATATG  >  CP009273/3837028‑3837326

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: ‑