Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 187,676 T→A 6.0% E325D (GAA→GAT)  glnD ← PII uridylyltransferase/uridylyl removing enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913187,6760TA 6.0% 211.7 / 3.1 66E325D (GAA→GAT) glnDPII uridylyltransferase/uridylyl removing enzyme
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base A (2/2);  total (31/35)
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.09e-01

GCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCA  >  NC_000913/187527‑187824
                                                                                                                                                     |                                                                                                                                                    
gattCCGGATGGCGCATAAATAGTGTTTAATAACGCAGGTCGATTAGCGTACCTAGTAGCTGAAACTCATCGTAGATTGGAAGTGGTTTTTCGTCGGCGGGAAGTGCGAGGATGGCTTCATCGAACAGTTTAAGCAGCATCTGGTTGAGt                                                                                                                                                      >  2:682382/3‑150 (MQ=255)
 cTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGtt                                                                                                                                                     >  1:614852/1‑150 (MQ=255)
 cTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGtt                                                                                                                                                     >  2:496383/1‑150 (MQ=255)
 cTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGtt                                                                                                                                                     >  1:237770/1‑150 (MQ=255)
 cTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGtt                                                                                                                                                     >  2:776869/1‑150 (MQ=255)
     cGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACt                                                                                                                                                 >  1:350987/1‑150 (MQ=255)
          ggCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGAcgc                                                                                                                                            >  2:370145/1‑150 (MQ=255)
            cgcATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACgcgg                                                                                                                                          >  2:850587/1‑150 (MQ=255)
              cATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACgcggcg                                                                                                                                        >  1:793775/1‑150 (MQ=255)
              cATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACgcggcg                                                                                                                                        <  2:191958/150‑1 (MQ=255)
                 aaaTAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCgtgt                                                                                                                                     <  1:213159/150‑1 (MQ=255)
                        gtTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCgg                                                                                                                              <  2:744567/150‑1 (MQ=255)
                           tcatcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAg                                                                                                                           <  1:332584/150‑1 (MQ=255)
                           tcatcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGGTGAGTTCACTGACGCGGCGTGTAACGCGGAAg                                                                                                                           <  1:474974/150‑1 (MQ=255)
                            catcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGAGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGt                                                                                                                          >  2:502070/1‑150 (MQ=255)
                               caCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAAt                                                                                                                       >  1:481363/1‑150 (MQ=255)
                                      gTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGTTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTtcat                                                                                                                <  2:149755/150‑1 (MQ=255)
                                         gATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTcatcat                                                                                                             >  2:209593/1‑150 (MQ=255)
                                                gTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCg                                                                                                      <  1:644588/150‑1 (MQ=255)
                                                   ccGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGAcc                                                                                                   >  1:59056/1‑150 (MQ=255)
                                                           cTGAAACTCTTCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGT‑‑‑‑‑‑AGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTc                                                                                     <  2:307511/150‑1 (MQ=255)
                                                            tGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTa                                                                                          >  1:517336/1‑150 (MQ=255)
                                                            tGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTa                                                                                          <  1:690876/150‑1 (MQ=255)
                                                               aaCTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCt                                                                                       <  1:496383/150‑1 (MQ=255)
                                                               aaCTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCt                                                                                       <  1:370145/150‑1 (MQ=255)
                                                                aCTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCTAGTATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTTTAACGCGTAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCtt                                                                                      <  1:682382/150‑1 (MQ=255)
                                                                   cATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTcac                                                                                   >  1:197073/1‑150 (MQ=255)
                                                                    aTCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTcacc                                                                                  >  1:587/1‑150 (MQ=255)
                                                                    aTCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTcacc                                                                                  >  2:454888/1‑150 (MQ=255)
                                                                    aTCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTcacc                                                                                  <  2:481363/150‑1 (MQ=255)
                                                                      cgtcgATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTcaccac                                                                                >  1:35336/1‑150 (MQ=255)
                                                                         cgATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGt                                                                             <  2:350987/150‑1 (MQ=255)
                                                                            ttGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAAt                                                                          <  1:1048799/150‑1 (MQ=255)
                                                                               gACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCGCCACTGTAATTCa                                                                       >  2:1201324/1‑150 (MQ=255)
                                                                                   tGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACg                                                                   >  2:13656/1‑150 (MQ=255)
                                                                                     gTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCt                                                                 >  1:981882/1‑150 (MQ=255)
                                                                                      tttttCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTg                                                                <  2:586555/150‑1 (MQ=255)
                                                                                        tttCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGATCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGAc                                                                                                    <  1:683731/112‑1 (MQ=255)
                                                                                        tttCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGATCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGAc                                                                                                    >  2:683731/1‑112 (MQ=255)
                                                                                        tttCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGATCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGAc                                                                                                    >  2:688735/1‑112 (MQ=255)
                                                                                        tttCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGATCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGAc                                                                                                    <  1:688735/112‑1 (MQ=255)
                                                                                         ttCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGc                                                             <  2:1350115/150‑1 (MQ=255)
                                                                                             tcgGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACg                                                         <  2:35336/150‑1 (MQ=255)
                                                                                                    gAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCt                                                  >  1:424683/1‑150 (MQ=255)
                                                                                                      aGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTgg                                                <  1:502070/150‑1 (MQ=255)
                                                                                                       gggCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGc                                               <  1:209593/150‑1 (MQ=255)
                                                                                                       gggCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGc                                               >  1:544783/1‑150 (MQ=255)
                                                                                                             ggATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGa                                         >  1:555544/1‑150 (MQ=255)
                                                                                                              gATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGaa                                        <  2:517336/150‑1 (MQ=255)
                                                                                                              gATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGaa                                        >  2:754431/1‑150 (MQ=255)
                                                                                                                tGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATa                                      <  1:454888/150‑1 (MQ=255)
                                                                                                                tGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATa                                      <  1:41196/150‑1 (MQ=255)
                                                                                                                     tCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCAGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGCTCGAATAACAgg                                 >  1:51005/1‑150 (MQ=255)
                                                                                                                      cATCGAACAGTTGCAGCAGAATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGc                                <  2:119678/150‑1 (MQ=255)
                                                                                                                       aTCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCg                               >  1:1363115/1‑150 (MQ=255)
                                                                                                                          gAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGAtt                            >  2:1346023/1‑150 (MQ=255)
                                                                                                                            aCAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTAt                          <  2:555544/150‑1 (MQ=255)
                                                                                                                            aCAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTAt                          <  1:487722/150‑1 (MQ=255)
                                                                                                                                  gcagcagcaTCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAAc                    <  1:1036026/150‑1 (MQ=255)
                                                                                                                                   cagcagcaTCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACg                   >  2:99310/1‑150 (MQ=255)
                                                                                                                                     gcagcaTCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGc                 <  1:850587/150‑1 (MQ=255)
                                                                                                                                         caTCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTgac             >  1:782760/1‑150 (MQ=255)
                                                                                                                                            cTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTgacgac          <  1:13656/150‑1 (MQ=255)
                                                                                                                                              ggTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCa        >  1:997131/1‑150 (MQ=255)
                                                                                                                                               gTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAg       >  1:1068785/1‑150 (MQ=255)
                                                                                                                                               gTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAg       >  2:965866/1‑150 (MQ=255)
                                                                                                                                                   aGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCg                        <  2:1001352/129‑1 (MQ=255)
                                                                                                                                                   aGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCg                        >  1:1001352/1‑129 (MQ=255)
                                                                                                                                                   aGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGc   <  2:197073/150‑1 (MQ=255)
                                                                                                                                                    gTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCa  <  2:59056/150‑1 (MQ=255)
                                                                                                                                                    gTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCa  <  1:776869/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCACCACTGTAATTCAGACGCTGGGCGACGCTAAGCTGGCGATCGAATAACAGGCGATTATCGTAACGGCTGACGACCAGATGCA  >  NC_000913/187527‑187824

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