Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 580,102 G→T 13.7% intergenic (+216/+34) appY → / ← ompT global transcriptional activator; DLP12 prophage/DLP12 prophage; outer membrane protease VII (outer membrane protein 3b)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,1020GT13.7% 760.7 / 57.5 223intergenic (+216/+34)appY/ompTglobal transcriptional activator; DLP12 prophage/DLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (93/98);  new base T (12/20);  total (105/118)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

CTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTATTACC  >  CP009273/579955‑580249
                                                                                                                                                   |                                                                                                                                                   
cTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGAt                                                                                                                                                   >  2:4623960/1‑150 (MQ=255)
cTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGAt                                                                                                                                                   <  1:2174423/150‑1 (MQ=255)
 tGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTAtt                                                                                                                                                  <  1:1813879/150‑1 (MQ=255)
      ccATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCAc                                                                                                                                             <  2:3038424/150‑1 (MQ=255)
       cATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCAcc                                                                                                                                            <  1:1197424/150‑1 (MQ=255)
        aTTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCt                                                                                                                                           >  2:2237260/1‑150 (MQ=255)
        aTTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCt                                                                                                                                           >  1:2426114/1‑150 (MQ=255)
         ttAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTc                                                                                                                                          <  2:4049116/150‑1 (MQ=255)
         ttAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTc                                                                                                                                          >  1:2013234/1‑150 (MQ=255)
          tAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCg                                                                                                                                         >  1:2860906/1‑150 (MQ=255)
          tAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCg                                                                                                                                         <  1:2506082/150‑1 (MQ=255)
            aTTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCggg                                                                                                                                       >  1:169759/1‑150 (MQ=255)
               tGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGaa                                                                                                                                    <  2:151838/150‑1 (MQ=255)
                gTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGaaa                                                                                                                                   >  1:722147/1‑150 (MQ=255)
                gTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGaaa                                                                                                                                   <  1:1847719/150‑1 (MQ=255)
                          tGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTgg                                                                                                                         >  1:975670/1‑150 (MQ=255)
                           gCTGTCGCCTATTTTATACATACTATCATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGc                                                                                                                        <  2:3235098/150‑1 (MQ=255)
                              gTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGtt                                                                                                                     >  2:1717404/1‑150 (MQ=255)
                                cGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCt                                                                                                                   >  2:1676766/1‑150 (MQ=255)
                                  ccTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTa                                                                                                                 <  1:3843645/150‑1 (MQ=255)
                                    tATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCTGGGAAATTTTAGTTGGCGTTCTTaaa                                                                                                               >  2:565003/1‑150 (MQ=255)
                                        ttATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAAtgt                                                                                                           >  1:2495014/1‑150 (MQ=255)
                                            acatacTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTAc                                                                                                       >  2:2491397/1‑150 (MQ=255)
                                                acTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTaa                                                                                                   >  2:4274660/1‑150 (MQ=255)
                                                        tGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACcagca                                                                                           <  1:4352606/150‑1 (MQ=255)
                                                        tGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACcagca                                                                                           >  1:254168/1‑150 (MQ=255)
                                                         gATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACcagcag                                                                                          <  1:2768593/150‑1 (MQ=255)
                                                          aTGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCagt                                                                                         >  1:1219190/1‑150 (MQ=255)
                                                          aTGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCagt                                                                                         <  1:88699/150‑1 (MQ=255)
                                                           tGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCagta                                                                                        <  2:1578357/150‑1 (MQ=255)
                                                            ggTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCagtag                                                                                       >  2:576112/1‑150 (MQ=255)
                                                            ggTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCagtag                                                                                       <  1:1244347/150‑1 (MQ=255)
                                                              ttttCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGtg                                                                                     <  2:2931102/150‑1 (MQ=255)
                                                              ttttCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGtg                                                                                     <  1:3464788/150‑1 (MQ=255)
                                                                 tCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGtgatg                                                                                  <  1:3951379/150‑1 (MQ=255)
                                                                   tATGTGATTTATTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCGAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGaa                                                                                <  2:4497686/150‑1 (MQ=255)
                                                                   tATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGaa                                                                                <  2:3043629/150‑1 (MQ=255)
                                                                   tATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGaa                                                                                >  2:1161687/1‑150 (MQ=255)
                                                                     tgtgATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGt                                                                              >  1:1405357/1‑150 (MQ=255)
                                                                      gtgATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGtt                                                                             >  2:4301706/1‑150 (MQ=255)
                                                                         aTTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTTTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTtata                                                                          <  2:50628/150‑1 (MQ=255)
                                                                         aTTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTtata                                                                          <  2:3016634/150‑1 (MQ=255)
                                                                             aGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGttt                                                                      >  2:867077/1‑150 (MQ=255)
                                                                             aGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGttt                                                                      <  1:3337233/150‑1 (MQ=255)
                                                                              gTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGtttt                                                                     <  2:3461475/150‑1 (MQ=255)
                                                                               ttAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTc                                                                    <  1:2942239/150‑1 (MQ=255)
                                                                                taataaCCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCt                                                                   >  2:2835351/1‑150 (MQ=255)
                                                                                 aataaCCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCta                                                                  >  1:3068178/1‑150 (MQ=255)
                                                                                   taaCCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCtata                                                                >  1:3770162/1‑150 (MQ=255)
                                                                                     aCCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATAcc                                                              <  2:3247681/150‑1 (MQ=255)
                                                                                          cTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTc                                                         >  1:1733373/1‑150 (MQ=255)
                                                                                            gggTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCa                                                       >  1:500829/1‑150 (MQ=255)
                                                                                            gggTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCa                                                       <  2:2831176/150‑1 (MQ=255)
                                                                                             ggTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCAt                                                      <  1:4534573/150‑1 (MQ=255)
                                                                                              gTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCAtt                                                     >  2:3507733/1‑150 (MQ=255)
                                                                                               tttATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCAttt                                                    <  2:1082721/150‑1 (MQ=255)
                                                                                               tttATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCAttt                                                    >  1:708633/1‑150 (MQ=255)
                                                                                                ttATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCAtttt                                                   >  1:476183/1‑150 (MQ=255)
                                                                                                ttATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCAtttt                                                   >  1:462173/1‑150 (MQ=255)
                                                                                                  aTTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTg                                                 >  2:4152578/1‑150 (MQ=255)
                                                                                                  aTTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTg                                                 <  2:685068/150‑1 (MQ=255)
                                                                                                    tttAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCt                                               <  2:2013248/150‑1 (MQ=255)
                                                                                                       aaGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTa                                            <  1:707868/150‑1 (MQ=255)
                                                                                                       aaGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTa                                            <  2:3934685/150‑1 (MQ=255)
                                                                                                        aGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAg                                           <  1:4314054/150‑1 (MQ=255)
                                                                                                          ggTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTc                                         <  1:805550/150‑1 (MQ=255)
                                                                                                           gTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCt                                        >  2:2663300/1‑150 (MQ=255)
                                                                                                            ttAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTg                                       >  1:3670578/1‑150 (MQ=255)
                                                                                                                 tGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAgtg                                  >  2:886288/1‑150 (MQ=255)
                                                                                                                 tGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAgtg                                  >  2:238292/1‑150 (MQ=255)
                                                                                                                   ttACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGtt                                <  1:3623809/150‑1 (MQ=255)
                                                                                                                    tACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGtta                               >  1:320289/1‑150 (MQ=255)
                                                                                                                      cATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGttatt                             <  1:841962/150‑1 (MQ=255)
                                                                                                                       aTTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGttatta                            <  2:2053360/150‑1 (MQ=255)
                                                                                                                        ttGAAATGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTtattat                           <  2:2082247/150‑1 (MQ=255)
                                                                                                                         tGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTattatt                          <  1:2721507/150‑1 (MQ=255)
                                                                                                                          gAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATtg                         <  1:645039/150‑1 (MQ=255)
                                                                                                                          gAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATtg                         >  1:29659/1‑150 (MQ=255)
                                                                                                                           aaaTGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATtgt                        <  2:3600099/150‑1 (MQ=255)
                                                                                                                            aaTGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATtgtg                       <  1:4129204/150‑1 (MQ=255)
                                                                                                                            aaTGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATtgtg                       <  1:2833847/150‑1 (MQ=255)
                                                                                                                             aTGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGa                      <  1:2146532/150‑1 (MQ=255)
                                                                                                                             aTGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGa                      >  1:1939959/1‑150 (MQ=255)
                                                                                                                              tGGCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGAt                     >  1:3296198/1‑150 (MQ=255)
                                                                                                                               ggCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATc                    <  2:4076559/150‑1 (MQ=255)
                                                                                                                               ggCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATc                    >  1:587831/1‑150 (MQ=255)
                                                                                                                                gCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCa                   <  1:4124182/150‑1 (MQ=255)
                                                                                                                                gCTAGTTATTCCCCGGGGCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCa                   >  2:1319489/1‑150 (MQ=255)
                                                                                                                                 cTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCAt                  >  1:1096386/1‑150 (MQ=255)
                                                                                                                                   aGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATaa                >  2:4142640/1‑150 (MQ=255)
                                                                                                                                   aGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATaa                >  1:4245918/1‑150 (MQ=255)
                                                                                                                                    gTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATaaa               >  1:4344846/1‑150 (MQ=255)
                                                                                                                                     ttATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAg              <  1:3712152/150‑1 (MQ=255)
                                                                                                                                     ttATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAg              <  2:708633/150‑1 (MQ=255)
                                                                                                                                          ccccGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAg         >  2:3053671/1‑150 (MQ=255)
                                                                                                                                              ggggCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTAtt     >  2:1296577/1‑150 (MQ=255)
                                                                                                                                              ggggCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTAtt     >  1:1821737/1‑150 (MQ=255)
                                                                                                                                              ggggCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTAtt     >  2:3733779/1‑150 (MQ=255)
                                                                                                                                                ggCTATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTATTAc   >  2:3851360/1‑150 (MQ=255)
                                                                                                                                                 gCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTATTAcc  <  1:2237260/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGAACCATTAATTTGTAATTGGGTTGCTGTCGCCTATTTTATACATACTATAATTGATGGTTTTCTATGTGATTTAGTTAATAACCTTCTGGGTTTATTTTAAGGGTTAATTGTTACATTGAAATGGCTAGTTATTCCCCGGGGCGATTTTCACCTCGGGGAAATTTTAGTTGGCGTTCTTAAAATGTGTACTTAAGACCAGCAGTAGTGATGAAGTTATAGTTTTCTATACCTGCTCCATTTTTGCTGTAGTCTGAAGTGTTATTATTGTGATCATAAAGTGAAGTATTACC  >  CP009273/579955‑580249

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: