Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 280,618 T→C 5.8% pseudogene (118/309 nt) yagB ← pseudogene, CP4‑6 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913280,6180TC 5.8% 202.3 / 3.3 71pseudogene (118/309 nt)yagBpseudogene, CP4‑6 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (34/31);  new base C (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

CACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTT  >  NC_000913/280471‑280765
                                                                                                                                                   |                                                                                                                                                   
cACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGt                                                                                                                                                     >  3:302454/1‑148 (MQ=255)
 aCCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTg                                                                                                                                                    <  4:251960/148‑1 (MQ=255)
 aCCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTg                                                                                                                                                    <  1:319447/148‑1 (MQ=255)
  ccAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTg                                                                                                                                                    >  1:142261/1‑147 (MQ=255)
    aGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAAt                                                                                                                                                 <  7:115241/148‑1 (MQ=255)
    aGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAAt                                                                                                                                                 >  8:115241/1‑148 (MQ=255)
           ccGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGGGAATCCAGccc                                                                                                                                          <  6:448785/148‑1 (MQ=255)
              ttACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGt                                                                                                                                       <  4:3158/148‑1 (MQ=255)
               tACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGt                                                                                                                                       >  3:179089/1‑147 (MQ=255)
               tACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGt                                                                                                                                       >  3:172606/1‑147 (MQ=255)
                 cGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGGCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCa                                                                                                                                     <  5:1301/147‑1 (MQ=255)
                   gTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGc                                                                                                                                   >  1:179034/1‑147 (MQ=255)
                      cAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGa                                                                                                                               <  6:293861/148‑1 (MQ=255)
                          gtgACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAAt                                                                                                                           >  7:208790/1‑148 (MQ=255)
                              cacaGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATg                                                                                                                          >  5:242511/1‑145 (MQ=255)
                              cacaGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAg                                                                                                                       <  7:263190/148‑1 (MQ=255)
                               acaGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGCGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGa                                                                                                                      >  4:421570/1‑148 (MQ=255)
                                   tGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGAt                                                                                                                  <  2:263447/148‑1 (MQ=255)
                                       cgTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGAcc                                                                                                              <  8:208790/148‑1 (MQ=255)
                                            cTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCt                                                                                                         >  1:180374/1‑148 (MQ=255)
                                             tGGTATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATTCGGCAAGGATTCATCCAGTTGTAGAGTCTGCACCTCCATGAATGAGCCGGATAATCCAGCCCGGTCAGCCTGAAAATGCAGACGATTACATTCCTg                                                                                                        >  5:263850/1‑148 (MQ=255)
                                             tGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCCCCTCACTGAATGAGCCGGTGAATCCAGCCCGGACAGCCAGAAAATGCAGACGATGACCTTCCTg                                                                                                        >  5:215469/1‑148 (MQ=255)
                                               ggATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCa                                                                                                      >  1:96577/1‑148 (MQ=255)
                                                  ttCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCa                                                                                                   <  2:142261/148‑1 (MQ=255)
                                                   tCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGAAGACCTTCCTGCACCAg                                                                                                  >  5:59191/1‑148 (MQ=255)
                                                      gCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCg                                                                                               >  3:432973/1‑148 (MQ=255)
                                                          acaGGAGAGCAGAATCAGCTCCAGATCTGCGACAAAATGCGGGAAGGCTTCATCCAGTAGTAGAGTCTGCACCTCACTCAATGAGCCCGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATCACCTTCCTGCACCAGTCGTGCt                                                                                           >  1:450539/1‑148 (MQ=255)
                                                            aGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAg                                                                                                                                    >  6:399700/1‑105 (MQ=255)
                                                            aGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAg                                                                                                                                    <  5:399700/105‑1 (MQ=255)
                                                             ggAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCAGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTcc                                                                                         >  4:376076/1‑147 (MQ=255)
                                                              gagaGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGaaa                                                                                                                             <  8:197284/110‑1 (MQ=255)
                                                              gagaGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGaaa                                                                                                                             >  7:197284/1‑110 (MQ=255)
                                                              gagaGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGGCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGa                                                                                       <  2:157596/148‑1 (MQ=255)
                                                                     gcaTCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACgg                                                                                <  4:412481/148‑1 (MQ=255)
                                                                         cAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGACACACGgtgt                                                                             >  5:115910/1‑147 (MQ=255)
                                                                         cAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGgtgt                                                                             >  4:249093/1‑147 (MQ=255)
                                                                          aGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGtgtg                                                                            <  2:96577/147‑1 (MQ=255)
                                                                             tCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGAACTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGatat                                                                        >  8:381402/1‑148 (MQ=255)
                                                                                    gtgCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGt                                                                                                                                       >  1:179521/1‑78 (MQ=255)
                                                                                    gtgCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGt                                                                                                                                       <  2:179521/78‑1 (MQ=255)
                                                                                     tgCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGttt                                                                 <  3:249093/147‑1 (MQ=255)
                                                                                     tgCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGttt                                                                 <  6:242511/147‑1 (MQ=255)
                                                                                       cGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGAAGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAg                                                              <  6:59191/148‑1 (MQ=255)
                                                                                       cGACAAAATGCAGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTaa                                                               <  3:376076/147‑1 (MQ=255)
                                                                                         aCAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGCGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGgtgt                                                                             >  5:209766/1‑131 (MQ=255)
                                                                                         aCAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGCGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGgtgt                                                                             <  6:209766/131‑1 (MQ=255)
                                                                                              aTGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCaa                                                       <  4:172606/148‑1 (MQ=255)
                                                                                              aTGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCaa                                                       <  4:179089/148‑1 (MQ=255)
                                                                                                 cGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGt                                                     <  2:180374/147‑1 (MQ=255)
                                                                                                   ggAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGAGATCTCGTTTAAGCCGCCAAGttt                                                   >  6:235602/1‑147 (MQ=255)
                                                                                                     aaGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAg                                                <  2:179034/148‑1 (MQ=255)
                                                                                                      aGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAg                                                >  6:183415/1‑147 (MQ=255)
                                                                                                      aGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGt                                               <  5:294712/148‑1 (MQ=255)
                                                                                                          ttCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGa                                                                           <  4:74248/116‑1 (MQ=255)
                                                                                                          ttCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGa                                                                           >  3:74248/1‑116 (MQ=255)
                                                                                                           tCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGt                                          >  1:150027/1‑148 (MQ=255)
                                                                                                              tCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCtt                                       >  7:88836/1‑148 (MQ=255)
                                                                                                                cAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCttt                                      <  8:81251/147‑1 (MQ=255)
                                                                                                                cAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGGCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTa                                     <  6:215469/148‑1 (MQ=255)
                                                                                                                    tGTAGAGTCTGCACCTCACTGAATGCGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGCCGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCa                                 >  2:24808/1‑148 (MQ=255)
                                                                                                                         agTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGCCGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAAtt                            >  4:232105/1‑148 (MQ=255)
                                                                                                                            cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTg                                                  >  6:227610/1‑123 (MQ=255)
                                                                                                                            cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTg                                                  <  5:227610/123‑1 (MQ=255)
                                                                                                                            cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTa                         >  4:20140/1‑148 (MQ=255)
                                                                                                                            cTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTa                         >  8:201533/1‑148 (MQ=255)
                                                                                                                               cACCTCACTGAATGAGCCGGTGAATCCAGCCCGGGCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTAt                        <  5:351088/146‑1 (MQ=255)
                                                                                                                                 ccTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCt                    >  8:314405/1‑148 (MQ=255)
                                                                                                                                  cTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTg                   <  2:150027/148‑1 (MQ=255)
                                                                                                                                      cTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTa                >  6:43694/1‑147 (MQ=255)
                                                                                                                                             aGCCGGCGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTAtt                       >  2:214328/1‑133 (MQ=255)
                                                                                                                                             aGCCGGCGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTAtt                       <  1:214328/133‑1 (MQ=255)
                                                                                                                                                cGGTGAATCCAGCCCGGTAACCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCAt     >  8:85305/1‑148 (MQ=255)
                                                                                                                                                   tGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAAtttt                          <  4:110831/124‑1 (MQ=255)
                                                                                                                                                   tGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAAtttt                          >  3:110831/1‑124 (MQ=255)
                                                                                                                                                   tGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGtt  >  4:413688/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACCAGTCAGCCCGTTACGGTACAGTGTGACACAGTGCGCGTACCTGGGATTCAGCTCACAGGAGAGCAGCATCAGCTCCAGATGTGCGACAAAATGCGGGAAGGCTTCATCCAGTTGTAGAGTCTGCACCTCACTGAATGAGCCGGTGAATCCAGCCCGGTCAGCCAGAAAATGCAGACGATGACCTTCCTGCACCAGTCGTGCTCCGAAACACGGTGTGATATCTCGTTTAAGCCCCCAAGTTTGAGTAGAGTCTTTATTCATAATTTTATTCCTGTTAAAGAAGACCATGTT  >  NC_000913/280471‑280765

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