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.

ATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑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)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑g                                                                                              >  1:3479015/1‑100 (MQ=255)
       cGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑g                                                                                              >  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)
        gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑a                                                                                          >  1:2580956/1‑100 (MQ=255)
        actaaCCGTTTCGGCACTACTGACCTGTTACGGTTCACGATTTCAGACATGGCTTTTGCCACTGCCCCCTTCCTGATCGAGGCGCTGAATCAgagcat                                                                                                >  2:2010629/5‑96 (MQ=255)
         cggACCGTTTCGGAACTGCTGACCTGTTACCATTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑At                                                                                         >  1:2825382/3‑100 (MQ=255)
         cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑At                                                                                         <  1:2989847/100‑1 (MQ=255)
         cTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑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)
           gACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGc                                                                                       >  1:2502016/1‑100 (MQ=255)
            aCCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCa                                                                                      <  1:664757/100‑1 (MQ=255)
              cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgataatg                                                                                        <  1:2503911/100‑9 (MQ=39)
              cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACg                                                                                    <  1:1581985/100‑1 (MQ=255)
               gTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACgg                                                                                   >  2:2986077/1‑100 (MQ=255)
                 ttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCg                                                                                 <  1:1486778/100‑1 (MQ=255)
                   cGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTa                                                                               >  1:1805333/1‑100 (MQ=255)
                      tcttgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTAttt                                                                            <  2:819564/97‑1 (MQ=255)
                      cACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTAttt                                                                            >  1:136446/1‑100 (MQ=255)
                       aCTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTg                                                                           <  1:1962568/100‑1 (MQ=255)
                         tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGc                                                                         >  2:4186204/1‑100 (MQ=255)
                         tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGc                                                                         <  2:3893737/100‑1 (MQ=255)
                         tgctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGc                                                                         <  2:3081967/100‑1 (MQ=255)
                            tgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTAc                                                                      <  2:376160/100‑1 (MQ=255)
                            tgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTAc                                                                      >  1:235408/1‑100 (MQ=255)
                            tgACCTGTTACCGTTCACGATTTCAGACATGGATTTTCCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTAc                                                                      >  2:2243507/1‑100 (MQ=255)
                             gACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACa                                                                     >  1:1163195/1‑100 (MQ=255)
                              aCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAg                                                                    <  1:3573835/100‑1 (MQ=255)
                               ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGc                                                                   >  1:2189191/1‑100 (MQ=255)
                               ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGc                                                                   >  1:1026421/1‑100 (MQ=255)
                               ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGc                                                                   <  1:2534908/100‑1 (MQ=255)
                                cTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGcc                                                                  <  1:2243512/100‑1 (MQ=255)
                                 tGTTACCGTTCACGATTTCAAACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCg                                                                 >  2:2568271/1‑100 (MQ=255)
                                   ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCt                                                               <  1:3215708/100‑1 (MQ=255)
                                   ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCt                                                               <  1:839688/100‑1 (MQ=255)
                                      ccGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGa                                                            >  1:2890937/1‑100 (MQ=255)
                                       cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGaa                                                           >  2:1476945/1‑100 (MQ=255)
                                       cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGaa                                                           >  2:3179501/1‑100 (MQ=255)
                                            aCGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACg                                                      >  1:3840165/1‑100 (MQ=255)
                                               aTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATg                                                   <  2:964340/100‑1 (MQ=255)
                                                   cAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGtt                                               >  2:4147056/1‑100 (MQ=255)
                                                   cAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGtt                                               <  1:2521472/100‑1 (MQ=255)
                                                     gACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtc                                             >  2:2454198/1‑100 (MQ=255)
                                                       cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                           >  1:1449640/1‑100 (MQ=255)
                                                       cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                           >  1:1554553/1‑100 (MQ=255)
                                                       cATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtctc                                           >  2:2696162/1‑100 (MQ=255)
                                                        aTGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcg                                          >  2:940237/1‑100 (MQ=255)
                                                        aTGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTcg                                          >  1:2211987/1‑100 (MQ=255)
                                                           gATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCgg                                       <  1:399304/100‑1 (MQ=255)
                                                            aTTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGc                                      <  1:1164271/100‑1 (MQ=255)
                                                            aTTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGc                                      <  2:281371/100‑1 (MQ=255)
                                                               ttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGCGTTTCTCGCGGCTAt                                   >  2:2658959/1‑100 (MQ=255)
                                                                 gCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGTAAAACGATGAGTTTCTCGCGGCTAtt                                  >  1:790598/1‑99 (MQ=255)
                                                                 gCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTg                                 >  2:4007983/1‑100 (MQ=255)
                                                                   cACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                               <  1:2205046/100‑1 (MQ=255)
                                                                    aCTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGccc                              <  1:1090852/100‑1 (MQ=255)
                                                                      tGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCAc                            <  1:1889674/100‑1 (MQ=255)
                                                                      tGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCAc                            <  1:2733899/100‑1 (MQ=255)
                                                                         ccccTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                         <  2:1554550/100‑1 (MQ=255)
                                                                         ccccTGCATTAGCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                         <  2:1016023/100‑1 (MQ=255)
                                                                          cccTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGt                        <  2:3840140/100‑1 (MQ=255)
                                                                          atcTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCTCTgga                        >  2:211408/3‑99 (MQ=255)
                                                                           ccTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtt                       >  2:4193909/1‑100 (MQ=255)
                                                                           ccTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtt                       <  1:4147085/100‑1 (MQ=255)
                                                                            cTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                      <  1:65420/100‑1 (MQ=255)
                                                                            cTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                      <  1:524171/100‑1 (MQ=255)
                                                                             gccgtgcatcaGGCGCTGATTCAGCGCCTCGATAATGCAGGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTgg                         >  1:3518953/11‑100 (MQ=12)
                                                                              gCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTc                    >  1:536767/1‑100 (MQ=255)
                                                                              ccgttcatcaGGCGCTGATTCAGCGCCT‑G‑‑‑ATGCCCGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                     >  2:608471/10‑99 (MQ=38)
                                                                                gtgcatcaGGCGCTGATTCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCa                  >  1:3208582/8‑100 (MQ=255)
                                                                                  gcatcaGGCGCTGCTTCCGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTATCGCGCGGCTATCGCCAACAGGTTTCCCAc                 >  2:261208/6‑99 (MQ=255)
                                                                                   aTCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCATGGCGTATTTGGCTACAGCCGCTGTAAAAACGCTGCGTGTCTCGCGGCTATTGCCCACTGGTTTTc                    <  2:2224192/95‑1 (MQ=255)
                                                                                   aTCGAGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAccc               >  2:3914521/1‑100 (MQ=255)
                                                                                    tCGAGGCGCTGAATCAGCGCAT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAACCGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa              <  1:2879543/100‑1 (MQ=255)
                                                                                    atcaGGCGCTGATTCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa              >  1:4010062/4‑100 (MQ=255)
                                                                                       aGGCGCTGATTCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCa           >  2:1480262/1‑100 (MQ=255)
                                                                                       aGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCa           <  1:1710056/100‑1 (MQ=255)
                                                                                       aGGCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATAAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCCGc            >  2:2079024/1‑99 (MQ=255)
                                                                                        ggCGCTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGACTATTGCCCACTGGTTTTCCACCCAGCAt          >  2:1675058/1‑100 (MQ=255)
                                                                                         gcgcTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt         <  2:4100561/100‑1 (MQ=255)
                                                                                          cgcTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa        <  1:1718232/100‑1 (MQ=255)
                                                                                          cgcTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa        >  1:3029045/1‑100 (MQ=255)
                                                                                           gcTGAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTac       >  2:2079510/1‑100 (MQ=255)
                                                                                              gAATCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAcc    >  2:313196/1‑100 (MQ=255)
                                                                                                aTCAGCGCCT‑G‑‑‑ATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGc  <  2:104432/100‑1 (MQ=255)
                                                                                                |                                                                                                       
ATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCT‑G‑‑‑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: ‑