Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,330 G→C 6.5% intergenic (+41/‑81) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3300GC 6.5% 342.6 / 11.0 117intergenic (+41/‑81)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base G (56/53);  new base C (4/4);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 118 total aligned reads displayed.

GGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGG  >  NC_000913/4612183‑4612478
                                                                                                                                                   |                                                                                                                                                    
ggTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCg                                                                                                                                                    >  1:22780/1‑150 (MQ=255)
 gTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCgg                                                                                                                                                   >  2:1267071/1‑150 (MQ=255)
    tggtggTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                                                                                >  2:1154054/1‑150 (MQ=255)
      gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGCCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTgg                                                                                                                                              >  1:350222/1‑150 (MQ=255)
         gtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTTTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCtt                                                                                                                                           <  2:1755722/150‑1 (MQ=255)
           cTGGCCGGTACCGCTGCCGGCGCAGCTAAAATTGGCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTCATCCAAAGCCAGTCCAGCGGACTGGCtttt                                                                                                                                         <  1:721832/150‑1 (MQ=255)
           cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCtttg                                                                                                                                         <  1:1279789/150‑2 (MQ=255)
            tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTTTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCtttgg                                                                                                                                        <  2:1050411/150‑3 (MQ=255)
            tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCtttgg                                                                                                                                        <  2:1192039/150‑3 (MQ=255)
             ggCCGGTACCGCTGCAGGCGCAGCTTAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTGGc                                                                                                                                       >  1:239407/1‑150 (MQ=255)
              gCCGGTACCGCTGCCGGCGCTGCTAAACTTGTCTTTGTCTTCGGGAGTTTTCTTATCCTGGTGAGTTTGTTCATGGGCGGAACACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCg                                                                                                                                      <  1:1710422/150‑1 (MQ=255)
              gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGGCGTCGGGATTATTGTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGGCCAGCGGACTGGCTTTTGCg                                                                                                                                      <  2:1662412/150‑1 (MQ=255)
               ccGGTACCGCTGCAGGCTAAGCTAAAATTGTCTTTGTCGTCGTGATTATTCTGTTCCTGTTGAGTTTGTTCATGGGCCTAAAACGACCATAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                     <  2:1724107/150‑1 (MQ=255)
                  gTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCATTCCAGCGGACTGGCTTTTGCGGttt                                                                                                                                  <  2:302589/150‑1 (MQ=255)
                  gTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGttt                                                                                                                                  <  1:640874/150‑1 (MQ=255)
                  gTACCGCTGCAGGCGCAGCTAAAATTGACTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGCCCCTAGATTTCAGACCATTAAGATAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGtgt                                                                                                                                  >  1:1610079/1‑148 (MQ=255)
                   tACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGtttt                                                                                                                                 <  1:450124/150‑1 (MQ=255)
                        cTGCAGGCTCAGGTAAAATTTTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTATATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCTGACTTGCTTTTGCGGTTTTAGCGa                                                                                                                            <  2:541588/150‑1 (MQ=255)
                         tGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaa                                                                                                                           <  2:1606736/150‑1 (MQ=255)
                          gCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaat                                                                                                                          >  2:20828/1‑150 (MQ=255)
                             ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                       >  2:368998/1‑150 (MQ=255)
                              gcgcAGCTAAAAGTGTCTTTGTCGTCGGGATTATTCTGGTCCTGGTGAGTTTGTTCATGGGCCGAACACGACCCTAGATTTCATACCATTAAGAAAATTTTTTCCAATGCCAGTCCAGCGGACTGGCTTTTTCGGTTTTAGCGAataata                                                                                                                      <  1:1254720/150‑1 (MQ=255)
                                    cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCGTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                                >  2:622829/1‑150 (MQ=255)
                                      aaaaTTGTCTTTGTCGTCGGGATTATTCTGATACTCGTGAGTTTGTTCATGGGCCGAAACCCACCCTAGATTTACGGCAATTAAGAAAATTTTAACCAAAGCCAGTCCAGCGGACTGGATTTTGCGGATTTAGCGAATAATAAATTGCGt                                                                                                              >  1:1306615/1‑150 (MQ=255)
                                      aaaaATGTCTATGTCCTCGGGAATATTCTGTTCCAGGTCAGATTCTTCAAGTGCCGAAAACCACCCTAGATTTCAGACCATTAAGAAAATATAATACAAAGCCACACCACCGCACTGGCTTTTCCGGTATTAGCGAATAATAAAttccct                                                                                                              >  1:1306608/1‑146 (MQ=255)
                                       aaaTTGTCTTTGTCGTCGGGATTATTCTGGTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGtt                                                                                                             <  2:1095856/150‑1 (MQ=255)
                                           tGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATCGCGTTACtt                                                                                                         <  1:849536/150‑1 (MQ=255)
                                             tCTTTGTGGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTg                                                                                                       <  2:1781212/150‑1 (MQ=255)
                                             tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTg                                                                                                       <  1:1484093/150‑1 (MQ=255)
                                              cTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                      >  1:372110/1‑150 (MQ=255)
                                               tttGTCGTCGGGATTATTCTGTTTCTGGTGATTTTGTTTAAGGGCCGTAAAAGAACCTATTTTTTATTCCATTAAGTAAATTTTATCCAAAGCCATTCCCGGGGTCTGGCTTTTGGGGTTTTATGGAATAATAAATTGCGTTACTTTGTc                                                                                                     <  2:1039296/150‑1 (MQ=255)
                                                ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCa                                                                                                    <  1:474249/150‑1 (MQ=255)
                                                     gtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCtt                                                                                               <  1:844304/150‑1 (MQ=255)
                                                     gtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATAAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCtt                                                                                               >  1:847873/1‑150 (MQ=255)
                                                     gtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGAACATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCtt                                                                                               <  1:80404/150‑1 (MQ=255)
                                                       cgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                             >  1:775384/1‑150 (MQ=255)
                                                         ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCTTTACTTTGTCATTCTTTcaa                                                                                           >  2:1544212/1‑149 (MQ=255)
                                                           attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                         <  1:1337978/150‑1 (MQ=255)
                                                             tattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaaca                                                                                       >  2:315984/1‑150 (MQ=255)
                                                               ttCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaa                                                                                     >  1:1003494/1‑150 (MQ=255)
                                                                   gTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                 <  2:385766/150‑1 (MQ=255)
                                                                   gTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                 >  2:548794/1‑150 (MQ=255)
                                                                     tcatGGGGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCATAAGATTAATTAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTAACTTTGTCATTCTTTCACAACAAAAAGAgg                                                                               <  2:51191/147‑1 (MQ=255)
                                                                        tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                            <  2:467807/150‑1 (MQ=255)
                                                                            gAGTTTGTTCATGGGCCGCAAACGCCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                                                        >  1:47723/1‑150 (MQ=255)
                                                                            gAGTTTGTTCATGGGACGAAAACGACCCAAGATTTCAGACCATTATGAAAATTTTATCCAACGCCAGTCCAGCGGACTGGCTGTTGCGGGTGTAGCGAATAATAACTTGCGTTACTTTGTCATTCTTTCCCAACAAACACAGGACGGCag                                                                        >  2:1232822/1‑150 (MQ=255)
                                                                              gTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                      >  1:632293/1‑150 (MQ=255)
                                                                                 tGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACACCAAAAACAGGAAGGCAGAGGTg                                                                   >  2:918041/1‑150 (MQ=255)
                                                                                 tGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                   <  1:1267071/150‑1 (MQ=255)
                                                                                   ttCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTggg                                                                 <  2:714360/150‑1 (MQ=255)
                                                                                     cctgcgCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCCGCGGACGGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTCTCATTCAATCACAACAAAAACAGGAAGCCAGAGGTGGGGc                                                               <  2:294307/145‑1 (MQ=255)
                                                                                        gggCTGTAAAAGACCCTATATTTGAGACCATTAAGAAAATTTTATGCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGATAAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                            <  1:336958/150‑1 (MQ=255)
                                                                                        gggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                            >  2:1664135/1‑150 (MQ=255)
                                                                                         ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                           >  2:851126/1‑150 (MQ=255)
                                                                                         ggCCGAAAACGACCCGAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCCGCGGTCTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACATCAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                           >  1:1146750/1‑150 (MQ=255)
                                                                                          gCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACGTTGTCATGCGTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGa                                                          >  1:35710/1‑150 (MQ=255)
                                                                                           ccGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGaa                                                         >  1:1784846/1‑150 (MQ=255)
                                                                                             gAAAACGACCCTAGTTTTCAGACCATTAAGAAAATTTTATCCAAAGCCTGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                       <  2:1630661/150‑1 (MQ=255)
                                                                                             gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                       <  1:1237283/150‑1 (MQ=255)
                                                                                              aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                      <  2:604971/150‑1 (MQ=255)
                                                                                               aaaCGACCCTAGATTTCAGACCATTAAGAAACTTTTATCCAAAGCCAGGCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAATCAGGAAGGCAGAGGTGGGGCAGCTAATAcc                                                     <  1:185518/150‑1 (MQ=255)
                                                                                               aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCCGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAcc                                                     <  1:766786/150‑1 (MQ=255)
                                                                                                aaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                    <  1:20828/150‑1 (MQ=255)
                                                                                                 aCGAACCTAGATTTCTGACCTTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTg                                                   <  2:912621/150‑1 (MQ=255)
                                                                                                    aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTaa                                                >  1:1403261/1‑150 (MQ=255)
                                                                                                        tAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                            >  1:1084942/1‑150 (MQ=255)
                                                                                                            tttCAGAGCATTAAGAATTTTTTATCCAAAGCGTGTCCAGCGGACTGGCTTTTGCGGTTTTAGGGAATATTTAATTGCGTTAATTTGTCATTCTTTCTCAACATAAACAGGAAGGCAGAGGTGGGGGAGCGAATTCCTGTAACGCTTGGt                                        <  2:1306608/150‑1 (MQ=255)
                                                                                                            tttCAGACCATTAAGAAAATTTTATCCAAAGCGAGTCCAGCGGACTTGCTTTTGCGGTTTTATCGAATAATAAATTGCGTTACTTTGTCATTCTATCACAACAAAAACAGGAAGGCAGAGGTGTGGCAGCGAATACCTGTAACGCTTGGt                                        <  2:1306615/150‑1 (MQ=255)
                                                                                                            tttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGt                                        >  2:409544/1‑150 (MQ=255)
                                                                                                              tCAGAACTTTAAGAAAATTTTATCCAAAGCCAGTCCAGCTGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTTCTTTGTCATTCTTTCACAACAAAAACAGGATGGGAGAGGTGGGGCAGCGAATACCTGGAACGCTTGGTaa                                      <  2:900542/150‑1 (MQ=255)
                                                                                                               cAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAat                                     >  2:698617/1‑150 (MQ=255)
                                                                                                                aGACGATTAAGAAAATTTTATCCAAATCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAata                                    <  1:276135/150‑1 (MQ=255)
                                                                                                                   acATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                                 <  2:817174/149‑1 (MQ=255)
                                                                                                                    cATTAAGAAAATTTTATCCAAAGCCAGTCCAGCTGACTGGTTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACTACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                <  1:316326/150‑1 (MQ=255)
                                                                                                                    cATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGATAGGTAATATTgc                                >  1:1133978/1‑150 (MQ=255)
                                                                                                                     aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                               >  1:493487/1‑150 (MQ=255)
                                                                                                                       tAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGcc                             >  2:504477/1‑150 (MQ=255)
                                                                                                                       tAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGAGGTTTTAGAGATTAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGcc                             >  2:168656/1‑150 (MQ=255)
                                                                                                                         aGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTACCGCTTGGTAATATTGCGCCGt                           >  2:1203923/1‑150 (MQ=255)
                                                                                                                         aGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt                           <  1:1773844/150‑1 (MQ=255)
                                                                                                                          gAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                          >  1:997101/1‑150 (MQ=255)
                                                                                                                              aTTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGTCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCg                      <  2:1289459/150‑1 (MQ=255)
                                                                                                                               ttttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCCGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGc                     <  1:94840/150‑1 (MQ=255)
                                                                                                                                   cgcaaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg                 >  1:1409346/5‑150 (MQ=255)
                                                                                                                                      cAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc              >  2:1816527/1‑150 (MQ=255)
                                                                                                                                       aaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg             >  2:1569123/1‑150 (MQ=255)
                                                                                                                                         aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt           >  1:565646/1‑150 (MQ=255)
                                                                                                                                          gacAGTTCAGAGTACTGGCTTTTTCTGTTTTATCTAATAATAAATTTCTTTACTTTGTTATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACATGTATCGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtta          <  1:1198256/148‑2 (MQ=255)
                                                                                                                                          gCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTc          >  2:1790964/1‑150 (MQ=255)
                                                                                                                                          gCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTc          <  1:1203923/150‑1 (MQ=255)
                                                                                                                                           gcAGTCCTGCTGACTGGCTTTTGGGGTTTTAGCGAATAATAAATTTCGTTACTTTGTAATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTcc         <  2:606131/149‑1 (MQ=255)
                                                                                                                                           ccAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCTCAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTcc         >  2:705088/1‑150 (MQ=255)
                                                                                                                                              gTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGc      <  2:1778544/150‑1 (MQ=255)
                                                                                                                                               tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGcc     >  1:426312/1‑150 (MQ=255)
                                                                                                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCt    >  1:1122767/1‑150 (MQ=255)
                                                                                                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCt    >  1:761247/1‑150 (MQ=255)
                                                                                                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCt    >  2:1561637/1‑150 (MQ=255)
                                                                                                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGGGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGAGCCGTTGTCGCAAAGGCCGTTCCAGCCt    >  2:1695107/1‑150 (MQ=255)
                                                                                                                                                 cAGCGGACTGGCTTTTGCGGCTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCTTTCCAGCCTg   <  1:1283143/150‑1 (MQ=255)
                                                                                                                                                 cAGCGGACTGGCTTTTGAGGATTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg   <  1:1702344/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGG  >  NC_000913/4612183‑4612478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: