Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,910,230 G→C 6.2% C152W (TGC→TGG mazG ← nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,910,2300GC 6.2% 382.2 / 11.4 113C152W (TGC→TGGmazGnucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (62/44);  new base C (5/2);  total (67/46)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 115 total aligned reads displayed.

GGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAAC  >  NC_000913/2910086‑2910378
                                                                                                                                                |                                                                                                                                                    
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)
         cagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGCGCCAACGttggcttcg                                                                                                                                        <  1:463046/150‑8 (MQ=255)
           ggTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGAttt                                                                                                                                      <  1:1030196/150‑1 (MQ=255)
               ccccATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTg                                                                                                                                  >  2:808451/1‑150 (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:2822/150‑1 (MQ=255)
                                    aGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTaa                                                                                                             >  1:739475/1‑150 (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)
                                                     ccTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGt                                                                                            >  1:1227070/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)
                                                              cTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCg                                                                                   >  1:947395/1‑150 (MQ=255)
                                                              cTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCg                                                                                   <  1:48776/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)
                                                                                     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)
                                                                                           tAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCg                                                      >  1:1495707/1‑150 (MQ=255)
                                                                                           tAGACTTTATCGACTACCGGACCAAGCGGCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCg                                                      <  2:1379532/150‑1 (MQ=255)
                                                                                               cTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGttt                                                  <  2:457696/150‑1 (MQ=255)
                                                                                                  tATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGa                                               <  2:437944/150‑1 (MQ=255)
                                                                                                  tATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTATGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGa                                               >  2:1077003/1‑150 (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:1075358/1‑150 (MQ=255)
                                                                                                          aCCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTcc                                       >  1:1080702/1‑150 (MQ=255)
                                                                                                           ccGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTccc                                      >  1:418602/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)
                                                                                                                         gtcCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACt                        >  2:144276/1‑150 (MQ=255)
                                                                                                                           ccAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTc                      >  2:412855/1‑150 (MQ=255)
                                                                                                                            cAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACt                        <  1:1361978/147‑1 (MQ=255)
                                                                                                                            cAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACt                        >  2:1361978/1‑147 (MQ=255)
                                                                                                                             aaTCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCac                    >  1:653052/1‑150 (MQ=255)
                                                                                                                             aaTCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCac                    >  1:1492551/1‑150 (MQ=255)
                                                                                                                               tCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGc                           >  1:1148927/1‑141 (MQ=255)
                                                                                                                               tCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGc                           <  2:1148927/141‑1 (MQ=255)
                                                                                                                               tCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCacta                  >  1:1370998/1‑150 (MQ=255)
                                                                                                                                 gAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCactact                >  1:1090601/1‑150 (MQ=255)
                                                                                                                                 gAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCactact                >  1:1192366/1‑150 (MQ=255)
                                                                                                                                   cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt              >  1:948474/5‑150 (MQ=255)
                                                                                                                                   cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt              >  2:568558/5‑150 (MQ=255)
                                                                                                                                   cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt              >  2:1534819/5‑150 (MQ=255)
                                                                                                                                   cagaaaCGTTGCGCCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt              >  1:846929/5‑150 (MQ=255)
                                                                                                                                   aGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAg                            >  1:1414970/1‑136 (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        <  1:461772/150‑1 (MQ=255)
                                                                                                                                         cGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCgg        <  2:583496/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)
                                                                                                                                               cgcAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAAc  <  2:857408/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGCAAACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCGAAGCCAACGTTGGCGCAACGTTTCTGGATTTTTTGCGCACGCATTAAAGCCGGTAAACTACGAGGAATATCGTCCAGCGCCGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAAC  >  NC_000913/2910086‑2910378

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: