Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,108 A→T 15.1% N51I (AAT→ATT)  malY → PLP‑dependent beta‑cystathionase and maltose regulon regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,1080AT15.1% 318.1 / 41.5 119N51I (AAT→ATT) malYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base A (50/51);  new base T (8/10);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 121 total aligned reads displayed.

ATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGC  >  NC_000913/1701012‑1701207
                                                                                                |                                                                                                       
aTTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCa                                                                                                      >  1:2992878/1‑100 (MQ=255)
 ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAg                                                                                                     <  1:2656682/100‑1 (MQ=255)
   aTGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAgcg                                                                                                   <  1:2254334/100‑1 (MQ=255)
     gTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCTCTGCATTATCGAGGCGCTGAATCAgcg                                                                                                   >  2:2932817/1‑98 (MQ=255)
     gTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGcc                                                                                                 >  2:3223678/1‑100 (MQ=255)
     gTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGcc                                                                                                 >  2:2183004/1‑100 (MQ=255)
     gTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTGCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGcc                                                                                                 <  2:2830144/100‑1 (MQ=255)
      tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCt                                                                                                <  1:912138/100‑1 (MQ=255)
      tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCt                                                                                                <  1:853814/100‑1 (MQ=255)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                               <  2:3620551/100‑2 (MQ=255)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                               <  2:3478026/100‑2 (MQ=255)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTg                                                                                              >  1:3479015/1‑100 (MQ=255)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTg                                                                                              >  1:1815424/1‑100 (MQ=255)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCAATATCGAGGCGCTGATTCAGCGCCtc                                                                                               <  2:157790/100‑2 (MQ=255)
        gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcg                                                                                              <  2:226991/100‑3 (MQ=255)
        gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcg                                                                                              <  1:2428570/100‑3 (MQ=255)
        gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcg                                                                                              <  1:1086069/100‑3 (MQ=255)
        gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑a                                                                                          >  1:2580956/1‑100 (MQ=255)
        actaaCCGTTTCGGCACTACTGACCTGTTACGGTTCACGATTTCAGACATGGCTTTTGCCACTGCCCCCTTCCTGATCGAGGCGCTGAATCAgagcat                                                                                                >  2:2010629/5‑96 (MQ=255)
         cggACCGTTTCGGAACTGCTGACCTGTTACCATTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑At                                                                                         >  1:2825382/3‑100 (MQ=255)
         cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑At                                                                                         <  1:2989847/100‑1 (MQ=255)
         cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑At                                                                                         >  1:1256656/1‑100 (MQ=255)
          tGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgat                                                                                            <  1:3999186/100‑5 (MQ=255)
          tGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgat                                                                                            <  1:3713819/100‑5 (MQ=255)
           gACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgata                                                                                           <  2:1725255/100‑6 (MQ=255)
           gACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGc                                                                                       >  1:2502016/1‑100 (MQ=255)
            aCCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCa                                                                                      <  1:664757/100‑1 (MQ=255)
              cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgataatg                                                                                        <  1:2503911/100‑9 (MQ=39)
              cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACg                                                                                    <  1:1581985/100‑1 (MQ=255)
               gTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACgg                                                                                   >  2:2986077/1‑100 (MQ=255)
                 ttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCg                                                                                 <  1:1486778/100‑1 (MQ=255)
                   cGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTa                                                                               >  1:1805333/1‑100 (MQ=255)
                      tcttgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTAttt                                                                            <  2:819564/97‑1 (MQ=255)
                      cACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTAttt                                                                            >  1:136446/1‑100 (MQ=255)
                       aCTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTg                                                                           <  1:1962568/100‑1 (MQ=255)
                         tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGc                                                                         >  2:4186204/1‑100 (MQ=255)
                         tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGc                                                                         <  2:3893737/100‑1 (MQ=255)
                         tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGc                                                                         <  2:3081967/100‑1 (MQ=255)
                            tgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTAc                                                                      <  2:376160/100‑1 (MQ=255)
                            tgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTAc                                                                      >  1:235408/1‑100 (MQ=255)
                            tgACCTGTTACCGTTCACGATTTCAGACATGGATTTTCCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTAc                                                                      >  2:2243507/1‑100 (MQ=255)
                             gACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACa                                                                     >  1:1163195/1‑100 (MQ=255)
                              aCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAg                                                                    <  1:3573835/100‑1 (MQ=255)
                               ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGc                                                                   >  1:2189191/1‑100 (MQ=255)
                               ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGc                                                                   >  1:1026421/1‑100 (MQ=255)
                               ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGc                                                                   <  1:2534908/100‑1 (MQ=255)
                                cTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGcc                                                                  <  1:2243512/100‑1 (MQ=255)
                                 tGTTACCGTTCACGATTTCAAACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCg                                                                 >  2:2568271/1‑100 (MQ=255)
                                   ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCt                                                               <  1:3215708/100‑1 (MQ=255)
                                   ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCt                                                               <  1:839688/100‑1 (MQ=255)
                                      ccGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGa                                                            >  1:2890937/1‑100 (MQ=255)
                                       cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGaa                                                           >  2:1476945/1‑100 (MQ=255)
                                       cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGaa                                                           >  2:3179501/1‑100 (MQ=255)
                                            aCGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACg                                                      >  1:3840165/1‑100 (MQ=255)
                                               aTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATg                                                   <  2:964340/100‑1 (MQ=255)
                                                   cAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGtt                                               >  2:4147056/1‑100 (MQ=255)
                                                   cAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGtt                                               <  1:2521472/100‑1 (MQ=255)
                                                     gACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtc                                             >  2:2454198/1‑100 (MQ=255)
                                                       cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                           >  1:1449640/1‑100 (MQ=255)
                                                       cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                           >  1:1554553/1‑100 (MQ=255)
                                                       cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                           >  2:2696162/1‑100 (MQ=255)
                                                        aTGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcg                                          >  2:940237/1‑100 (MQ=255)
                                                        aTGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcg                                          >  1:2211987/1‑100 (MQ=255)
                                                           gATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCgg                                       <  1:399304/100‑1 (MQ=255)
                                                            aTTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGc                                      <  1:1164271/100‑1 (MQ=255)
                                                            aTTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGc                                      <  2:281371/100‑1 (MQ=255)
                                                               ttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGCGTTTCTCGCGGCTAt                                   >  2:2658959/1‑100 (MQ=255)
                                                                 gCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGTAAAACGATGAGTTTCTCGCGGCTAtt                                  >  1:790598/1‑99 (MQ=255)
                                                                 gCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTg                                 >  2:4007983/1‑100 (MQ=255)
                                                                   cACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                               <  1:2205046/100‑1 (MQ=255)
                                                                    aCTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                              <  1:1090852/100‑1 (MQ=255)
                                                                      tGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCAc                            <  1:1889674/100‑1 (MQ=255)
                                                                      tGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCAc                            <  1:2733899/100‑1 (MQ=255)
                                                                         ccccTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                         <  2:1554550/100‑1 (MQ=255)
                                                                         ccccTGCATTAGCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                         <  2:1016023/100‑1 (MQ=255)
                                                                          cccTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGt                        <  2:3840140/100‑1 (MQ=255)
                                                                          atcTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCTCTgga                        >  2:211408/3‑99 (MQ=255)
                                                                           ccTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtt                       >  2:4193909/1‑100 (MQ=255)
                                                                           ccTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtt                       <  1:4147085/100‑1 (MQ=255)
                                                                            cTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                      <  1:65420/100‑1 (MQ=255)
                                                                            cTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                      <  1:524171/100‑1 (MQ=255)
                                                                             gccgtgcatcaGGCGCTGATTCAGCGCCTCGATAATGCAGGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                         >  1:3518953/11‑100 (MQ=12)
                                                                              gCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTc                    >  1:536767/1‑100 (MQ=255)
                                                                              ccgttcatcaGGCGCTGATTCAGCGCCTG‑‑‑ATGCCCGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                     >  2:608471/10‑99 (MQ=38)
                                                                                gtgcatcaGGCGCTGATTCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCa                  >  1:3208582/8‑100 (MQ=255)
                                                                                  gcatcaGGCGCTGCTTCCGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTATCGCGCGGCTATCGCCAACAGGTTTCCCAc                 >  2:261208/6‑99 (MQ=255)
                                                                                   aTCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCATGGCGTATTTGGCTACAGCCGCTGTAAAAACGCTGCGTGTCTCGCGGCTATTGCCCACTGGTTTTc                    <  2:2224192/95‑1 (MQ=255)
                                                                                   aTCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAccc               >  2:3914521/1‑100 (MQ=255)
                                                                                    tCGAGGCGCTGAATCAGCGCATG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAACCGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa              <  1:2879543/100‑1 (MQ=255)
                                                                                    atcaGGCGCTGATTCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa              >  1:4010062/4‑100 (MQ=255)
                                                                                       aGGCGCTGATTCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCa           >  2:1480262/1‑100 (MQ=255)
                                                                                       aGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCa           <  1:1710056/100‑1 (MQ=255)
                                                                                       aGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATAAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCCGc            >  2:2079024/1‑99 (MQ=255)
                                                                                        ggCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGACTATTGCCCACTGGTTTTCCACCCAGCAt          >  2:1675058/1‑100 (MQ=255)
                                                                                         gcgcTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt         <  2:4100561/100‑1 (MQ=255)
                                                                                          cgcTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa        <  1:1718232/100‑1 (MQ=255)
                                                                                          cgcTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa        >  1:3029045/1‑100 (MQ=255)
                                                                                           gcTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTac       >  2:2079510/1‑100 (MQ=255)
                                                                                              gAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAcc    >  2:313196/1‑100 (MQ=255)
                                                                                                aTCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGc  <  2:104432/100‑1 (MQ=255)
                                                                                                |                                                                                                       
ATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTG‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGC  >  NC_000913/1701012‑1701207

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