Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,910,228 G→C 6.2% A153G (GCC→GGC) ‡ mazG ← nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,910,2280GC 6.2% 382.1 / 12.1 112A153G (GCC→GGC) ‡mazGnucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (62/43);  new base C (5/2);  total (67/45)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 116 total aligned reads displayed.

TAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAA  >  NC_000913/2910079‑2910377
                                                                                                                                                     |                                                                                                                                                     
tAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTgg                                                                                                                                                       >  2:1077457/1‑150 (MQ=255)
 aaCCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGgc                                                                                                                                                      >  1:1279291/1‑150 (MQ=255)
       ggCAAACATCAGGACCACCATTTCCACCTCAAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGCCCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACg                                                                                                                                                >  2:1414609/1‑150 (MQ=255)
         cAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGtt                                                                                                                                              >  2:917466/1‑150 (MQ=255)
            aCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCt                                                                                                                                           >  2:419694/1‑150 (MQ=255)
            aCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCt                                                                                                                                           >  2:284299/1‑150 (MQ=255)
             cagcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTg                                                                                                                                          <  1:1224888/150‑1 (MQ=255)
             aagcagGTCCCCAATTTCCTCCTCAAGTTTAGCCTGGTCGACAAAAGCCTGCCGAGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTg                                                                                                                                          <  2:191087/149‑1 (MQ=255)
              agcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTGTCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTgg                                                                                                                                         >  1:1424242/1‑150 (MQ=255)
               gcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGCGCCAACGttggcttc                                                                                                                                        <  1:1442090/150‑7 (MQ=255)
                cagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGCGCCAACGttggcttcg                                                                                                                                       <  1:463046/150‑8 (MQ=255)
                  ggTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGAttt                                                                                                                                     <  1:1030196/150‑1 (MQ=255)
                                 cctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGAtttttt                                                                                                                                  >  2:253613/1‑138 (MQ=255)
                                 cctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGAtttttt                                                                                                                                  <  1:253613/138‑1 (MQ=255)
                                 cctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTa                                                                                                                      >  1:1097088/1‑150 (MQ=255)
                                   tccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTaaa                                                                                                                    >  1:634811/1‑150 (MQ=255)
                                           aGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTaa                                                                                                            >  1:739475/1‑150 (MQ=255)
                                           aGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTaa                                                                                                            <  1:2822/150‑1 (MQ=255)
                                            gCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTGTTCGTAGACTTTATCGACTACCGGCCCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTATGGATTTTTTGCGCACGCATTAAAGCCGGTaaa                                                                                                           <  2:814050/150‑1 (MQ=255)
                                            gCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTaaa                                                                                                           >  2:469152/1‑150 (MQ=255)
                                              cTGGTCGACAACAGCCTGCCGCGCTTCGTTCATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACt                                                                                                         >  1:1147646/1‑150 (MQ=255)
                                              cTGGTCGACAACAGCCTGCCGCGCTTCGTACATCCCCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACt                                                                                                         >  2:1110316/1‑150 (MQ=255)
                                              cTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACt                                                                                                         >  1:857408/1‑150 (MQ=255)
                                                ggTCAACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTAc                                                                                                       >  2:112680/1‑150 (MQ=255)
                                                 gTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACg                                                                                                      <  2:522463/150‑1 (MQ=255)
                                                    gACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAgg                                                                                                   <  1:160760/150‑1 (MQ=255)
                                                     acaacaGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGa                                                                                                  >  1:1006953/1‑150 (MQ=255)
                                                     acaacaGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGa                                                                                                  >  2:159924/1‑150 (MQ=255)
                                                      caacaGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGaa                                                                                                 >  2:1282382/1‑150 (MQ=255)
                                                      caacaGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGaa                                                                                                 >  1:437944/1‑150 (MQ=255)
                                                       aacaGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAat                                                                                                <  2:86760/150‑1 (MQ=255)
                                                        acaGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAata                                                                                               >  2:1174548/1‑150 (MQ=255)
                                                          aGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATc                                                                                             >  2:1255988/1‑150 (MQ=255)
                                                           gCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCg                                                                                            >  1:446156/1‑150 (MQ=255)
                                                              tGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTcc                                                                                         <  1:558740/150‑1 (MQ=255)
                                                                ccGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAg                                                                                       <  1:284299/150‑1 (MQ=255)
                                                                   cgcTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAgcgc                                                                                    <  1:317479/150‑1 (MQ=255)
                                                                     cTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGGCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCg                                                                                  <  1:1119495/150‑1 (MQ=255)
                                                                       tCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGaa                                                                                >  2:461772/1‑150 (MQ=255)
                                                                       tCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGaa                                                                                >  2:464065/1‑150 (MQ=255)
                                                                         gTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAtg                                                                              >  2:514655/1‑150 (MQ=255)
                                                                          tACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAtgc                                                                             >  1:1132677/1‑150 (MQ=255)
                                                                           aCATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAtgct                                                                            >  1:534915/1‑150 (MQ=255)
                                                                             aTCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAtgctgc                                                                          >  1:1443541/1‑150 (MQ=255)
                                                                             aTCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAAtgctgc                                                                          >  1:1448955/1‑150 (MQ=255)
                                                                               cACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTgcgc                                                                        <  1:419694/150‑1 (MQ=255)
                                                                               cACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTgcgc                                                                        <  2:897069/150‑1 (MQ=255)
                                                                                 ccTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCtt                                                                      <  2:1006953/150‑1 (MQ=255)
                                                                                   tcgtcgATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTc                                                                    <  2:601022/150‑1 (MQ=255)
                                                                                   tcgtcgATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTc                                                                    <  2:687598/150‑1 (MQ=255)
                                                                                    cgtcgATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCt                                                                   >  2:1075150/1‑150 (MQ=255)
                                                                                    cgtcgATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCt                                                                   >  2:1038325/1‑150 (MQ=255)
                                                                                           ctctTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTgcgcgcg                                                            <  2:1132677/150‑1 (MQ=255)
                                                                                            tctTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTgcgcgcgc                                                           >  2:899385/1‑150 (MQ=255)
                                                                                               tCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGctct                                                        >  1:112723/1‑150 (MQ=255)
                                                                                                cGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCtt                                                       >  2:1019064/1‑150 (MQ=255)
                                                                                                 gTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTc                                                      >  2:485000/1‑150 (MQ=255)
                                                                                                  tAGACTTTATCGACTACCGGACCAAGCGGCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCg                                                     <  2:1379532/150‑1 (MQ=255)
                                                                                                   aGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCgg                                                    <  2:1313202/150‑1 (MQ=255)
                                                                                                         tATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGa                                              <  1:78462/150‑1 (MQ=255)
                                                                                                         tATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGa                                              <  2:437944/150‑1 (MQ=255)
                                                                                                          aTCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGAt                                             <  1:176894/150‑1 (MQ=255)
                                                                                                            cGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGAttt                                           >  1:1411109/1‑150 (MQ=255)
                                                                                                            cGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGAttt                                           <  2:115825/150‑1 (MQ=255)
                                                                                                             gACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTg                                          >  2:1052522/1‑150 (MQ=255)
                                                                                                             gACTACCGGACCAAGCGGCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTg                                          <  1:1506470/150‑1 (MQ=255)
                                                                                                               cTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCt                                        >  1:606558/1‑150 (MQ=255)
                                                                                                                 aCCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTcc                                      <  2:1158341/150‑1 (MQ=255)
                                                                                                                 aCCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTcc                                      >  1:1080702/1‑150 (MQ=255)
                                                                                                                 aCCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTcc                                      >  1:1075358/1‑150 (MQ=255)
                                                                                                                       ccAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACgg                                >  1:1303421/1‑150 (MQ=255)
                                                                                                                        cAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGAGCGCTCTTCGGTTTTGATTTGCTCCCAACggg                               >  1:687011/1‑150 (MQ=255)
                                                                                                                            cgtcgtCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAg                           <  2:1028824/150‑1 (MQ=255)
                                                                                                                            cgtcgtCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAg                           <  2:1147646/150‑1 (MQ=255)
                                                                                                                             gtcgtcCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGc                          >  1:582278/1‑150 (MQ=255)
                                                                                                                              tcgtcCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCa                         <  1:125705/150‑1 (MQ=255)
                                                                                                                                gtcCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACt                       >  2:144276/1‑150 (MQ=255)
                                                                                                                                  ccAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTc                     >  2:412855/1‑150 (MQ=255)
                                                                                                                                   cAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACt                       <  1:1361978/147‑1 (MQ=255)
                                                                                                                                   cAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCa                    <  2:446156/150‑1 (MQ=255)
                                                                                                                                    aaTCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCac                   >  1:1492551/1‑150 (MQ=255)
                                                                                                                                      tCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGc                          >  1:1148927/1‑141 (MQ=255)
                                                                                                                                      tCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCacta                 >  1:1370998/1‑150 (MQ=255)
                                                                                                                                        gAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCactact               >  1:1090601/1‑150 (MQ=255)
                                                                                                                                        gAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCactact               >  1:1192366/1‑150 (MQ=255)
                                                                                                                                          cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt             >  2:1534819/5‑150 (MQ=255)
                                                                                                                                          cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt             >  1:948474/5‑150 (MQ=255)
                                                                                                                                          cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt             >  2:568558/5‑150 (MQ=255)
                                                                                                                                          cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt             >  1:846929/5‑150 (MQ=255)
                                                                                                                                          aGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAg                           <  2:1414970/136‑1 (MQ=255)
                                                                                                                                          aGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt             >  2:898145/1‑150 (MQ=255)
                                                                                                                                           agaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACAAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGtt            >  1:62680/4‑150 (MQ=255)
                                                                                                                                             cAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGtttt          >  2:690528/1‑150 (MQ=255)
                                                                                                                                              aaCGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTc         <  1:485000/150‑1 (MQ=255)
                                                                                                                                                cGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCgg       <  2:583496/150‑1 (MQ=255)
                                                                                                                                                cGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCgg       <  1:461772/150‑1 (MQ=255)
                                                                                                                                                 gTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGc      <  2:739475/150‑1 (MQ=255)
                                                                                                                                                 gTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGc      <  1:721930/150‑1 (MQ=255)
                                                                                                                                                   tGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAg    <  2:1103552/150‑1 (MQ=255)
                                                                                                                                                     gcgcAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGaa  >  1:1493673/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAACCGTGGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAA  >  NC_000913/2910079‑2910377

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