Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 232,196 G→T G91C (GGC→TGC)  yafE → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913232,1960GT96.3% 438.3 / 5.2 135G91C (GGC→TGC) yafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/5);  new base T (62/68);  total (62/73)
Fisher's exact test for biased strand distribution p-value = 6.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 135 total aligned reads displayed.

GATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACG  >  NC_000913/232100‑232292
                                                                                                |                                                                                                
gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                               <  2:1012648/100‑1 (MQ=255)
gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                               <  2:2693381/100‑1 (MQ=255)
    aCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAg                                                                                           >  2:2646414/1‑100 (MQ=255)
     cGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAgg                                                                                          <  2:1048757/100‑1 (MQ=255)
         tttGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGa                                                                                      <  2:2972760/100‑1 (MQ=255)
          ttGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAt                                                                                     >  1:2209012/1‑100 (MQ=255)
          ttGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAt                                                                                     <  2:2493725/100‑1 (MQ=255)
           tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGACTTCACTAGAATCTTGAAACCTTGCGGTAGGCTGa                                                                                      >  1:3285511/1‑98 (MQ=255)
           tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAtt                                                                                    <  1:2155133/100‑1 (MQ=255)
            gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtg                                                                                   <  1:1048462/100‑1 (MQ=255)
              tatTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtgtg                                                                                 >  2:1248730/1‑100 (MQ=255)
               atTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGa                                                                                <  2:3271845/100‑1 (MQ=255)
                ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGa                                                                                <  2:2945815/99‑1 (MQ=255)
                ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                               <  2:3173781/100‑1 (MQ=255)
                ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                               >  1:2022757/1‑100 (MQ=255)
                ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                               >  2:1403640/1‑100 (MQ=255)
                 tGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATg                                                                              <  1:949816/100‑1 (MQ=255)
                 tGTTATAAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCCCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGa                                                                                <  2:509703/98‑1 (MQ=255)
                  gTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATgg                                                                             >  2:1684206/1‑100 (MQ=255)
                   ttATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGa                                                                            <  2:1204189/100‑1 (MQ=255)
                     aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGa                                                                            >  2:1287497/1‑98 (MQ=255)
                       cAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTa                                                                        <  1:1397473/100‑1 (MQ=255)
                       cAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTa                                                                        >  1:707029/1‑100 (MQ=255)
                       cAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTa                                                                        >  2:743702/1‑100 (MQ=255)
                        agccggaagtCTGCCCATCTCTGGGATGATTGGGTTTCAGCACTGCGATAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACg                                                                          <  1:3357671/88‑1 (MQ=255)
                        aGCCGTTATTCTGCCCCTCACTGGCATGATGTTGGTGCAGCACCGACACAAGTGAATAGGATATCGAAACCTTGCGGTAGGCTGCTTGTGACGGACCTac                                                                       >  1:3205354/1‑99 (MQ=255)
                        aGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTaa                                                                       <  1:3232740/100‑1 (MQ=255)
                         gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGTACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      <  1:3176594/100‑1 (MQ=255)
                         gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      >  2:2170184/1‑100 (MQ=255)
                          ccGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGGTATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                     >  2:733602/1‑100 (MQ=255)
                          ccGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                     <  2:1732191/100‑1 (MQ=255)
                           cGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                     >  1:1085126/1‑99 (MQ=255)
                           cGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGt                                                                    >  2:2656739/1‑100 (MQ=255)
                             ttATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGtct                                                                  <  1:226544/100‑1 (MQ=255)
                              tATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGtctc                                                                 >  2:194066/1‑100 (MQ=255)
                               aTTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTcc                                                                <  2:915146/100‑1 (MQ=255)
                                 tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCg                                                               >  1:2549708/1‑99 (MQ=255)
                                  cTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCggg                                                             <  2:330054/100‑1 (MQ=255)
                                     cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCa                                                          >  2:2082071/1‑100 (MQ=255)
                                      ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAc                                                         >  1:490815/1‑100 (MQ=255)
                                      ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAc                                                         >  2:1685429/1‑100 (MQ=255)
                                       cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTACTGTCTCCGGGTCAcc                                                        >  2:3424950/1‑100 (MQ=255)
                                        aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAccc                                                       <  2:1266162/100‑1 (MQ=255)
                                          cACTGTCATGATGTTGGTGCCGCACTTCGAGAAGTGAATAGGAGATTTAAACCTTGCGGTAGGCTGATTGTGTTGGCCGTATTGTCTCCGGGTCAccc                                                       <  1:3274061/98‑1 (MQ=255)
                                            cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   >  2:3413554/1‑100 (MQ=255)
                                            cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   <  1:2239489/100‑1 (MQ=255)
                                            cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   >  2:1803921/1‑100 (MQ=255)
                                            cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   >  1:1139555/1‑100 (MQ=255)
                                            cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   <  2:3112452/100‑1 (MQ=255)
                                              ggCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTgcg                                                 >  2:2126599/1‑100 (MQ=255)
                                              gacATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTgcg                                                 >  1:2069041/3‑100 (MQ=255)
                                                cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                               <  2:2750374/100‑1 (MQ=255)
                                                cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                               <  1:1684183/100‑1 (MQ=255)
                                                cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                               <  2:2513332/100‑1 (MQ=255)
                                                cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                               <  2:1297805/100‑1 (MQ=255)
                                                cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                               >  1:2688172/1‑100 (MQ=255)
                                                cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                               >  1:952007/1‑100 (MQ=255)
                                                  tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                             <  1:1370837/100‑1 (MQ=255)
                                                  tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                             <  1:1370847/100‑1 (MQ=255)
                                                   gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACa                                            >  2:2567865/1‑100 (MQ=255)
                                                    atgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACAt                                           >  1:2405196/1‑100 (MQ=255)
                                                     tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                          <  2:2169974/100‑1 (MQ=255)
                                                     tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                          <  1:2798753/100‑1 (MQ=255)
                                                         ggTGCCGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                      >  1:1206871/1‑100 (MQ=255)
                                                         ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                      <  2:1510305/100‑1 (MQ=255)
                                                           tGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGtt                                     <  1:2384804/99‑1 (MQ=255)
                                                           tGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTa                                    >  1:938800/1‑100 (MQ=255)
                                                           tGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTATCCGGGTCACGCAGTGCGCGACATCTGGTTa                                    >  1:109873/1‑100 (MQ=255)
                                                            gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                   >  1:2595354/1‑100 (MQ=255)
                                                             cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                  <  1:117986/100‑1 (MQ=255)
                                                             cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                  <  1:2698002/100‑1 (MQ=255)
                                                               gcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGa                                >  1:1661719/1‑100 (MQ=255)
                                                                caCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGAc                               <  1:194063/100‑1 (MQ=255)
                                                                 aCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                              >  1:304645/1‑100 (MQ=255)
                                                                 aCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                              <  2:2529521/100‑1 (MQ=255)
                                                                  cggCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                             <  2:2415111/98‑1 (MQ=255)
                                                                  cTGCGAGGAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                             <  1:2646738/100‑1 (MQ=255)
                                                                  cTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgn                             <  1:2480370/100‑2 (MQ=255)
                                                                  cTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                             >  2:1456378/1‑100 (MQ=255)
                                                                    gCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                           >  1:2881557/1‑100 (MQ=255)
                                                                     cGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                          >  1:1487266/1‑100 (MQ=255)
                                                                     cGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                          >  2:2939933/1‑100 (MQ=255)
                                                                     cGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                          >  2:3017466/1‑100 (MQ=255)
                                                                      gagaAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGa                         <  2:547684/100‑1 (MQ=255)
                                                                       agaAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGaa                        <  2:394300/100‑1 (MQ=255)
                                                                        gaAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                       <  1:546570/100‑1 (MQ=255)
                                                                          aGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                     <  1:1799559/100‑1 (MQ=255)
                                                                             gAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                  >  2:3039675/1‑100 (MQ=255)
                                                                              ttccgatctttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                 <  2:9539/91‑1 (MQ=255)
                                                                              aaTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                 <  2:2595389/100‑1 (MQ=255)
                                                                                 aGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgcg              <  2:1278665/100‑1 (MQ=255)
                                                                                   gATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGACGCATTACGCGAt            <  2:2282120/100‑1 (MQ=255)
                                                                                   gATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt            <  2:506767/100‑1 (MQ=255)
                                                                                   gATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt            <  2:1033753/100‑1 (MQ=255)
                                                                                        tGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc       >  1:642324/1‑100 (MQ=255)
                                                                                        tGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc       >  1:466841/1‑100 (MQ=255)
                                                                                           aaCCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCa    <  2:2738204/100‑1 (MQ=255)
                                                                                            aCCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCAc   >  2:66884/1‑100 (MQ=255)
                                                                                             ccTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCtctc     <  2:204894/97‑1 (MQ=255)
                                                                                             ccTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg  >  1:1766750/1‑100 (MQ=255)
                                                                                                |                                                                                                
GATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACG  >  NC_000913/232100‑232292

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