Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,415,523 G→C 9.3% E417Q (GAG→CAG) ‡ aidB → DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,415,5230GC 9.3% 378.5 / 28.7 129E417Q (GAG→CAG) ‡aidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (56/60);  new base C (8/4);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 139 total aligned reads displayed.

ATGCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATCGCGCG  >  NC_000913/4415374‑4415672
                                                                                                                                                     |                                                                                                                                                     
 tGCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAACGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGga                                                                                                                                                      <  2:16323/150‑1 (MQ=255)
  gCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAATAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgag                                                                                                                                                     >  2:1641504/1‑150 (MQ=255)
  gCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgag                                                                                                                                                     >  1:499627/1‑150 (MQ=255)
  gCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgag                                                                                                                                                     <  1:569900/150‑1 (MQ=255)
    cAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGAt                                                                                                                                                   >  1:1018328/1‑150 (MQ=255)
    cAAAGAAGCCCTGCGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGAt                                                                                                                                                   <  1:1307652/150‑1 (MQ=255)
      aagaagCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGc                                                                                                                                                 >  1:168102/1‑150 (MQ=255)
      aagaagCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGc                                                                                                                                                 >  1:746294/1‑150 (MQ=255)
        gaagCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCg                                                                                                                                               >  1:573481/1‑150 (MQ=255)
         aagCCCTGTGGGCGCGTTTATTTACGCCGGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCgg                                                                                                                                              <  1:857017/150‑1 (MQ=255)
          atcccTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGGGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGt                                                                                                                                             <  1:925215/148‑1 (MQ=255)
           gCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTa                                                                                                                                            >  2:690989/1‑150 (MQ=255)
              cTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaag                                                                                                                                         <  1:1414547/150‑2 (MQ=255)
               tgtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                        <  1:1642673/150‑3 (MQ=255)
               tgtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGGTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                        <  1:667606/150‑3 (MQ=255)
               tgtgGGCGCGTCTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGGGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                        <  1:450358/150‑3 (MQ=255)
                gtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAg                                                                                                                                       <  1:1024666/150‑1 (MQ=255)
                  gggCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCTGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTa                                                                                                                                     >  2:956551/1‑150 (MQ=255)
                  gggCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTTCCGGGAGATGCCGGTAAACAGTa                                                                                                                                     >  2:341192/1‑150 (MQ=255)
                  gggCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCAGTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTa                                                                                                                                     >  2:1463685/1‑150 (MQ=255)
                       cgTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGCATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTgg                                                                                                                                <  1:1251831/150‑1 (MQ=255)
                         tttatttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGCAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGa                                                                                                                              >  1:1668279/1‑150 (MQ=255)
                         tttatttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGa                                                                                                                              >  2:49263/1‑150 (MQ=255)
                            atttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAgg                                                                                                                           <  1:1875660/150‑1 (MQ=255)
                            atttaCGCCTGCGGCAAAATTTCTGATCTCCAAAAGAGGTATGCCGCTTGTGGCCGAAGCGAGGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAgg                                                                                                                           <  2:252138/150‑1 (MQ=255)
                                 cGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAACCAGGATTTGGGATGGTTCCg                                                                                                                      >  2:86968/1‑150 (MQ=255)
                                  gCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCgg                                                                                                                     <  2:1485376/150‑1 (MQ=255)
                                  cccTGCGGCAAAATTTGTGATCTTCAAAAGCGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCGGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCgg                                                                                                                     <  2:648397/149‑1 (MQ=255)
                                  cacTGCGGAAAAATTTTGTATCTGCAAAAGAGGTATGACGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGTGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGTGAAGGTTCCgg                                                                                                                     <  2:655053/148‑1 (MQ=255)
                                   ccTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGTCATTTGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGc                                                                                                                    >  1:379452/1‑150 (MQ=255)
                                      gCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAat                                                                                                                 <  1:1286694/150‑1 (MQ=255)
                                             aaTTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtgt                                                                                                          >  2:1565215/1‑150 (MQ=255)
                                             aaTTTGTGATATGCAAAAGAGGTATGCCGTTTGTGGCCGTAGCGTTGGAGGTGCTGTGCTGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGTCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtgt                                                                                                          <  1:322077/150‑1 (MQ=255)
                                             aaTTTGGGATCTGCAAATGAGGTATGCCGGTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtgt                                                                                                          <  1:629458/150‑1 (MQ=255)
                                               tttGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGc                                                                                                        <  2:1455280/150‑1 (MQ=255)
                                                  gtgATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTg                                                                                                     >  1:398977/1‑150 (MQ=255)
                                                    gATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGa                                                                                                   >  2:1852605/1‑150 (MQ=255)
                                                        tGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtgtg                                                                                               >  1:1593015/1‑150 (MQ=255)
                                                        tGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtgtg                                                                                               >  1:844370/1‑150 (MQ=255)
                                                           aaaaGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTg                                                                                            <  2:1093888/150‑1 (MQ=255)
                                                            aaaGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCTAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTgc                                                                                           >  2:96208/1‑150 (MQ=255)
                                                             aaGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTgcg                                                                                          <  1:389269/150‑1 (MQ=255)
                                                               tagGTATGCCGTTGGTGGCCGAAGCGATGGATGTGCTGGGCGTGAATGGTTATTGCGATGAGAGTCAGCTGCCGCGTCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGcgcg                                                                                        <  1:330443/149‑1 (MQ=255)
                                                                  gTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtc                                                                                     <  2:973888/150‑1 (MQ=255)
                                                                   tATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtct                                                                                    >  1:1186/1‑150 (MQ=255)
                                                                   tATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtct                                                                                    <  1:1386253/150‑1 (MQ=255)
                                                                    ctGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtctc                                                                                   <  2:1838279/149‑1 (MQ=255)
                                                                     tgtcGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCa                                                                                  <  1:14252/147‑1 (MQ=255)
                                                                     tGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTATACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCa                                                                                  <  1:1045643/150‑1 (MQ=255)
                                                                     tGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCa                                                                                  >  1:825660/1‑150 (MQ=255)
                                                                       ccGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGCGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAAt                                                                                >  2:1190293/1‑150 (MQ=255)
                                                                       ccGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAAt                                                                                >  1:1516711/1‑150 (MQ=255)
                                                                           ttGTGGCCGCTGCGATCGAGGTGCTGGGCGGGATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCAGGTAAACATTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagc                                                                            <  2:1251760/150‑1 (MQ=255)
                                                                            tgtgGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagca                                                                           >  2:1722523/1‑150 (MQ=255)
                                                                            tgtgGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGAGAGCGGCCGCGGCGTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagca                                                                           <  1:794501/150‑1 (MQ=255)
                                                                               ggCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagcaagc                                                                        >  2:1646332/1‑150 (MQ=255)
                                                                                gCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCg                                                                       <  1:1483098/150‑1 (MQ=255)
                                                                                 ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTTCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCgg                                                                      >  1:1108063/1‑150 (MQ=255)
                                                                                 ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCgg                                                                      >  1:1269667/1‑150 (MQ=255)
                                                                                   gAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGGGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGc                                                                    >  2:666371/1‑150 (MQ=255)
                                                                                   gAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAAGAAGCAAGCGGGc                                                                    >  1:1169340/1‑150 (MQ=255)
                                                                                    aaGCGAGGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGGCTCAATAAGCAAGCGGGCg                                                                   >  2:1479964/1‑150 (MQ=255)
                                                                                         aTGGAGGTGCTGGGCGGCATTGGGTATTGCGAGGAGAGGGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGGTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATAc                                                              <  2:755070/150‑1 (MQ=255)
                                                                                           ggAGGTGCAGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGa                                                            >  2:1853149/1‑150 (MQ=255)
                                                                                              ggTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACtt                                                         <  2:932799/150‑1 (MQ=255)
                                                                                               gTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTa                                                        <  1:607299/150‑1 (MQ=255)
                                                                                                tGCTGGGCGGCATTGGTTATTGCGAGGAGTGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTAt                                                       <  2:1513151/150‑1 (MQ=255)
                                                                                                 gCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTAtt                                                      >  2:1074577/1‑150 (MQ=255)
                                                                                                  ccgggCGGCATTGGTTATTGCGAGGAGAGTGAGCAGCCGCGGCTTTACCTTTATATGCCGGTAAACAGTATTTGGGAATGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTg                                                     <  2:38295/148‑1 (MQ=255)
                                                                                                        ggcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGaa                                               <  1:1241358/150‑1 (MQ=255)
                                                                                                                 tATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTtgtg                                      <  1:1070799/150‑1 (MQ=255)
                                                                                                                  aTTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTgg                                     >  2:1148998/1‑150 (MQ=255)
                                                                                                                     gCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAg                                  >  1:1773168/1‑150 (MQ=255)
                                                                                                                       gaggagAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTAGGGAAGGTTCCGGCAATATTATGTGCGTGGATGTGTTGCGCGTTCTCAATCAGCAAGCGGGCGTATACGCCTTATTGTCGGAAGCATTTGTGGAAGTg                                >  2:1131425/1‑150 (MQ=255)
                                                                                                                           agagTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAg                            <  2:573481/150‑1 (MQ=255)
                                                                                                                             agTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAggg                          >  1:546762/1‑150 (MQ=255)
                                                                                                                              gTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagca                                                                           <  2:405068/100‑1 (MQ=255)
                                                                                                                              gTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagca                                                                           >  1:405068/1‑100 (MQ=255)
                                                                                                                               tGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCa                        >  2:194169/1‑150 (MQ=255)
                                                                                                                                 aGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAgg                      >  1:889163/1‑150 (MQ=255)
                                                                                                                                     gCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCg                  >  2:1363515/1‑150 (MQ=255)
                                                                                                                                     gCCGCGGCGTTACCGGGAGCTGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCg                  <  2:64514/150‑1 (MQ=255)
                                                                                                                                      ccGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGGGAAAGGGCAGGCTCGc                 <  2:1269667/150‑1 (MQ=255)
                                                                                                                                        gcgGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTa               >  1:1464969/1‑150 (MQ=255)
                                                                                                                                          ggCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt             <  2:398977/150‑1 (MQ=255)
                                                                                                                                          ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGATGTGAAAGGGCAGGATCGCTAtt             >  1:1091218/4‑150 (MQ=255)
                                                                                                                                          ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt             >  2:1340802/4‑150 (MQ=255)
                                                                                                                                           tgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAttt            >  2:1048422/3‑150 (MQ=255)
                                                                                                                                           gcgttACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAttt            <  2:60844/147‑1 (MQ=255)
                                                                                                                                            gtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATTTTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTTTTGTCGGAAGCATTTGTGGAAGTGAAAGGGCTGGATCGCTAtttt           >  1:371962/2‑150 (MQ=255)
                                                                                                                                            gtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtttt           >  1:87322/2‑150 (MQ=255)
                                                                                                                                            cTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtttt           <  2:1018328/150‑1 (MQ=255)
                                                                                                                                             tttACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTg          <  1:327519/150‑1 (MQ=255)
                                                                                                                                                agccggAGCTGCCCGTAAACAGGATGTGGGAAGCCGCCGGCAATATTAAGTCCCGGGATGTGAAGCGCGTTCTCAATAAGCAATCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGGGAAAGGGCAGGATCGCTATTATGATc       <  2:354320/146‑1 (MQ=255)
                                                                                                                                                aCCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCTATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCTAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATc       <  2:1353400/150‑1 (MQ=255)
                                                                                                                                                aCCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATc       <  2:408918/150‑1 (MQ=255)
                                                                                                                                                aCCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATc       <  1:1565215/150‑1 (MQ=255)
                                                                                                                                                 ccGGGAGATGCCGGTAAACAGTATTTGGGCAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTTGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcg      <  2:499627/150‑1 (MQ=255)
                                                                                                                                                 ccGGGAGATGCCGGTAAACAGTAATTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcg      >  2:899406/1‑150 (MQ=255)
                                                                                                                                                     gagaTGCCGGTAAACAGTATATGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcgcgcg  <  2:605968/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATCGCGCG  >  NC_000913/4415374‑4415672

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