Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,159,287 C→A Y473* (TAC→TAA ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,2870CA98.8% 281.4 / NA 83Y473* (TAC→TAAptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (44/37);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGG  >  NC_000913/1159157‑1159432
                                                                                                                                  |                                                                                                                                                 
cTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATcc                                                                                                                                 <  4:11240/149‑1 (MQ=255)
    tgtgTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTaa                                                                                                                             <  1:119003/149‑1 (MQ=255)
     gtgtCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGt                                                                                                                                             <  5:39082/132‑1 (MQ=255)
     gtgtCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGt                                                                                                                                             >  6:39082/1‑132 (MQ=255)
     gtgtCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGAAGCGGGCGTAGTGGTGGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGAGGAGTAAATCCGTAACCACTAATCCGTaa                                                                                                                             <  6:356728/148‑1 (MQ=255)
          tAATGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGtt                                                                                                                       >  4:5460/1‑149 (MQ=255)
          tAACGTGGATCAGGCCGGCCGGAAGAAACGGGCCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGtt                                                                                                                       <  6:374037/149‑1 (MQ=255)
            aaGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCAc                                                                                                                                       <  1:31382/131‑1 (MQ=255)
            aaGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCAc                                                                                                                                       >  2:31382/1‑131 (MQ=255)
            aaGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTgg                                                                                                                     >  5:214623/1‑149 (MQ=255)
             aGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTggg                                                                                                                    <  5:132322/149‑1 (MQ=255)
              gTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTGCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTggg                                                                                                                    >  2:152800/1‑148 (MQ=255)
              gTGGATCAGGCCGGCCTGAAGAAACTGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTTGGga                                                                                                                  <  7:155896/149‑1 (MQ=255)
                 gATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCGGTAACCACTAATCCGTAAGACGTTGGGgaga                                                                                                                <  2:215524/149‑1 (MQ=255)
                  aTCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGAc                                                                                                               >  1:75072/1‑149 (MQ=255)
                     aGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTaa                                                                                                            <  3:151044/149‑1 (MQ=255)
                      ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAg                                                                                                           <  3:25554/149‑1 (MQ=255)
                      ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAg                                                                                                           >  3:188196/1‑149 (MQ=255)
                       gccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTa                                                                                                                                            <  6:108148/115‑1 (MQ=255)
                       gccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTa                                                                                                                                            >  5:108148/1‑115 (MQ=255)
                       gccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAgg                                                                                                          >  2:227965/1‑149 (MQ=255)
                           gccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGc                                                                                                      <  2:305148/149‑1 (MQ=255)
                           gccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGc                                                                                                      >  6:218699/1‑149 (MQ=255)
                             cTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGc                                                                                                      >  7:198043/1‑147 (MQ=255)
                              tGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGc                                                                                                      >  5:39558/1‑146 (MQ=255)
                                aagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                 <  6:39558/149‑1 (MQ=255)
                                aagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                 >  1:270866/1‑149 (MQ=255)
                                    aaCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCt                                                                                             <  2:65803/149‑1 (MQ=255)
                                    aaCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCt                                                                                             <  8:198043/149‑1 (MQ=255)
                                         ggCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtt                                                                                        >  4:102960/1‑149 (MQ=255)
                                               gCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                  <  1:200653/149‑1 (MQ=255)
                                                 gggCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCa                                                                                                        <  7:91693/125‑1 (MQ=255)
                                                 gggCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCa                                                                                                        >  8:91693/1‑125 (MQ=255)
                                                 gggCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                >  1:335636/1‑149 (MQ=255)
                                                 gggCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                <  1:227965/149‑1 (MQ=255)
                                                   gCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTATATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgt                                                                              <  1:152800/149‑1 (MQ=255)
                                                      tAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCTTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                 <  5:5935/143‑1 (MQ=255)
                                                      tAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCTTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAAGGCTGCCTTTTTTACa                                                                                 >  6:5935/1‑143 (MQ=255)
                                                        gTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTa                                                                          <  2:401755/148‑1 (MQ=255)
                                                         tGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTaat                                                                        <  5:317419/149‑3 (MQ=255)
                                                             tGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAg                                                                                                   >  2:332747/1‑118 (MQ=255)
                                                             tGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAg                                                                                                   <  1:332747/118‑1 (MQ=255)
                                                               cTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGaa                                                                   <  3:5460/148‑1 (MQ=255)
                                                               cTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                  <  3:102960/149‑1 (MQ=255)
                                                                tGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAATCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAtt                                                                 >  5:191173/1‑149 (MQ=255)
                                                                 ggTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACGGTGTTATTCAGAATTGa                                                               >  6:131522/1‑149 (MQ=255)
                                                                   ttCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAGGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                              <  4:240504/149‑1 (MQ=255)
                                                                   ttCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                              <  5:218699/149‑1 (MQ=255)
                                                                     cTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATAc                                                            >  7:197812/1‑149 (MQ=255)
                                                                        gtgtTCAGGCGATTTTCGGTACTAAATCCGATAATCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATAc                                                            <  6:191173/146‑1 (MQ=255)
                                                                        gtgtTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTg                                                         >  2:303601/1‑149 (MQ=255)
                                                                           ttCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCg                                                      <  6:214623/149‑1 (MQ=255)
                                                                             cAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                     <  1:303601/148‑1 (MQ=255)
                                                                             cAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                     >  5:303181/1‑148 (MQ=255)
                                                                              aGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTg                                                                                                                      >  2:347351/1‑82 (MQ=255)
                                                                              aGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTg                                                                                                                      <  1:347351/82‑1 (MQ=255)
                                                                              aGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                    >  3:147776/1‑148 (MQ=255)
                                                                               ggCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTaa                                                  >  2:130190/1‑149 (MQ=255)
                                                                                         ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAAt                                         >  5:17203/1‑148 (MQ=255)
                                                                                           tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCTGTAATGCTGAAATTAc                                      >  7:270133/1‑149 (MQ=255)
                                                                                              tAAATCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgcg                                   >  1:282615/1‑149 (MQ=255)
                                                                                               aaaTCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCgg                                  >  7:64572/1‑149 (MQ=255)
                                                                                                 aTCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGc                                                                                           >  5:23496/1‑90 (MQ=255)
                                                                                                 aTCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGc                                                                                           <  6:23496/90‑1 (MQ=255)
                                                                                                 aTCCGATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGgtg                                >  5:162700/1‑149 (MQ=255)
                                                                                                     gATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                            >  6:140626/1‑149 (MQ=255)
                                                                                                     gATAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                            >  7:120272/1‑149 (MQ=255)
                                                                                                      aTAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTATGCGGTGTGCCg                           >  1:137894/1‑149 (MQ=255)
                                                                                                       tAACCTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCg                           >  2:352708/1‑148 (MQ=255)
                                                                                                           cTGAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTata                       >  7:160737/1‑146 (MQ=255)
                                                                                                             gAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                    >  4:92762/1‑149 (MQ=255)
                                                                                                             gAAAACCGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                    >  4:14784/1‑149 (MQ=255)
                                                                                                                   cGAGATGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc              >  7:182135/1‑149 (MQ=255)
                                                                                                                    gagaTGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAcc             <  2:282615/149‑1 (MQ=255)
                                                                                                                     agaTGGATGAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt            <  5:140626/149‑1 (MQ=255)
                                                                                                                       atggatgAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTATGCGGTGTGCCGTAGACGATAGAACCTTc          <  2:137894/149‑1 (MQ=255)
                                                                                                                         ggatgAGTAAATCCGTAACCACTAATCGGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGACCCTTCCa        >  6:297791/1‑149 (MQ=255)
                                                                                                                           atgAGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg      >  6:69519/1‑149 (MQ=255)
                                                                                                                              aGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTg   <  6:303181/149‑1 (MQ=255)
                                                                                                                              aGTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTg   <  2:335636/149‑1 (MQ=255)
                                                                                                                               gTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTACGGCAGCCAGAGGGCTGCCTTTTTTACAGGTGGTATTGAGAATTGAGACGTGCCGGTAATGCTGAAATTGCGCGGTGTGCCGTAGACGAAGGACGCTTCCAGGTTgg  >  6:269468/1‑149 (MQ=255)
                                                                                                                               gTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg  >  3:208560/1‑149 (MQ=255)
                                                                                                                               gTAAATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg  >  5:56196/1‑149 (MQ=255)
                                                                                                                                  |                                                                                                                                                 
CTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGG  >  NC_000913/1159157‑1159432

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: