Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,016 T→C 5.4% intergenic (+32/‑28) secA → / → mutT preprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,0160TC 5.4% 568.6 / 11.1 173intergenic (+32/‑28)secA/mutTpreprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (87/76);  new base C (3/7);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

TCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAAT  >  NC_000913/110868‑111162
                                                                                                                                                    |                                                                                                                                                  
tctGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                                                                    >  2:468100/1‑149 (MQ=255)
tctGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAtt                                                                                                                                                    >  4:104472/1‑149 (MQ=255)
  tGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCt                                                                                                                                                  >  2:630100/1‑149 (MQ=255)
  tGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCt                                                                                                                                                  >  1:1668125/1‑149 (MQ=255)
    cAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                                                                <  2:291267/149‑1 (MQ=255)
    cAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                                                                <  4:837227/149‑1 (MQ=255)
      gCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcgc                                                                                                                                              >  4:752522/1‑149 (MQ=255)
       ccccAGCTGCACTGGCGGCGCAAACCTGAGAGCACAAAGTAGCACGTAACGATCCTTGCCCGTGCCGTTCTGGTAAAAAATACAAGCAGTGCTATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGcc                                                                                                                                             >  4:1383438/4‑149 (MQ=255)
        cGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                            <  4:793253/149‑1 (MQ=255)
         gCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtt                                                                                                                                           <  4:1095252/149‑1 (MQ=255)
             cTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCtttttt                                                                                                                                       <  4:552051/149‑1 (MQ=255)
               gCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtat                                                                                                                                     >  4:292408/1‑149 (MQ=255)
                cACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtata                                                                                                                                    <  3:289959/149‑1 (MQ=255)
                cACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtata                                                                                                                                    <  3:590181/149‑1 (MQ=255)
                 aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                   <  4:1432159/149‑1 (MQ=255)
                 aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                   <  4:126968/149‑1 (MQ=255)
                   tGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttt                                                                                                                                          <  2:718146/140‑1 (MQ=255)
                   tGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttt                                                                                                                                          >  1:718146/1‑140 (MQ=255)
                          gcAAACCGGAGGGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                                                                          >  3:1009047/1‑149 (MQ=255)
                          gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                                                                          >  4:1638362/1‑149 (MQ=255)
                          gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                                                                          >  4:342465/1‑149 (MQ=255)
                            aaaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                                                                   >  3:725563/1‑122 (MQ=255)
                             aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAAt                                                                                                                       <  1:1082356/149‑1 (MQ=255)
                                  gagaGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaa                                                                                                                  >  4:227766/1‑149 (MQ=255)
                                   agagCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaaa                                                                                                                 >  1:858093/1‑149 (MQ=255)
                                     agCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGc                                                                                                               <  2:1668125/149‑1 (MQ=255)
                                       cgcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTg                                                                                                             <  1:1163436/149‑1 (MQ=255)
                                          aaaGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTaa                                                                                                                             <  1:1534167/130‑1 (MQ=255)
                                           aaGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGATAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCaaa                                                                                                         <  3:1196735/149‑1 (MQ=255)
                                                ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                    <  3:1320087/149‑1 (MQ=255)
                                                ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                    <  3:240109/149‑1 (MQ=255)
                                                    gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                                                                <  3:388932/149‑1 (MQ=255)
                                                    gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtac                                                                                                <  2:1491181/149‑2 (MQ=255)
                                                    gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                            >  1:1022017/1‑105 (MQ=255)
                                                    gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                            <  2:1022017/105‑1 (MQ=255)
                                                         gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatt                                                                                           >  4:1561171/1‑149 (MQ=255)
                                                          aTCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatta                                                                                          >  2:153005/1‑149 (MQ=255)
                                                             cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                       <  2:311382/149‑1 (MQ=255)
                                                              ttGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                   >  3:53706/1‑104 (MQ=255)
                                                              ttGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGc                                                                                                               <  1:201512/124‑1 (MQ=255)
                                                                gCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCa                                                                                    >  2:1538089/1‑149 (MQ=255)
                                                                       cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGaa                                                                             <  4:898303/149‑1 (MQ=255)
                                                                        ggTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAAc                                                                            >  2:901445/1‑149 (MQ=255)
                                                                         gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCtttt                                                                                                                                         <  3:1665684/87‑1 (MQ=255)
                                                                          ttCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACaa                                                                          >  2:860887/1‑149 (MQ=255)
                                                                             tGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGa                                                                       >  4:901797/1‑149 (MQ=255)
                                                                               gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                                     >  3:300872/1‑149 (MQ=255)
                                                                                tAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                       >  2:1604752/1‑130 (MQ=255)
                                                                                 aaaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGTAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATc                                                                   <  3:1575064/149‑1 (MQ=255)
                                                                                  aaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCt                                                                  >  4:650712/1‑149 (MQ=255)
                                                                                   aaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCtt                                                                 <  1:1418079/149‑1 (MQ=255)
                                                                                    aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                >  4:591620/1‑149 (MQ=255)
                                                                                    aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                <  2:1075144/149‑1 (MQ=255)
                                                                                     aaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                               <  3:1561171/149‑1 (MQ=255)
                                                                                      aTACAAGCAGTGCCATTGCCGCCTGCAATGAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat                                                              <  3:366277/149‑1 (MQ=255)
                                                                                      aTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGTAATCTTtat                                                              <  3:1300261/149‑1 (MQ=255)
                                                                                      aTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAACGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                       <  2:627051/124‑1 (MQ=255)
                                                                                      aTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAACGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                       >  1:627051/1‑124 (MQ=255)
                                                                                        aCAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                    >  1:678970/1‑109 (MQ=255)
                                                                                            gCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgcg                                                        >  2:332026/1‑149 (MQ=255)
                                                                                              aGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTc                                                      >  1:993462/1‑149 (MQ=255)
                                                                                               gTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcg                                                     >  3:496214/1‑149 (MQ=255)
                                                                                                    aTGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcag                                                <  3:227766/149‑1 (MQ=255)
                                                                                                        ccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGa                                            <  2:751126/149‑1 (MQ=255)
                                                                                                            cTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgc                                                         >  3:653270/1‑132 (MQ=255)
                                                                                                                aaTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGcaca                                    >  3:1650882/1‑149 (MQ=255)
                                                                                                                  tAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATg                                  >  1:218350/1‑149 (MQ=255)
                                                                                                                   aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATgg                                 >  3:5654/1‑149 (MQ=255)
                                                                                                                    aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACTAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACGTGGc                                <  4:696568/149‑1 (MQ=255)
                                                                                                                    aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACaa                                                                          <  2:963659/107‑1 (MQ=255)
                                                                                                                    aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACaa                                                                          >  1:963659/1‑107 (MQ=255)
                                                                                                                     aaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                                                                <  4:504515/36‑1 (MQ=255)
                                                                                                                           aatGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATaaa                         >  1:1109478/3‑149 (MQ=255)
                                                                                                                            cTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTGTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAAc                        <  3:1406631/149‑1 (MQ=255)
                                                                                                                              gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTg                      >  2:1292169/1‑149 (MQ=255)
                                                                                                                               ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTgg                     <  2:215134/149‑1 (MQ=255)
                                                                                                                               ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTgg                     >  1:1661948/1‑149 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCATGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  2:1033694/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAGAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  1:519165/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  2:1345127/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  4:1486649/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  3:752522/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  1:901445/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  2:1187860/149‑1 (MQ=255)
                                                                                                                                tGAAGTAAAAGACGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  3:1059093/149‑1 (MQ=255)
                                                                                                                                  aaGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGt                  <  1:1555533/149‑1 (MQ=255)
                                                                                                                                    gTAAAAGGCGCAGGATTCTGCGCCCTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATCATTCGCAACGAGAACAATGAAATCTATATAACGCGTCGCGCAGCAGGTGCGCACATGGCGAATAAGCTGGAGttt                >  3:1455631/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTc               >  1:1178317/2‑149 (MQ=255)
                                                                                                                                       aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc             >  3:1213185/1‑149 (MQ=255)
                                                                                                                                       aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc             >  4:262222/1‑149 (MQ=255)
                                                                                                                                       aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc             >  2:1719387/1‑149 (MQ=255)
                                                                                                                                       aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc             >  1:1320167/1‑149 (MQ=255)
                                                                                                                                        aaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcg            >  4:1706454/1‑149 (MQ=255)
                                                                                                                                          ggCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggc          <  3:591620/149‑1 (MQ=255)
                                                                                                                                            cgcAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGc                                >  1:1598348/1‑125 (MQ=255)
                                                                                                                                             gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                <  2:1174562/92‑1 (MQ=255)
                                                                                                                                             gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                >  1:1174562/1‑92 (MQ=255)
                                                                                                                                             gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTc                                                      <  3:231429/102‑1 (MQ=255)
                                                                                                                                             gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGt       >  3:81054/1‑149 (MQ=255)
                                                                                                                                                   ttCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggc          <  4:232723/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAAT  >  NC_000913/110868‑111162

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