Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,875,281 C→A 6.4% intergenic (‑121/‑157) dgoR ← / → yidX D‑galactonate catabolism operon transcriptional repressor/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,2810CA 6.4% 437.9 / 12.5 145intergenic (‑121/‑157)dgoR/yidXD‑galactonate catabolism operon transcriptional repressor/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (76/59);  new base A (2/8);  total (78/67)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 148 total aligned reads displayed.

CAATGCGATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGG  >  NC_000913/3875133‑3875430
                                                                                                                                                    |                                                                                                                                                     
cAATGCGATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACt                                                                                                                                                      >  1:1615143/1‑150 (MQ=255)
   tGCGATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTtat                                                                                                                                                   >  1:1682913/1‑150 (MQ=255)
     cGATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTatat                                                                                                                                                 >  2:953024/1‑150 (MQ=255)
     cGATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTatat                                                                                                                                                 >  1:671791/1‑150 (MQ=255)
      gATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATAtt                                                                                                                                                <  1:1408479/150‑1 (MQ=255)
           gTTTTATTGTGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTTCCCCGCTCGATGCTTAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGGTTGCCGTGCGTGTAATGTAGTACAAACTTATATtgttgt                                                                                                                                           <  1:234305/150‑1 (MQ=255)
           gTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACTCTGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATtgttgt                                                                                                                                           <  2:921374/150‑1 (MQ=255)
           gTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATtgttgt                                                                                                                                           >  2:1623506/1‑150 (MQ=255)
             tttATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACt                                                                                                                                                      >  2:100507/1‑137 (MQ=255)
             tttATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACt                                                                                                                                                      <  1:100507/137‑1 (MQ=255)
              ttATTGAGAGTCATGGAGAGTACTTGTGCTCTTTATGTTATTCTCTACTTTACCACGATCTATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGtact                                                                                                                                        <  2:458768/150‑1 (MQ=255)
               tATTGAGATTCATGGAGAGTCCTTGTGCTCTGTATTTTCTTCTCTACTTTACCCAGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGtacta                                                                                                                                       <  2:853740/150‑1 (MQ=255)
               tATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGtacta                                                                                                                                       <  2:726388/150‑1 (MQ=255)
                aTTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAACAATATAAGTTTGtactac                                                                                                                                      <  2:1633504/150‑1 (MQ=255)
                aTTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGtactac                                                                                                                                      >  1:760411/1‑150 (MQ=255)
                 ttGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTTCTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACa                                                                                                                                     <  1:841260/150‑1 (MQ=255)
                 ttGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTACTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAACAATATAAGTTTGTACTACa                                                                                                                                     <  2:1564583/150‑1 (MQ=255)
                  tGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACaa                                                                                                                                    >  2:1074392/1‑150 (MQ=255)
                  tGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACaa                                                                                                                                    >  1:1365549/1‑150 (MQ=255)
                   gagaGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGATATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAAt                                                                                                                                   >  1:870013/1‑150 (MQ=255)
                    agcgTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTAACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACAACCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTAGAAACTTATATTTTTGTACTACAAtt                                                                                                                                  >  1:284774/4‑150 (MQ=255)
                    agagTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAAtt                                                                                                                                  >  2:124684/1‑150 (MQ=255)
                    agagTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAAtt                                                                                                                                  >  1:288230/1‑150 (MQ=255)
                     gagTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAACAATATAAGTTTGTACTACattac                                                                                                                                 <  2:41795/150‑5 (MQ=255)
                     gagTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAAttt                                                                                                                                 >  2:503903/1‑150 (MQ=255)
                       gTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAACAATATAAGTTTGTACTACattacac                                                                                                                               <  2:1268430/150‑7 (MQ=38)
                       gTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAACAATATAAGTTTGTACTACattacac                                                                                                                               <  1:1441686/150‑7 (MQ=38)
                       gTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAg                                                                                                                               <  2:532037/150‑1 (MQ=255)
                        tCATGGAGAGTCCTTGTGCTCTGTATGTTATTATATAATTTTATACGATCTATTCATATAGCTTCAACTTTTTAGTATCATTGTGATGACGTTGGTGTTTTGCGTGTGTTTATTGTAGTACAAACTTATATTGTTGTACTACAATTTAGa                                                                                                                              <  2:39993/150‑1 (MQ=255)
                         cATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGAt                                                                                                                             >  1:404089/1‑150 (MQ=255)
                           tGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATca                                                                                                                           >  2:150625/1‑150 (MQ=255)
                            ggAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATcac                                                                                                                          <  2:1466583/150‑1 (MQ=255)
                                  gccTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAg                                                                                                                    <  2:490256/149‑1 (MQ=255)
                                   ccTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGa                                                                                                                   <  1:150625/150‑1 (MQ=255)
                                         gCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATg                                                                                                             <  1:1515831/150‑1 (MQ=255)
                                          ctctGTATGTTCTTCTCTACTTTACCCCGTTCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTACGCGTTTGCCGGGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCCCAAAATGAACAATGc                                                                                                            >  2:637856/1‑150 (MQ=255)
                                              gTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATaa                                                                                                        <  2:749450/150‑1 (MQ=255)
                                               tATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATaaa                                                                                                       >  2:247251/1‑150 (MQ=255)
                                               tATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATaaa                                                                                                       >  2:1607926/1‑150 (MQ=255)
                                                     cttctCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATTAAAAAATg                                                                                                 <  1:960682/150‑1 (MQ=255)
                                                       tctcTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGAc                                                                                               >  2:664576/1‑150 (MQ=255)
                                                             cTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTACGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGt                                                                                         >  1:1186004/1‑150 (MQ=255)
                                                                 aCCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCggg                                                                                     <  1:1536625/150‑1 (MQ=255)
                                                                  ccccGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGc                                                                                    >  1:814283/1‑150 (MQ=255)
                                                                  ccccGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGc                                                                                    >  2:1289170/1‑150 (MQ=255)
                                                                  ccccGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGc                                                                                    >  1:1655075/1‑150 (MQ=255)
                                                                   cccGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCa                                                                                   >  2:681673/1‑150 (MQ=255)
                                                                   cccGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCa                                                                                   <  1:1074392/150‑1 (MQ=255)
                                                                   accGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCa                                                                                   <  2:535898/149‑1 (MQ=255)
                                                                    ccGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAg                                                                                  >  1:1904/1‑150 (MQ=255)
                                                                     cgatcgatGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGa                                                                                 <  2:595676/150‑1 (MQ=255)
                                                                     cgatcgatGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGa                                                                                 <  1:124684/150‑1 (MQ=255)
                                                                       ttcgatGCATAACGCGGCAACTTTGTAGTACCTGCGTGATGACGTTCGCGTTTGCCGTGCGTGTATTGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGaaa                                                                               <  2:263995/149‑1 (MQ=255)
                                                                          gatGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCt                                                                            >  1:1699156/1‑150 (MQ=255)
                                                                           atGCATAACGCTGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTg                                                                           <  1:608936/150‑1 (MQ=255)
                                                                           atGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTg                                                                           >  2:467895/1‑150 (MQ=255)
                                                                            tGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGa                                                                          >  1:204364/1‑150 (MQ=255)
                                                                            tGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACAAGCGTCGGGCAGAAATCTGa                                                                          >  2:1043597/1‑150 (MQ=255)
                                                                             gCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGaa                                                                         >  1:99479/1‑150 (MQ=255)
                                                                                tAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAg                                                                      >  1:94129/1‑150 (MQ=255)
                                                                                    gcgGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGat                                                                  >  1:772362/1‑150 (MQ=255)
                                                                                    gcgGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGat                                                                  >  2:1411608/1‑150 (MQ=255)
                                                                                     cgGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGata                                                                 <  2:979821/150‑1 (MQ=255)
                                                                                      ggCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGatat                                                                >  2:1636042/1‑150 (MQ=255)
                                                                                         aaCTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAg                                                             >  2:1362341/1‑150 (MQ=255)
                                                                                          aCTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAgg                                                            <  2:259999/150‑1 (MQ=255)
                                                                                           cTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGgc                                                           >  1:1610345/1‑150 (MQ=255)
                                                                                             ttGTAGTAACAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGgcgc                                                         <  2:1302565/150‑1 (MQ=255)
                                                                                               gtagtaCCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTa                                                       >  2:1246732/1‑150 (MQ=255)
                                                                                                 agtaCCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTaaa                                                     >  1:1009127/1‑150 (MQ=255)
                                                                                                  gtaCCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAAc                                                    >  2:1485508/1‑150 (MQ=255)
                                                                                                    aCCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACag                                                  >  2:1491450/1‑150 (MQ=255)
                                                                                                    aCCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACag                                                  <  2:1614814/150‑1 (MQ=255)
                                                                                                         cGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAggg                                             >  2:1540944/1‑150 (MQ=255)
                                                                                                         cGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAggg                                             <  2:1699156/150‑1 (MQ=255)
                                                                                                                 cGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAagaaga                                     <  2:814283/150‑1 (MQ=255)
                                                                                                                  gTTCGCGTTTGCCTTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAagag                                    <  1:29283/150‑1 (MQ=255)
                                                                                                                  gTTAGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAagag                                    <  1:527081/150‑1 (MQ=255)
                                                                                                                    tCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGTGCTAAACAGGAGGGAAAGAAGAGTa                                  <  2:1669346/150‑1 (MQ=255)
                                                                                                                           tGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGgcg                                                          >  1:1662509/1‑119 (MQ=255)
                                                                                                                            gCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCa                          <  2:1112963/150‑1 (MQ=255)
                                                                                                                              cgtgcgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAAc                        >  2:202526/1‑150 (MQ=255)
                                                                                                                               gtgcgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACg                       <  2:204364/150‑1 (MQ=255)
                                                                                                                                tgcgtgTAATGTAGTACAAACTTATATTTTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACgg                      <  1:1179255/150‑1 (MQ=255)
                                                                                                                                  tctaaattGTAGTACAACAATATAAGTTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCt                    >  1:537609/8‑150 (MQ=38)
                                                                                                                                   gtgtAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCtt                   <  2:1619371/150‑1 (MQ=255)
                                                                                                                                     gtAATGTAGTACAAAATTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAg                 <  1:1623506/150‑1 (MQ=255)
                                                                                                                                      tAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCttt                            >  2:985007/1‑138 (MQ=255)
                                                                                                                                      tAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCttt                            <  1:985007/138‑1 (MQ=255)
                                                                                                                                      tAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGc                <  1:681673/150‑1 (MQ=255)
                                                                                                                                      tAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGc                >  2:1077021/1‑150 (MQ=255)
                                                                                                                                      tAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGc                <  2:13910/150‑1 (MQ=255)
                                                                                                                                       aaTGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCt               <  1:664576/150‑1 (MQ=255)
                                                                                                                                        aTGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGATAGAAGAGTATGCTTTCAACGGCTTAggtt              >  2:1152602/1‑147 (MQ=255)
                                                                                                                                         tGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTAc             <  1:223708/150‑1 (MQ=255)
                                                                                                                                           tagtaCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCAAAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGAATAACAGGCGATAAACAGGAGGTAAAGAAGAGTATGCTTTCAACGGCTATCCTACTc           >  1:730131/1‑150 (MQ=255)
                                                                                                                                                cAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGt         >  1:367970/1‑147 (MQ=255)
                                                                                                                                                cAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTa      <  1:1289170/150‑1 (MQ=255)
                                                                                                                                                  aaCTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTaaa    >  1:1419938/1‑150 (MQ=255)
                                                                                                                                                    cTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAgg  <  2:1170315/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATGCGATCGGTTTTATTGAGAGTCATGGAGAGTCCTTGTGCTCTGTATGTTCTTCTCTACTTTACCCCGATCGATGCATAACGCGGCAACTTTGTAGTACCAGCGTGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGG  >  NC_000913/3875133‑3875430

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