Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,973,111 G→A 17.7% Q77* (CAG→TAG)  cheW ← purine‑binding chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,973,1110GA17.7% 187.4 / 33.9 79Q77* (CAG→TAG) cheWpurine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base A (8/6);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CGTCACGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTA  >  NC_000913/1972965‑1973255
                                                                                                                                                  |                                                                                                                                                
cGTCACGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGc                                                                                                                                                <  2:473110/149‑1 (MQ=255)
   cACGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGa                                                                                                                                             <  1:676185/149‑1 (MQ=255)
    aCGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGaa                                                                                                                                            <  2:30805/149‑1 (MQ=255)
    aCGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGaa                                                                                                                                            <  4:522257/149‑1 (MQ=255)
    aCGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGaa                                                                                                                                            >  2:287386/1‑149 (MQ=255)
     cGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAAc                                                                                                                                           <  3:71276/149‑1 (MQ=255)
       gCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACtt                                                                                                                                         >  4:309531/1‑149 (MQ=255)
            ttCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTc                                                                                                                                    <  3:624009/149‑1 (MQ=255)
                ggTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTaa                                                                                                                                <  1:244334/149‑1 (MQ=255)
                ggTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTaa                                                                                                                                <  1:375910/149‑1 (MQ=255)
                ggTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTaa                                                                                                                                <  3:217963/149‑1 (MQ=255)
                ggTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTaa                                                                                                                                <  3:309531/149‑1 (MQ=255)
                ggTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTaa                                                                                                                                <  1:168097/149‑1 (MQ=255)
                     cGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTcaa                                                                                                                           >  1:693612/1‑149 (MQ=255)
                      ggACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAAt                                                                                                                                      <  4:504163/137‑1 (MQ=255)
                      ggACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAAt                                                                                                                                      >  3:504163/1‑137 (MQ=255)
                          gAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATc                                                                                                                      <  3:339127/149‑1 (MQ=255)
                          gAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATc                                                                                                                      >  2:198735/1‑149 (MQ=255)
                             tttGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGc                                                                                                                   <  1:638711/149‑1 (MQ=255)
                             tttGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGc                                                                                                                   >  2:541827/1‑149 (MQ=255)
                              ttGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGca                                                                                                                  <  3:241769/149‑1 (MQ=255)
                               tGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGcac                                                                                                                 <  2:588526/149‑1 (MQ=255)
                                  tCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAAt                                                                                                              >  1:634137/1‑149 (MQ=255)
                                           aaTGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCcgcg                                                                                                     >  2:215961/1‑149 (MQ=255)
                                            aTGAAAGCACGTCTGAGACGCCGTAAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATGACGCCgcgc                                                                                                    >  3:583193/1‑149 (MQ=255)
                                             tGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCa                                                                                                   >  4:164024/1‑149 (MQ=255)
                                                aaGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGAt                                                                                                >  3:535897/1‑149 (MQ=255)
                                                      gTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGa                                                                                          >  4:311348/1‑149 (MQ=255)
                                                      gTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGa                                                                                          <  2:693612/149‑1 (MQ=255)
                                                                ccGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGAt                                                                                >  1:638065/1‑149 (MQ=255)
                                                                ccGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGAt                                                                                >  1:308540/1‑149 (MQ=255)
                                                                    cAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAAc                                                                            <  2:277615/149‑1 (MQ=255)
                                                                     aaCCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTAGTGACGCCTTTGATAAACg                                                                           >  4:470975/1‑149 (MQ=255)
                                                                      aCCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGc                                                                          <  2:308540/149‑1 (MQ=255)
                                                                               ccGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAAc                                                                            >  4:468392/1‑138 (MQ=255)
                                                                               ccGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAAc                                                                            <  3:468392/138‑1 (MQ=255)
                                                                               ccGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGtt                                                                 <  3:45208/149‑1 (MQ=255)
                                                                               ccGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGtt                                                                 <  3:311348/149‑1 (MQ=255)
                                                                                cGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTc                                                                <  3:541476/149‑1 (MQ=255)
                                                                                 gACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCg                                                               >  1:626177/1‑149 (MQ=255)
                                                                                 gACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCg                                                               >  3:411873/1‑149 (MQ=255)
                                                                                 gACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCttt                                                                                   <  2:576952/129‑1 (MQ=255)
                                                                                 gACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCttt                                                                                   >  1:576952/1‑129 (MQ=255)
                                                                                  accaccCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGc                                                              >  1:426962/1‑149 (MQ=255)
                                                                                    caccCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCaa                                                            <  1:198735/149‑1 (MQ=255)
                                                                                     accCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGAt                                                                                >  2:233441/1‑128 (MQ=255)
                                                                                     accCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGAt                                                                                <  1:233441/128‑1 (MQ=255)
                                                                                       ccGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCttt                                                                                   <  4:37222/123‑1 (MQ=255)
                                                                                       ccGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCttt                                                                                   >  3:37222/1‑123 (MQ=255)
                                                                                       ccGCTGTCAGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACCCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATcc                                                         <  4:369166/149‑1 (MQ=255)
                                                                                               cGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGACTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTAcc                                                 <  2:634137/149‑1 (MQ=255)
                                                                                                gagaTTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCt                                                >  2:591459/1‑149 (MQ=255)
                                                                                                gagaTTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCt                                                >  4:631891/1‑149 (MQ=255)
                                                                                                gagaTTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCt                                                >  2:139934/1‑149 (MQ=255)
                                                                                                    ttCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATc                                            <  4:535897/149‑1 (MQ=255)
                                                                                                       aGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTa                                         <  1:411128/149‑1 (MQ=255)
                                                                                                         gACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGc                                       >  4:199153/1‑149 (MQ=255)
                                                                                                                aCTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGGTCGTAGCCACGGAt                                <  3:470975/149‑1 (MQ=255)
                                                                                                                aCTACCGTGTTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGAt                                >  4:321305/1‑149 (MQ=255)
                                                                                                                  tACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGAtct                              >  3:41459/1‑149 (MQ=255)
                                                                                                                  tACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGAtct                              >  3:525183/1‑149 (MQ=255)
                                                                                                                  tACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGAtct                              >  4:358678/1‑149 (MQ=255)
                                                                                                                  tACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTCTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGAtct                              >  1:119069/1‑149 (MQ=255)
                                                                                                                      gtgtTGTCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTg                          <  1:287386/149‑1 (MQ=255)
                                                                                                                         ttgtCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTc                                                                >  4:29595/1‑108 (MQ=255)
                                                                                                                         ttgtCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTc                                                                <  3:29595/108‑1 (MQ=255)
                                                                                                                         ttgtCGTTATAGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCAc                       >  4:49580/1‑149 (MQ=255)
                                                                                                                           gtCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACtt                     >  3:372871/1‑149 (MQ=255)
                                                                                                                              gTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGCTCGTAGCCACGGATCTCCTGCACTTTCa                  <  1:312890/149‑1 (MQ=255)
                                                                                                                               ttATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAg                 >  3:361413/1‑149 (MQ=255)
                                                                                                                               ttATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAg                 <  3:203670/149‑1 (MQ=255)
                                                                                                                                   aGTCCACATCCACCTAGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGata             <  1:541827/149‑1 (MQ=255)
                                                                                                                                    gTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGatat            <  3:276192/149‑1 (MQ=255)
                                                                                                                                       cacaTCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCaa         >  1:95417/1‑149 (MQ=255)
                                                                                                                                          aTCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATAc      >  3:206078/1‑149 (MQ=255)
                                                                                                                                           tCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAg                 <  2:112527/137‑1 (MQ=255)
                                                                                                                                           tCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAg                 >  1:112527/1‑137 (MQ=255)
                                                                                                                                              aCCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTa  >  3:355574/1‑149 (MQ=255)
                                                                                                                                              aCCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTa  <  2:626177/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGTCACGGCAAATTCCGGTGCCGGACGAATTTGCTCCGCCGTCAATGAAAGCACGTCTGAGACGCCGTCAACCACGATGCCGACCACCCGCTGTCCGAGATTCAGGACGATAACTACCGTGTTGTCGTTATAGTCCACATCCACCTGGCTGAACTTAATTCGTAAGTCAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTA  >  NC_000913/1972965‑1973255

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