Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,212,482 G→T 5.6% intergenic (‑178/+58) ymgK ← / ← ymgL protein YmgK/protein YmgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,4820GT 5.6% 511.8 / 11.2 150intergenic (‑178/+58)ymgK/ymgLprotein YmgK/protein YmgL
Reads supporting (aligned to +/- strand):  ref base G (67/74);  new base T (2/7);  total (69/81)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 154 total aligned reads displayed.

AATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTC  >  NC_000913/1212333‑1212631
                                                                                                                                                     |                                                                                                                                                     
aatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACgg                                                                                                                                                       >  2:1096682/1‑150 (MQ=255)
  tGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGcc                                                                                                                                                     <  2:1310251/150‑1 (MQ=255)
    aaTATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgt                                                                                                                                                   >  2:1259520/1‑150 (MQ=255)
     atatAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtt                                                                                                                                                  <  2:476976/150‑1 (MQ=255)
     atatAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACAAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtt                                                                                                                                                  <  2:1003746/150‑1 (MQ=255)
      tataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATTATTTAAAATGTTCCTGGGTAAAATTATTTTTACCGGATTTTTAATAAATAAAAATAAATGTGCGTTAAGTATTTAAACGGCCGTGCCACACGTCCgttt                                                                                                                                                 <  2:669798/150‑2 (MQ=255)
       ataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttc                                                                                                                                                <  2:407462/150‑3 (MQ=255)
        taAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttca                                                                                                                                               <  1:597445/150‑4 (MQ=255)
        taAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgttgtt                                                                                                                                               >  1:590420/1‑150 (MQ=255)
         aaGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttt                                                                                                                                              <  2:1018519/150‑1 (MQ=255)
          aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                             <  1:665681/150‑1 (MQ=255)
          aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                             <  2:589899/150‑1 (MQ=255)
             atgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtat                                                                                                                                          <  1:1225558/150‑1 (MQ=255)
              tgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtata                                                                                                                                         <  1:990349/150‑1 (MQ=255)
               gacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaa                                                                                                                                        >  2:1068208/1‑150 (MQ=255)
                    gagaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATg                                                                                                                                   <  2:1069849/150‑1 (MQ=255)
                     agaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCATAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGt                                                                                                                                  <  1:733298/150‑1 (MQ=255)
                      gaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGtt                                                                                                                                 >  2:258029/1‑150 (MQ=255)
                        aTAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTaa                                                                                                                               <  1:948744/150‑1 (MQ=255)
                           aGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCg                                                                                                                            <  2:19133/150‑1 (MQ=255)
                            ggTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTTTAAATATGTTAACCGt                                                                                                                           <  2:1081184/150‑1 (MQ=255)
                              tttCAATTTTTGAGTTATATAGGAATGATTTAACCTGTGCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTta                                                                                                                         <  1:1007339/150‑1 (MQ=255)
                                 cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                                      <  1:48803/150‑1 (MQ=255)
                                  aaTTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaaa                                                                                                                     >  2:685017/1‑150 (MQ=255)
                                        tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                               <  2:708311/150‑1 (MQ=255)
                                            ttATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATTTTATAACggtta                                                                                                           <  1:328309/150‑5 (MQ=255)
                                             tatataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCa                                                                                                          <  1:413744/150‑1 (MQ=255)
                                             tatataGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCa                                                                                                          >  2:1275920/1‑150 (MQ=255)
                                                 taGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAg                                                                                                      <  1:960084/150‑1 (MQ=255)
                                                   ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                                    <  1:900548/150‑1 (MQ=255)
                                                   ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                                    >  1:627109/1‑150 (MQ=255)
                                                   ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                                    >  1:558092/1‑150 (MQ=255)
                                                   ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATAAATGTGCGTTAAGTATTGAAACGGACTTGTGGCAAGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGta                                                                                                    >  1:558093/1‑149 (MQ=255)
                                                       tGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGt                                                                                                <  1:1373880/150‑1 (MQ=255)
                                                        gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGGTGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                               >  2:1338303/1‑150 (MQ=255)
                                                         aTTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGa                                                                                              >  1:1400412/1‑150 (MQ=255)
                                                           ttAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACTTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATc                                                                                            <  2:825724/150‑1 (MQ=255)
                                                           ttAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATc                                                                                            <  1:1258214/150‑1 (MQ=255)
                                                            tAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATca                                                                                           <  1:1096682/150‑1 (MQ=255)
                                                                       tGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAata                                                                                <  2:859524/150‑1 (MQ=255)
                                                                        ggCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAatat                                                                               >  2:57365/1‑150 (MQ=255)
                                                                          cTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCa                                                                             >  2:382216/1‑150 (MQ=255)
                                                                            aaaaTACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAg                                                                           <  1:470130/150‑1 (MQ=255)
                                                                            aaaaTACATATAACCGGATGATGACTAAACCAAAATAAATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAg                                                                           <  1:382216/150‑1 (MQ=255)
                                                                                tACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGa                                                                       >  1:537823/1‑150 (MQ=255)
                                                                                  cATATAACCGGATGATGACTAAACAAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGATTTTATAAATATGTTAACAGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTCATa                                                                     >  1:444883/1‑150 (MQ=255)
                                                                                   atatAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAg                                                                    >  1:1420226/1‑150 (MQ=255)
                                                                                    tataACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGt                                                                   >  2:37227/1‑150 (MQ=255)
                                                                                    tataACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGt                                                                   >  1:873712/1‑150 (MQ=255)
                                                                                      taACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGt                                                                                                              <  1:224953/105‑1 (MQ=255)
                                                                                        aCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAATAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTa                                                               >  2:1239815/1‑150 (MQ=255)
                                                                                        aCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAata                                                                                <  2:319668/133‑1 (MQ=255)
                                                                                          cGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTAGCAAAAGTCAAGTGCTCCCAGTTCACATATCAAGTGGATAGTATTAgg                                                             >  1:1038032/1‑149 (MQ=255)
                                                                                              tgatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTgg                                                         <  1:872921/150‑1 (MQ=255)
                                                                                                atgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTg                                                       <  2:875030/150‑1 (MQ=255)
                                                                                                 tgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGa                                                      >  1:923035/1‑150 (MQ=255)
                                                                                                  gaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGat                                                     >  1:281141/1‑150 (MQ=255)
                                                                                                  gaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGat                                                     >  2:1169436/1‑150 (MQ=255)
                                                                                                   aCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGatt                                                    >  1:392526/1‑150 (MQ=255)
                                                                                                      aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                                 <  1:614947/150‑1 (MQ=255)
                                                                                                      aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                                 >  2:614947/1‑150 (MQ=255)
                                                                                                        aCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttattta                                               >  1:240873/1‑150 (MQ=255)
                                                                                                         ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                              <  2:627109/150‑1 (MQ=255)
                                                                                                          cAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATg                                             >  1:186459/1‑150 (MQ=255)
                                                                                                           aaaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                            >  2:1243041/1‑150 (MQ=255)
                                                                                                              aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTg                                         <  2:806164/150‑1 (MQ=255)
                                                                                                              aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTg                                         <  2:285906/150‑1 (MQ=255)
                                                                                                               tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTta                                                                                                                                           <  1:1374909/51‑1 (MQ=255)
                                                                                                               tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                            <  1:286361/146‑1 (MQ=255)
                                                                                                               tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                            >  2:286361/1‑146 (MQ=255)
                                                                                                                  aTGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACaa                                     >  1:889239/1‑150 (MQ=255)
                                                                                                                    gtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATa                                   <  1:280367/150‑1 (MQ=255)
                                                                                                                      gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                 >  2:236588/1‑150 (MQ=255)
                                                                                                                      gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                 >  1:812633/1‑150 (MQ=255)
                                                                                                                      gCGTTAAGTATTGAAACGGACGTGTGACACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                 >  2:1370076/1‑150 (MQ=255)
                                                                                                                             gTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGt                                                                                                                           >  2:663665/1‑53 (MQ=255)
                                                                                                                             gTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGt                                                                                                                           <  1:663665/53‑1 (MQ=255)
                                                                                                                             gTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAg                          <  2:537823/150‑1 (MQ=255)
                                                                                                                               aTTGAAACGGACGTGTGGCACTGCCGATGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGCTGATTATTTATGGGTGACAATAAAAAGACAGTa                        >  2:393848/1‑150 (MQ=255)
                                                                                                                               aTTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTa                        >  2:240490/1‑150 (MQ=255)
                                                                                                                                ttGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAt                       >  2:634213/1‑150 (MQ=255)
                                                                                                                                 tGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAtt                      >  2:1381993/1‑150 (MQ=255)
                                                                                                                                  gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                     <  1:1275920/150‑1 (MQ=255)
                                                                                                                                  gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                     >  1:544360/1‑150 (MQ=255)
                                                                                                                                  gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATAAGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                     >  2:252607/1‑150 (MQ=255)
                                                                                                                                   aaaCGGACGTGTGGCACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTaa                    <  2:1184913/150‑1 (MQ=255)
                                                                                                                                   aaaCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTaa                    <  2:889239/150‑1 (MQ=255)
                                                                                                                                     aCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatc                  >  2:517337/1‑150 (MQ=255)
                                                                                                                                      cgtaCGTGTGGCACGGCCGTTGTTTTTATAAATATTTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatca                 <  1:930247/147‑1 (MQ=255)
                                                                                                                                       ggACGTGTGGCCCGGCCGTTGTTTTTATAAATATCTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcat                >  2:451233/1‑150 (MQ=255)
                                                                                                                                         aCGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATcc              <  1:1338303/150‑1 (MQ=255)
                                                                                                                                          cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa             <  2:1365869/150‑1 (MQ=255)
                                                                                                                                          cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa             <  2:765795/150‑1 (MQ=255)
                                                                                                                                            tgtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATa           >  2:309436/1‑150 (MQ=255)
                                                                                                                                             gtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTATGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATag          >  2:547969/1‑150 (MQ=255)
                                                                                                                                               ggCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATagag        >  1:1002634/1‑150 (MQ=255)
                                                                                                                                                gCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATAAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAgaga       >  1:309020/1‑150 (MQ=255)
                                                                                                                                                 cACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGAt      >  2:13276/1‑150 (MQ=255)
                                                                                                                                                     gCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGtc  >  1:7888/1‑150 (MQ=255)
                                                                                                                                                     gCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGtc  <  2:946381/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTC  >  NC_000913/1212333‑1212631

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