Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,517 T→A 15.2% intergenic (+42/‑206) ygaM → / → nrdH putative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,5170TA15.2% 316.6 / 46.8 125intergenic (+42/‑206)ygaM/nrdHputative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (62/44);  new base A (9/10);  total (71/54)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

Only 100 of 128 total aligned reads displayed.

TACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGT  >  NC_000913/2800418‑2800616
                                                                                                   |                                                                                                   
tACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                                     >  1:1356718/1‑100 (MQ=255)
tACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                                     >  1:1426610/1‑100 (MQ=255)
 aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTg                                                                                                    <  1:1801819/100‑1 (MQ=255)
 aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCAAAGCCTg                                                                                                    <  1:2042028/100‑1 (MQ=255)
  cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGt                                                                                                   >  1:1386549/1‑100 (MQ=255)
    gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                 >  1:2733690/1‑100 (MQ=255)
    gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                 <  1:388468/100‑1 (MQ=255)
    gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                 >  1:476821/1‑100 (MQ=255)
    gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                 >  1:852012/1‑100 (MQ=255)
    gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                 >  2:591994/1‑100 (MQ=255)
    gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                 >  1:1603500/1‑100 (MQ=255)
      agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                               >  2:602139/1‑100 (MQ=255)
      agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                               >  1:1877040/1‑100 (MQ=255)
      agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAAAGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                               >  1:311341/1‑100 (MQ=255)
       gcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCACAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                              >  1:322455/1‑100 (MQ=255)
       gcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                              >  1:250943/1‑100 (MQ=255)
          gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgg                                                                                             <  1:1166082/98‑1 (MQ=255)
          gCGGGAGGGAATTTTATCGGCGCACTGTTAAGCATGCGCAAAAAAAAAAAAAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                           <  1:2138120/100‑1 (MQ=255)
           cGGTAGGTGTTCTTATCGGGGCACTGATAAGAATGCGGAAATCGTAGTGCAAAAATGATCATAAATACGGGTCTTTGACACAGAAGGCTGTCTTCggg                                                                                            <  1:1473671/98‑1 (MQ=255)
           cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                                          <  2:676857/100‑1 (MQ=255)
           cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                          <  1:2084233/100‑1 (MQ=255)
           cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                          <  2:1537645/100‑1 (MQ=255)
           cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                          <  1:431417/100‑1 (MQ=255)
           cGGTAGGTATTTTTATCGGCGCACGGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                          <  1:2441613/100‑1 (MQ=255)
            ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                         >  2:2255167/1‑100 (MQ=255)
            ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtc                                                                                         <  2:1403823/100‑2 (MQ=255)
            ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAAAAAAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtc                                                                                         <  1:2354764/100‑2 (MQ=255)
             gtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGttt                                                                                        >  2:900495/1‑100 (MQ=255)
               aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtcaaa                                                                                      <  2:670101/100‑5 (MQ=38)
                ggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                                     >  1:830838/1‑100 (MQ=255)
                ggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                                     >  2:2199528/1‑100 (MQ=255)
                ggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtcaaag                                                                                     <  2:1134794/100‑6 (MQ=38)
                  taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtttt                                                                                   >  1:1329730/1‑100 (MQ=255)
                    tttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                                 <  2:2739953/100‑1 (MQ=255)
                     ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTATTGCAAAAATGATAATAAATACGCGACTTTGACCCCGATGCCTGTCTTCGGAGTATCttttt                                                                                  >  1:1118102/1‑98 (MQ=255)
                        tATCGGCGCACTGTTAAGCATGCGCAAATCGTAGGGAAAAAATGATAATAGATACGCGGATTTGACGCCGAAGACTGTCTTCGGGGTTTCTTTGTGc                                                                                >  1:2852447/1‑97 (MQ=255)
                         aTCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                            >  1:27566/1‑100 (MQ=255)
                           cGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                          <  2:144224/100‑1 (MQ=255)
                              cgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAAt                                                                       >  2:2093224/1‑100 (MQ=255)
                               gcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTCGCCTGGTGAATc                                                                      <  1:523292/100‑1 (MQ=255)
                                 aCTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                    >  2:1592147/1‑100 (MQ=255)
                                 aCTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                                    <  1:2093217/100‑1 (MQ=255)
                                  cTGCTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcaca                                                                   >  1:2188991/1‑100 (MQ=255)
                                    gTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaa                                                                 <  1:472534/100‑1 (MQ=255)
                                     ttAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaa                                                                <  1:302574/100‑1 (MQ=255)
                                      tAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                               <  2:846539/100‑1 (MQ=255)
                                        aGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATc                                                             >  2:534608/1‑100 (MQ=255)
                                          cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                           >  1:2251004/1‑100 (MQ=255)
                                            tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccccc                                                         <  2:296467/100‑1 (MQ=255)
                                            tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccccc                                                         >  2:238747/1‑100 (MQ=255)
                                              cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccccc                                                         <  2:2813823/98‑1 (MQ=255)
                                              cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTa                                                       >  1:1260236/1‑100 (MQ=255)
                                               gcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                      >  1:2435063/1‑100 (MQ=255)
                                                      gTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                               >  2:2599849/1‑100 (MQ=255)
                                                         gTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGc                                            >  2:2139049/1‑100 (MQ=255)
                                                         gTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGc                                            <  2:823838/100‑1 (MQ=255)
                                                         gTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGc                                            >  2:2743831/1‑100 (MQ=255)
                                                         gTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGc                                            >  2:413641/1‑100 (MQ=255)
                                                          tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGGCTGTATTCAGGGTTTCTTTTTGACTGTTGGATCACAAAAATCCCCCTATTCCGTCagg                                             >  1:1912705/1‑96 (MQ=255)
                                                          tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                           >  1:2867521/1‑100 (MQ=255)
                                                           gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                          >  2:953494/1‑100 (MQ=255)
                                                           gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                          >  2:418523/1‑100 (MQ=255)
                                                            cAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                         >  2:2906324/1‑100 (MQ=255)
                                                            cAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                         <  2:2740445/100‑1 (MQ=255)
                                                            cAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                         <  2:673019/100‑1 (MQ=255)
                                                            cAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                         >  1:2189154/1‑100 (MQ=255)
                                                             aaaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCat                                        >  2:1457507/1‑100 (MQ=255)
                                                               aaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCatat                                      <  2:167998/100‑1 (MQ=255)
                                                                aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                     <  2:443540/100‑1 (MQ=255)
                                                                aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                     <  1:2193586/100‑1 (MQ=255)
                                                                 aTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATcc                                    <  2:2003662/100‑1 (MQ=255)
                                                                  tGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCa                                   >  1:2880257/1‑100 (MQ=255)
                                                                   gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                  <  2:1932520/100‑1 (MQ=255)
                                                                        taaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTaa                             >  1:1461659/1‑100 (MQ=255)
                                                                         aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAAt                            <  2:796708/100‑1 (MQ=255)
                                                                             aCGCGTTTTTGAACACGAAACCTGTATTCTGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCTCCTTCCCCGTCACGCTCATATCCAGGGTAATTTc                         <  1:1362001/99‑1 (MQ=255)
                                                                              cgcgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGa                       <  2:2172681/100‑1 (MQ=255)
                                                                               gcgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                      <  1:2033335/100‑1 (MQ=255)
                                                                                cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAcc                     <  2:875816/100‑1 (MQ=255)
                                                                                cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAcc                     <  1:1925842/100‑1 (MQ=255)
                                                                                       aaCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAttt              >  2:59551/2‑100 (MQ=255)
                                                                                        aCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg             >  2:1007345/1‑100 (MQ=255)
                                                                                        aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg             >  1:1787163/1‑100 (MQ=255)
                                                                                         tccgatctccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc            <  2:1069412/92‑1 (MQ=255)
                                                                                         ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc            <  2:291683/100‑1 (MQ=255)
                                                                                         ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc            <  1:101676/100‑1 (MQ=255)
                                                                                         ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc            >  2:171963/1‑100 (MQ=255)
                                                                                         ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc            >  1:2560822/1‑100 (MQ=255)
                                                                                          cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt           >  2:172/1‑100 (MQ=255)
                                                                                           ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCta          <  1:2414238/100‑1 (MQ=255)
                                                                                           ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCta          >  1:445125/1‑100 (MQ=255)
                                                                                           ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCta          >  2:793457/1‑100 (MQ=255)
                                                                                           ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCta          <  2:1603504/100‑1 (MQ=255)
                                                                                           ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAAACCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg             <  1:1993452/97‑1 (MQ=255)
                                                                                              aaGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtatat       <  2:1724808/100‑1 (MQ=255)
                                                                                                 ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGcagatcgga    >  1:1069410/1‑92 (MQ=255)
                                                                                                  cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtg   <  1:533972/100‑1 (MQ=255)
                                                                                                  cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtg   >  2:1987459/1‑100 (MQ=255)
                                                                                                  cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTCGCTatatat     >  1:1719768/1‑98 (MQ=255)
                                                                                                   tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgt  >  2:2268499/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGT  >  NC_000913/2800418‑2800616

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