Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,415,526 A→C 7.1% M418L (ATG→CTG) ‡ 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,5260AC 7.1% 287.2 / 11.2 99M418L (ATG→CTG) ‡aidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (47/44);  new base C (3/4);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

Only 100 of 106 total aligned reads displayed.

CCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATCGCGCGG  >  NC_000913/4415377‑4415673
                                                                                                                                                     |                                                                                                                                                   
ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgaga                                                                                                                                                     >  1:1059761/1‑150 (MQ=255)
ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgaga                                                                                                                                                     >  1:1188138/1‑150 (MQ=255)
ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgaga                                                                                                                                                     >  2:1748811/1‑150 (MQ=255)
ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgaga                                                                                                                                                     >  1:598135/1‑150 (MQ=255)
ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGgaga                                                                                                                                                     >  1:1653906/1‑150 (MQ=255)
ccAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGATATTGCGAGGAGAGTCAGCTGCAGCGGCTTTACCGGgaga                                                                                                                                                     >  2:137860/1‑150 (MQ=255)
 cAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGAt                                                                                                                                                    >  1:8884/1‑150 (MQ=255)
  aaaGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATg                                                                                                                                                   <  2:581311/150‑1 (MQ=255)
    agaagCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGcc                                                                                                                                                 <  1:1437710/150‑1 (MQ=255)
       agCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGt                                                                                                                                              <  2:1766731/150‑1 (MQ=255)
        gCCCTGTGGGCGCGTATATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTa                                                                                                                                             >  2:505845/1‑150 (MQ=255)
         cccTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTaa                                                                                                                                            >  1:1112642/1‑150 (MQ=255)
          ccTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCATAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTTTTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTaaa                                                                                                                                           >  2:138194/1‑150 (MQ=255)
           cTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAAc                                                                                                                                          <  1:1939674/150‑1 (MQ=255)
           cTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaag                                                                                                                                          <  1:783156/150‑2 (MQ=255)
            tgtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGCATCTCCCGGTaaagc                                                                                                                                         <  2:605114/150‑3 (MQ=255)
             gtgGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAg                                                                                                                                        <  2:1750099/150‑1 (MQ=255)
                   gcgTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTg                                                                                                                                  >  1:680433/1‑150 (MQ=255)
                    cgTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTgg                                                                                                                                 >  1:1437261/1‑150 (MQ=255)
                    cgTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTgg                                                                                                                                 >  1:319486/1‑150 (MQ=255)
                     gTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTggg                                                                                                                                <  1:611806/150‑1 (MQ=255)
                         atttaCGCCTGCGTCAATCTTTTTTCTCTGCAACACAGGTATGCCGTTTGTGGCCGAAGCGAGGGCGGTGCTGGTCGTCATTGGTTATTGCGAGGAGCGTGAGCTGACGCGGCTTTACCGGGAGATGCCGGTAAACAGAATTTGGGAAgg                                                                                                                            <  2:1534430/150‑1 (MQ=255)
                          tttaCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGt                                                                                                                           >  1:659870/1‑150 (MQ=255)
                            tccGCCTGCGGCAAAATTTGTGATCTGCACAAGAGGTATGCCGTTTGTGTCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTc                                                                                                                         <  1:394515/148‑1 (MQ=255)
                            taCGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTc                                                                                                                         <  2:1995232/150‑1 (MQ=255)
                               gCCGGCGGCAAAATTTGTGATCTGCACAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCGGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCTgg                                                                                                                      <  2:1362267/150‑1 (MQ=255)
                                ccTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGc                                                                                                                     <  2:960859/150‑1 (MQ=255)
                                 cTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCa                                                                                                                    >  1:262309/1‑150 (MQ=255)
                                   gCGGCAAAATTTGTGATCTGCAAAAGAGGTTTGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAat                                                                                                                  <  1:936581/150‑1 (MQ=255)
                                      gCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAAtatt                                                                                                               >  1:172278/1‑150 (MQ=255)
                                       cAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAAtatta                                                                                                              >  1:1306294/1‑150 (MQ=255)
                                       cAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAAtatta                                                                                                              >  2:571237/1‑150 (MQ=255)
                                         aaaTTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtg                                                                                                            <  2:919853/150‑1 (MQ=255)
                                         aaaTTGGTGATCTGCAAAAGAGGTATGCCGTTGGGGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtg                                                                                                            <  1:258745/150‑1 (MQ=255)
                                          aaTTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTAtgt                                                                                                           <  1:1972779/150‑1 (MQ=255)
                                            tttGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGc                                                                                                         <  1:1830290/150‑1 (MQ=255)
                                               gtgATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTg                                                                                                      <  2:1882967/150‑1 (MQ=255)
                                               gtgATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGACGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGACTg                                                                                                      >  2:358922/1‑150 (MQ=255)
                                                tgATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTgg                                                                                                     >  2:1846929/1‑150 (MQ=255)
                                                 gATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGa                                                                                                    >  1:1952828/1‑150 (MQ=255)
                                                 gATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGa                                                                                                    >  2:429064/1‑150 (MQ=255)
                                                 gATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGa                                                                                                    <  2:659870/150‑1 (MQ=255)
                                                  aTCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAt                                                                                                   >  1:866192/1‑150 (MQ=255)
                                                   tCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtg                                                                                                  <  1:675691/150‑1 (MQ=255)
                                                   tCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGAtg                                                                                                  >  1:1366671/1‑150 (MQ=255)
                                                       cAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGtt                                                                                              >  1:926738/1‑150 (MQ=255)
                                                        taaaGAGGTTTGCGGTTTGGCGCCGAAGCGATGGTGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTATACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTg                                                                                             <  1:329917/149‑1 (MQ=255)
                                                        aaaaGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCAGTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAGACAGTATTTCGGAAGGTTCCGGCAAAAATATGTGCCTGGAAgtgtat                                                                                             >  1:322585/1‑148 (MQ=255)
                                                          aaGACGTATGCCGTATGTGGCCGAAGCTATGGAGGCGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTATACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTgcg                                                                                           <  2:987887/150‑1 (MQ=255)
                                                               gTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtc                                                                                      >  1:563851/1‑150 (MQ=255)
                                                                tATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtct                                                                                     <  2:28681/150‑1 (MQ=255)
                                                                 ttGCCGGTGGTGGCCGTAGCGCTGGCGGGGGTGGGCGGGAGTGGGTTTTGCGAGGAGAGTGAGGTGCCGCGGGTTTACCGGGTGATGCCGGTAAACAGTATTTGGTATGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTtctc                                                                                    <  1:813059/149‑1 (MQ=255)
                                                                 attccGTGTGTGGCCGACTCGCGGGCGGTGCGGGGCGGCCTGGGGTATTGCGAGGAGATTGAGCTGCCGCGTCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGtctc                                                                                    <  2:1476721/147‑1 (MQ=255)
                                                                   gCCGTTTGGGGCCGAAGCGAGGGAGGTGCTGGGCGGCTTTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGTTGCCGGTAAACAGTATTTGGGAAGGATCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCaa                                                                                  <  1:1495313/150‑1 (MQ=255)
                                                                    ccGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAAt                                                                                 >  2:1556451/1‑150 (MQ=255)
                                                                       tttGTGGCCGACGCGATGGAGGGGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGACTGGATGTGTTGCGCGTTCTCAATaag                                                                              <  1:1919202/150‑1 (MQ=255)
                                                                       tttGTGGCCGAAGCGATGGAGGTGCTTGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaag                                                                              <  1:1977067/150‑1 (MQ=255)
                                                                       tttGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGTCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaag                                                                              <  1:2009514/150‑1 (MQ=255)
                                                                        ttGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagc                                                                             <  2:580024/150‑1 (MQ=255)
                                                                          gtgGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagcaa                                                                           >  2:828959/1‑150 (MQ=255)
                                                                            ggCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATaagcacgc                                                                         >  2:746787/1‑150 (MQ=255)
                                                                              ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCgg                                                                       >  2:663452/1‑150 (MQ=255)
                                                                              ccGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCgg                                                                       >  1:1766874/1‑150 (MQ=255)
                                                                                 aaGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCg                                                                    >  2:304948/1‑150 (MQ=255)
                                                                                   gCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGta                                                                  <  2:1059760/150‑1 (MQ=255)
                                                                                     gATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGtata                                                                <  1:875310/150‑1 (MQ=255)
                                                                                                   gcggcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCgg                                                  >  2:623596/1‑150 (MQ=255)
                                                                                                   gcggcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCgg                                                  >  1:886125/1‑150 (MQ=255)
                                                                                                    cggcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATGAGCAAGCGGGCGTATACGACTTATTGTCGGa                                                 <  2:1842569/150‑1 (MQ=255)
                                                                                                    cggcATTGGTTATTGCGAGGAGAGTGAGCGGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGa                                                 <  2:2044300/150‑1 (MQ=255)
                                                                                                     ggcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGaa                                                <  2:695978/150‑1 (MQ=255)
                                                                                                      gcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAg                                               >  1:1273498/1‑150 (MQ=255)
                                                                                                      gcATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTAGGAAg                                               >  2:703352/1‑150 (MQ=255)
                                                                                                      gcATTGGTTATTGAGAGTAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAg                                               <  1:137860/150‑1 (MQ=255)
                                                                                                         ttGGGTATTGCGAGGAGAGTGAGCTGCCGCGGCGTTACCGGGCGATGCCGGTAAACAGTAGTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCAt                                            <  2:551598/150‑1 (MQ=255)
                                                                                                          tGGTTATTGCGCGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTCCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCAtt                                           <  2:319486/150‑1 (MQ=255)
                                                                                                           ggTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCAttt                                          >  1:695209/1‑150 (MQ=255)
                                                                                                            gTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTAGCGGGAGATGCCGGTAAACAGTATTTGGGTAGGTTCCGGCAAGATTAAGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTtg                                         <  2:411640/150‑1 (MQ=255)
                                                                                                              tATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTtgtg                                       >  2:876494/1‑150 (MQ=255)
                                                                                                              tATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTtgtg                                       <  2:1915459/150‑1 (MQ=255)
                                                                                                                ttGCGAGGAGAGTGAGATGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGTGAAGGTTCCGGCTATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGa                                     >  2:1709525/1‑150 (MQ=255)
                                                                                                                  gCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAg                                   >  2:1666583/1‑150 (MQ=255)
                                                                                                                    gatgagaGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTTAACAGTATTTGGGAAGGTGCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAATCATTTGTGGAAGTg                                 <  2:559114/147‑1 (MQ=255)
                                                                                                                    gaggagAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGtt                                 >  2:1649696/1‑149 (MQ=255)
                                                                                                                      ggagAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGaa                               >  1:12122/1‑150 (MQ=255)
                                                                                                                           gTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGc                          >  1:1162877/1‑150 (MQ=255)
                                                                                                                             gAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAg                        >  2:1291169/1‑150 (MQ=255)
                                                                                                                                cTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTTTTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGGATTTGTGGAAGTGATAGGGGAGGAt                     >  1:369889/1‑150 (MQ=255)
                                                                                                                                  gCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCg                   <  2:598135/150‑1 (MQ=255)
                                                                                                                                  gCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCg                   >  2:853228/1‑150 (MQ=255)
                                                                                                                                   ccGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGGTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGc                  <  1:535369/150‑1 (MQ=255)
                                                                                                                                       ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAATGGGCTGGAGGGCTagt              >  2:118544/4‑148 (MQ=255)
                                                                                                                                       ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt              >  2:236232/4‑150 (MQ=255)
                                                                                                                                       ctgtttACCGGCATCTCCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtt              >  1:458656/4‑150 (MQ=255)
                                                                                                                                         cTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTAtttt            >  2:444503/1‑150 (MQ=255)
                                                                                                                                            tACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCTATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGAt         >  2:1338723/1‑150 (MQ=255)
                                                                                                                                            tACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGAt         >  2:946121/1‑150 (MQ=255)
                                                                                                                                            tACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGAGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGAt         >  2:957772/1‑150 (MQ=255)
                                                                                                                                              ccGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATcg       <  2:8884/150‑1 (MQ=255)
                                                                                                                                                   agaTGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGAGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATAGCTATTTTGATCGCGCgg  >  1:18314/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCAAAGAAGCCCTGTGGGCGCGTTTATTTACGCCTGCGGCAAAATTTGTGATCTGCAAAAGAGGTATGCCGTTTGTGGCCGAAGCGATGGAGGTGCTGGGCGGCATTGGTTATTGCGAGGAGAGTGAGCTGCCGCGGCTTTACCGGGAGATGCCGGTAAACAGTATTTGGGAAGGTTCCGGCAATATTATGTGCCTGGATGTGTTGCGCGTTCTCAATAAGCAAGCGGGCGTATACGACTTATTGTCGGAAGCATTTGTGGAAGTGAAAGGGCAGGATCGCTATTTTGATCGCGCGG  >  NC_000913/4415377‑4415673

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: