Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 8.5% A97P (GCA→CCA)  yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5890GC 8.5% 397.1 / 21.8 118A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (49/58);  new base C (5/5);  total (55/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 125 total aligned reads displayed.

CTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCA  >  NC_000913/252441‑252737
                                                                                                                                                    |                                                                                                                                                    
cTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                                                                      >  4:491757/1‑149 (MQ=255)
   tACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCAt                                                                                                                                                   <  3:706241/149‑1 (MQ=255)
    aCCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCAt                                                                                                                                                   <  1:356776/148‑1 (MQ=255)
      cTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                                                                                                                                                 >  1:119914/1‑148 (MQ=255)
       tGGCCATTAATCACATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                <  3:218748/148‑1 (MQ=255)
       tGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGTGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                >  4:218748/1‑148 (MQ=255)
        ggCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                                                              >  4:271815/1‑149 (MQ=255)
            aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                                                                           >  1:489382/1‑148 (MQ=255)
            aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                          <  3:770708/149‑1 (MQ=255)
            aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGGTTGATGAGCAAGCATGGTTGCTc                                                                                                                                           >  1:161634/1‑148 (MQ=255)
             ttAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                          >  2:211883/1‑148 (MQ=255)
             ttAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                          >  1:482484/1‑148 (MQ=255)
              tAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                                                                         <  1:169042/148‑1 (MQ=255)
                 tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGGGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                     >  4:548835/1‑149 (MQ=255)
                  cAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                     >  2:390869/1‑148 (MQ=255)
                  cAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                     >  1:140751/1‑148 (MQ=255)
                  cAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCa                                                                                                                                    >  3:367663/1‑149 (MQ=255)
                    aaTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAtt                                                                                                                                  <  4:367663/149‑1 (MQ=255)
                     aTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAtt                                                                                                                                  >  1:276288/1‑148 (MQ=255)
                      tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTc                                                                                                                                           <  3:396842/138‑1 (MQ=255)
                      tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTc                                                                                                                                           >  4:396842/1‑138 (MQ=255)
                             gAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCt                                                                                                                         <  3:902480/149‑1 (MQ=255)
                                ttGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTc                                                                                                                                           <  4:402603/128‑1 (MQ=255)
                                   cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                                                                   <  3:491757/149‑1 (MQ=255)
                                     cATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                  >  2:52039/1‑148 (MQ=255)
                                     cATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                                                                 <  4:75280/149‑1 (MQ=255)
                                       tatTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCAt                                                                                                                >  1:102533/1‑148 (MQ=255)
                                         ttCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaa                                                                                                              <  1:211883/148‑1 (MQ=255)
                                          tCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAAc                                                                                                            <  3:413419/149‑1 (MQ=255)
                                          tCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGCGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                                                                             >  2:518666/1‑148 (MQ=255)
                                           cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                                                                           >  3:313065/1‑149 (MQ=255)
                                            aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                          >  3:684546/1‑149 (MQ=255)
                                             tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                          >  1:143486/1‑148 (MQ=255)
                                             tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                          >  1:497330/1‑148 (MQ=255)
                                               tGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                        <  2:272575/148‑1 (MQ=255)
                                                       gagaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                                                                               <  4:192125/149‑1 (MQ=255)
                                                       gagaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATa                                                                                                >  1:525120/1‑148 (MQ=255)
                                                        agaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaa                                                                                               <  2:419970/148‑1 (MQ=255)
                                                           ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacac                                                                                            <  2:818351/148‑2 (MQ=255)
                                                           ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaa                                                                                            <  2:119914/148‑1 (MQ=255)
                                                             aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                                                                          >  1:86608/1‑148 (MQ=255)
                                                               tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                                                                       <  4:224097/149‑1 (MQ=255)
                                                               tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                                                                       >  3:645563/1‑149 (MQ=255)
                                                                  aTTTCCGCCACAGTTTCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTTGTATATTGTCAGGGGGTGTTTGATGAGCAATCATGGTTGCGCATTGCCATTCTGAAACCTAAACCTCATAAACTGGCTCGCTATaacaacagatcgt                                                                                    >  4:753007/1‑142 (MQ=255)
                                                                      ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                                                                                >  4:276222/1‑149 (MQ=255)
                                                                        tccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                                                                           >  2:771699/2‑88 (MQ=255)
                                                                        gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGCTGAGCAAGCATGGTTGCTc                                                                                                                                           <  1:771699/88‑1 (MQ=255)
                                                                         ccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaaa                                                                              >  1:947900/1‑148 (MQ=255)
                                                                          cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAt                                                                             <  2:544820/148‑1 (MQ=255)
                                                                            caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                                                                          <  4:84227/149‑1 (MQ=255)
                                                                            caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                           <  2:187012/148‑1 (MQ=255)
                                                                            caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                           <  1:670418/148‑1 (MQ=255)
                                                                                      aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                 >  2:951552/1‑148 (MQ=255)
                                                                                      aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa                                                                <  3:548835/149‑1 (MQ=255)
                                                                                            gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg                                                          >  4:22189/1‑149 (MQ=255)
                                                                                             cagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg                                                          >  1:390994/1‑148 (MQ=255)
                                                                                                aaCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt                                                       <  2:497330/148‑1 (MQ=255)
                                                                                                       gACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGc                                                >  1:881386/1‑148 (MQ=255)
                                                                                                        aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGt                                              >  3:404873/1‑149 (MQ=255)
                                                                                                           aaGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAGCTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGttt                                            >  2:274417/1‑148 (MQ=255)
                                                                                                           aaGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGttt                                            <  2:159198/148‑1 (MQ=255)
                                                                                                            aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtttt                                           <  1:787417/148‑1 (MQ=255)
                                                                                                            aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtttt                                           <  2:62593/148‑1 (MQ=255)
                                                                                                              cgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                 >  2:720117/1‑124 (MQ=255)
                                                                                                               gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGaa                                       <  4:37596/149‑1 (MQ=255)
                                                                                                               gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGaa                                       <  3:859080/149‑1 (MQ=255)
                                                                                                                cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtttt                                           >  3:503823/1‑144 (MQ=255)
                                                                                                                cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtttt                                           <  4:503823/144‑1 (MQ=255)
                                                                                                                cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAt                                      <  3:98496/149‑1 (MQ=255)
                                                                                                                     ggTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGACATTCTGTAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTtat                                  <  2:390994/148‑1 (MQ=255)
                                                                                                                      gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTatat                                >  3:619440/1‑149 (MQ=255)
                                                                                                                        atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt                                                       <  2:800992/124‑1 (MQ=255)
                                                                                                                        atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGttt                                                       >  1:800992/1‑124 (MQ=255)
                                                                                                                        atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatt                               <  2:525120/148‑1 (MQ=255)
                                                                                                                         tatTGTCTGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATATACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatta                              <  1:747870/148‑1 (MQ=255)
                                                                                                                         tatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatta                              <  2:160086/148‑1 (MQ=255)
                                                                                                                            tGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGa                           <  2:947900/148‑1 (MQ=255)
                                                                                                                             gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGaa                          <  1:565976/148‑1 (MQ=255)
                                                                                                                                 gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAAc                      <  1:167217/148‑1 (MQ=255)
                                                                                                                                 gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAAc                      <  1:871667/148‑1 (MQ=255)
                                                                                                                                  ggggCATTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAATGGCAGTAGCGTTTCGCATGCGATTTTGAATTTATATTATGAATAACATa                   <  4:304822/149‑1 (MQ=255)
                                                                                                                                   gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATa                   >  4:892045/1‑149 (MQ=255)
                                                                                                                                   gggCGTTTGATGAGCAAGCATCGTTGCACATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATCCATAAATTTGGCAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATAATATGAATAACATa                   <  4:750137/149‑1 (MQ=255)
                                                                                                                                     gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACa                 >  3:50056/1‑149 (MQ=255)
                                                                                                                                      cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACa                 >  1:875771/1‑148 (MQ=255)
                                                                                                                                       gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa               >  4:243138/1‑149 (MQ=255)
                                                                                                                                        tttGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa               >  1:416164/1‑148 (MQ=255)
                                                                                                                                           gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg           <  3:343297/149‑1 (MQ=255)
                                                                                                                                           aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc                                                      >  3:427476/2‑106 (MQ=255)
                                                                                                                                           aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc                                                      <  4:427476/105‑1 (MQ=255)
                                                                                                                                            atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGa          <  4:819144/149‑1 (MQ=255)
                                                                                                                                            atgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGa                           <  1:653862/132‑1 (MQ=255)
                                                                                                                                            atgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGa                           >  2:653862/1‑132 (MQ=255)
                                                                                                                                            atgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg           >  1:657223/1‑148 (MQ=255)
                                                                                                                                              gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGaaa        <  4:92611/149‑1 (MQ=255)
                                                                                                                                              gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGaa         >  2:696145/1‑148 (MQ=255)
                                                                                                                                               agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAAc       >  4:733757/1‑149 (MQ=255)
                                                                                                                                               agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTCGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGaaa        >  2:288195/1‑148 (MQ=255)
                                                                                                                                                  aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAt                                      >  3:170524/1‑115 (MQ=255)
                                                                                                                                                    gcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCa  >  3:416903/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCA  >  NC_000913/252441‑252737

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