Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,521,346 A→T 5.1% intergenic (+18/+247) lysV → / ← xapR tRNA‑Lys/transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,3460AT 5.1% 885.4 / 29.9 314intergenic (+18/+247)lysV/xapRtRNA‑Lys/transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base A (158/138);  new base T (6/10);  total (164/150)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 318 total aligned reads displayed.

ACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/2521198‑2521494
                                                                                                                                                    |                                                                                                                                                    
 ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                    <  2:1575603/150‑1 (MQ=255)
  cTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                                   <  2:4115549/150‑1 (MQ=255)
   tCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCCCCCAAAAgg                                                                                                                                                  >  2:3819139/1‑150 (MQ=255)
   tCCCTTACAAGGAGGGGGGCGGCGGTTCGATCCCGGCCGCACCCCCCACTGCGGGTCGTTAGCTCAGTTGGTAGCGCAGTTGCCGGGGAAGCCCTTGGTCGCAGGGGCGACTCCTGCACGACCCACCACTGTCCAAAAGCGCCCGAAAgg                                                                                                                                                  >  2:2736060/1‑150 (MQ=9)
    cccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                                                                                                                                 >  2:6278403/1‑150 (MQ=255)
     ccTTACTAGGTGGGGGTCGGGGGTTCGATCCCGTCATCACCCACCACTTGGGTTCGTTAGCTCAGTTGTTTTAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                                                                                                                                <  1:5535532/150‑1 (MQ=17)
      cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                               >  2:6341897/1‑150 (MQ=255)
         aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttt                                                                                                                                            <  2:179143/150‑1 (MQ=255)
         aCAAGGAGGGGGTCGGCGGTTCGAGCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGGTGCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGGGCttt                                                                                                                                            >  2:5899600/1‑150 (MQ=255)
           ccggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                          <  1:1328542/148‑1 (MQ=255)
           atggtggggggCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCTGGTTCGAATCGTGCACGACCGACCAATGTAAAAAAGCGCCCTAAAGGCGCttttt                                                                                                                                          <  1:6223098/145‑1 (MQ=18)
            aggtggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCTGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                         <  1:4308272/150‑1 (MQ=255)
            aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATACTGCACGACCCACCAAAGTAAAAAAGCGCCCAAAAGGCGCtttttt                                                                                                                                         >  1:1842653/1‑150 (MQ=255)
             ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGGGCTTTTTTa                                                                                                                                        <  1:3961202/150‑1 (MQ=255)
             ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                                                                                                        <  1:2754721/150‑1 (MQ=255)
                gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCTTTTTTAcat                                                                                                                                     <  2:3693996/150‑3 (MQ=255)
                  gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATc                                                                                                                                   >  2:642144/1‑150 (MQ=255)
                   ggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCt                                                                                                                                  <  2:4395313/150‑1 (MQ=255)
                    gTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTg                                                                                                                                 <  2:3586961/150‑1 (MQ=255)
                     tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTTTTACTATCTGc                                                                                                                                >  1:6028210/1‑150 (MQ=255)
                     tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTACTATCTGc                                                                                                                                >  1:1409129/1‑150 (MQ=255)
                     tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTACTATCTGc                                                                                                                                >  1:330404/1‑150 (MQ=255)
                         cggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGGTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAAAAAAGCGCCCTAATGGCGCTTTTTTAATATATGCGTTa                                                                                                                            >  1:3696358/1‑150 (MQ=17)
                          ggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                                                                                                           >  2:1719012/1‑150 (MQ=255)
                          ggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                                                                                                           <  2:5445930/150‑1 (MQ=255)
                           gTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                                                                                                          <  1:5179995/150‑1 (MQ=255)
                            ttCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCGAAAGGCGCTTTTTTACTATCTGCGATACTc                                                                                                                         >  1:1061942/1‑150 (MQ=255)
                              cGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                                                                                                       <  1:3273781/150‑1 (MQ=255)
                               gATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                                                                                                      <  1:642144/150‑1 (MQ=255)
                               gATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                                                                                                      <  1:4711229/150‑1 (MQ=255)
                                aTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAGACGCGAAg                                                                                                                     >  2:6676764/1‑150 (MQ=255)
                                 tCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                                                                                                                    >  1:3943086/1‑150 (MQ=255)
                                   ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                                  <  1:1425921/150‑1 (MQ=255)
                                   ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                                  >  1:6662212/1‑150 (MQ=255)
                                   ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                                  <  1:2626249/150‑1 (MQ=255)
                                   ccGTCATCACCCACCACGTCGGGTCGTGAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                                                                                                                  <  2:3710788/150‑1 (MQ=255)
                                    cGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                                                                                                                 <  1:2073404/150‑1 (MQ=255)
                                    cGTCATCACCCACCACTTCGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCCATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                                                                                                                 >  1:1174751/1‑150 (MQ=255)
                                        atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATATGCGATAGTCCAACATTCGAAc                                                                                                             >  1:1374531/1‑150 (MQ=255)
                                           acccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                          <  1:1907012/150‑1 (MQ=255)
                                            cccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgc                                                                                                         <  1:4860072/150‑1 (MQ=40)
                                             ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATTGTAGCAGGTTCGAATCCAGCACGACCCACCAAAGTAAAAAAGCGCCCGAAAGGCGCTTTTTTACTATCTGCGATACTCAAATATTAGAACCAgca                                                                                                        >  1:1974443/1‑150 (MQ=14)
                                             ccaccaCTTCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGTCTTTTAATCAATTGGTCGCAGGGTCGAATACTGCACCACCCACCAAAGTAAAAAAACGCCCAAAAGGCGCTTTTTTACTATCAGCGAAACAAAAAGATTCGAACCAgca                                                                                                        >  1:3052512/1‑150 (MQ=14)
                                              caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                       <  2:5952480/150‑1 (MQ=39)
                                              caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                       <  2:6338175/150‑1 (MQ=39)
                                                        ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGATTCCTGCACGACCCACCATTGTAATAAAGCGCCCTACAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                             <  1:5574052/150‑1 (MQ=255)
                                                          tCGTTAGCTCCGTTGGTAGAGCAGTTGACTTTTACTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                                                           <  1:6678883/150‑1 (MQ=17)
                                                                 cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                    <  1:3089536/150‑1 (MQ=32)
                                                                       tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                              <  1:1618143/150‑1 (MQ=17)
                                                                        ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                             >  1:2888130/1‑150 (MQ=32)
                                                                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCGAACATTCGACCCTCCCGCAGGTATGGCGTCCGCGCCGCGAaaca                                                                            >  2:1457294/1‑150 (MQ=11)
                                                                         gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGAAGGTTTGAGTTGAGCGCAGCGAaaca                                                                            >  2:5389369/1‑150 (MQ=17)
                                                                            gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGCAACAACg                                                                         >  1:4126479/1‑150 (MQ=17)
                                                                            gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTAGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                                                                         >  1:975554/1‑150 (MQ=17)
                                                                                    gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                                 >  1:745536/1‑150 (MQ=32)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              <  2:3276346/150‑1 (MQ=17)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  1:3715596/1‑150 (MQ=32)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  1:3490996/1‑150 (MQ=32)
                                                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                              >  2:3936251/1‑150 (MQ=32)
                                                                                         ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                            <  2:5284357/150‑1 (MQ=32)
                                                                                          tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTACAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                           <  2:1592669/150‑1 (MQ=17)
                                                                                            aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCATAAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa                                                         >  1:5923861/1‑150 (MQ=17)
                                                                                            aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGAGTTGAGCGCAGCGAAACAACGGATCAGATCGCGGCGa                                                         >  2:524562/1‑150 (MQ=255)
                                                                                              cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg                                                       >  1:4629991/1‑150 (MQ=32)
                                                                                               aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                      >  2:3393623/1‑150 (MQ=32)
                                                                                               aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                      >  2:6099559/1‑150 (MQ=32)
                                                                                               aaTTGGTCGCAGGTTCGAATCCTACACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                      >  1:2486645/1‑150 (MQ=17)
                                                                                                     tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAATAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAATCAAAGATTCGAACCTGCAGCAGCTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                                >  2:4547926/1‑150 (MQ=12)
                                                                                                          ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTGTGAGTTGAGCGCAGCGACCCAACGGAGCCGCTCGCGGCGACGGCCCGCAGGgcg                                           >  1:1793082/1‑150 (MQ=11)
                                                                                                           gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                          >  1:684363/1‑150 (MQ=32)
                                                                                                           gTTCGAAGCCTGCACCACCCACCAATGTAAAAACGCGCCCAAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGCGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                          <  2:1830082/150‑1 (MQ=11)
                                                                                                            ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                         >  2:150775/1‑150 (MQ=32)
                                                                                                            tccGAATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCTAATATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                         <  2:4071394/148‑1 (MQ=255)
                                                                                                             tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCAGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCGCGCGGCGACGGCCCGACGGgcgagc                                        >  1:5831219/1‑150 (MQ=255)
                                                                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                         <  2:2582681/98‑1 (MQ=33)
                                                                                                                aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                     >  2:2182239/1‑150 (MQ=32)
                                                                                                                  tCCGGCACGCCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                   <  2:5868360/150‑1 (MQ=17)
                                                                                                                    cggCACGACCCACCATTTTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                 <  2:4225650/148‑1 (MQ=12)
                                                                                                                       cACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTc                              <  2:2562164/150‑1 (MQ=17)
                                                                                                                        ccgccccACAAATGAAAAAACGCGCCCAAACGGCGCTTTTTAACTATATGCGATACTCAAACTTCCTACCCTGCCGCCGTTTTGTTTTGCGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                             <  1:5037501/146‑1 (MQ=2)
                                                                                                                          gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATc                           <  2:5583484/150‑1 (MQ=32)
                                                                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                      >  1:5225102/1‑150 (MQ=32)
                                                                                                                                 caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCTCg                    >  2:6585401/1‑150 (MQ=17)
                                                                                                                                 caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                    >  2:6502725/1‑150 (MQ=32)
                                                                                                                                     gTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc                <  2:4927904/150‑1 (MQ=32)
                                                                                                                                      tAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca               <  2:1326146/150‑1 (MQ=32)
                                                                                                                                       aaaaaaGCGCCTTTAGGGCGCTGTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac              >  1:4500317/1‑150 (MQ=255)
                                                                                                                                       aaaaaaGCGCCCTAACGGCTATTTTTTACTATCTGCGATACTCAAATATACGAACCTGCAGCTGTTTTGTTTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac              <  1:3428564/150‑1 (MQ=11)
                                                                                                                                         aaaaGCGCCCTCACGGCGATTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca            <  2:1974443/150‑1 (MQ=11)
                                                                                                                                         aaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca            >  1:474062/1‑150 (MQ=32)
                                                                                                                                           aaGCGCCCTAACGGCGATTTTTTACTATCTGCGATACTAAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt          <  1:3819139/150‑1 (MQ=255)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt          >  2:5749843/1‑150 (MQ=32)
                                                                                                                                           aaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt          >  2:2426490/1‑150 (MQ=32)
                                                                                                                                            aGCGCCCAAACGGCGTTTTTTTACTATCTGCTATCCACAAAGATGCCAACCTGCCGCCGTTTTGTGTTGCGCGCAGCCAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATg         <  1:1510427/150‑1 (MQ=11)
                                                                                                                                              cccccTACAGGCTATTTTTTACTATCTGCGATACTAAAAGTTTCGACCCTGCTGCCGTTTTGTGTTGCGCGCCGCCAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTa       <  1:5223330/148‑1 (MQ=11)
                                                                                                                                               gcCCTACCGGCGCTATTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa      <  2:3568561/150‑1 (MQ=12)
                                                                                                                                                  cTAAAGGCGCTTTTTTACTATCTGCGATACACAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGACGCGAGTCATCCTGCACGACCCACCAAAGTaaaa  >  2:4755051/1‑149 (MQ=12)
                                                                                                                                                   tAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGTTTCGAACCTGCAGCTGTTTTGAGTTGTGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa  <  2:7868/150‑1 (MQ=255)
                                                                                                                                                   tAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATGTaaaaaa  >  2:3099590/1‑150 (MQ=17)
                                                                                                                                                   tAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa  >  2:5149007/1‑150 (MQ=32)
                                                                                                                                                    |                                                                                                                                                    
ACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAA  >  NC_000913/2521198‑2521494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: