Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,715,439 T→C 55.6% L6S (TTA→TCA)  eamB → cysteine and O‑acetylserine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,4390TC55.6% 57.3 / 184.2 152L6S (TTA→TCA) eamBcysteine and O‑acetylserine exporter
Reads supporting (aligned to +/- strand):  ref base T (27/40);  new base C (39/45);  total (67/85)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 155 total aligned reads displayed.

ATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGT  >  NC_000913/2715340‑2715533
                                                                                                   |                                                                                              
atatCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCtttt                                                                                                <  1:3401936/100‑1 (MQ=255)
   tCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAg                                                                                             <  1:3334990/100‑1 (MQ=255)
   tCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAg                                                                                             <  2:1907266/100‑1 (MQ=255)
    cGAATCTAGAATATAAAAACATTCATTTTTTTTTTTTTTTTTTTTAGGGTACTGTCTAACAAAACATAGGAGATAACAAGTGACACCGACCTttttca                                                                                              >  2:1283474/1‑96 (MQ=255)
     gattCTATAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTCTCTACCAAAACAGAGGAGATAACAATTGACACCGACCCTTTCAAGTg                                                                                           <  2:3990333/97‑1 (MQ=255)
     gAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTg                                                                                           >  1:158059/1‑100 (MQ=255)
     gAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTg                                                                                           >  1:2582516/1‑100 (MQ=255)
        tCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCtt                                                                                        <  1:3449760/100‑1 (MQ=255)
         cTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCttt                                                                                       <  1:395439/100‑1 (MQ=255)
         cTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCttt                                                                                       >  1:607594/1‑100 (MQ=255)
         cTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCttt                                                                                       >  1:3195504/1‑100 (MQ=255)
          tAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCtttt                                                                                      >  1:3134297/1‑100 (MQ=255)
          tAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGATACTGTCTAGCAAAACAGAGGAGATTACAAGTGAGACCGACCCTTTCAAGAGCtttt                                                                                      >  2:2789206/1‑100 (MQ=255)
           aGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCttttt                                                                                     <  2:1703422/100‑1 (MQ=255)
            gAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTg                                                                                    <  2:1254271/100‑1 (MQ=255)
              atatAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGCTAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGa                                                                                  <  2:2218479/100‑1 (MQ=255)
               tataAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGAc                                                                                 >  1:2674769/1‑100 (MQ=255)
                ataAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACGCTCTACCAAAACAGAGGAGATAACAAGTGACACCGAACCTTTCAAGTGATTTTTGGAc                                                                                 <  2:1800117/99‑1 (MQ=255)
                   aaaaCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTac                                                                             >  2:1439581/1‑100 (MQ=255)
                   aaaaCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTac                                                                             >  2:2327374/1‑100 (MQ=255)
                    aaaCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTaca                                                                            <  1:407226/100‑1 (MQ=255)
                     aaCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTacac                                                                           <  1:3069458/100‑1 (MQ=255)
                      aCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTacatc                                                                          <  1:2809574/100‑3 (MQ=255)
                      aCATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACAcc                                                                          <  1:2393030/100‑1 (MQ=255)
                      aCATTCATTTTTTTAAATGTTCCGTATCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTAACACCGACCCTTTTAAGTGCTTTTTGGACTTacac                                                                           <  2:688502/99‑1 (MQ=255)
                       cattcattTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACAccc                                                                         <  2:1099803/100‑1 (MQ=255)
                       cattcattTTTTTAAATGTTACGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAGGTGCTTTTTGGACTTATAcc                                                                          >  2:3670472/1‑99 (MQ=255)
                         ttcattTGTTTAAAGGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGATATAACAAGTGACACCCACCCTTTCAAGTGCTTTTTGGACTTACACCCTg                                                                       <  2:2480329/100‑1 (MQ=255)
                          tcattTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAATTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGa                                                                      <  2:2280508/100‑1 (MQ=255)
                           cattTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGTAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGAt                                                                     >  1:1259738/1‑100 (MQ=255)
                           cattTTTATAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGa                                                                      <  2:2570379/99‑1 (MQ=255)
                            attTTTTTAAATGTTCTGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGAtt                                                                    >  2:2464115/1‑100 (MQ=255)
                            attTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGAtt                                                                    >  1:474146/1‑100 (MQ=255)
                             tttttttAAATGTTCCGTGTCGGGTACTGTCTACGAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTa                                                                   <  2:2177656/100‑1 (MQ=255)
                                 tttAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGc                                                               >  2:3024835/1‑100 (MQ=255)
                                 tttAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGc                                                               <  2:2244773/100‑1 (MQ=255)
                                 tttAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGc                                                               >  2:3814894/1‑100 (MQ=255)
                                   tAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTa                                                             >  1:846198/1‑100 (MQ=255)
                                     aaTGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATg                                                           >  2:857903/1‑100 (MQ=255)
                                     aaTGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGCGACACCGACCCTTTCAAGTGCTTGTTATCCTTACACCCTGGTTACCGGTATg                                                           >  2:4110966/1‑100 (MQ=255)
                                      aTGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGa                                                          <  2:1913730/100‑1 (MQ=255)
                                       tGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGAc                                                         >  2:3252497/1‑100 (MQ=255)
                                       tGTTCCGTCTCGGGTACTGTCTACCAAAACAGTGGAGATAACAAGTGACACCGACCCTTTTAAGCGCTTTTTGGACTTACACCCTGATTACCGCTATg                                                           <  2:2719604/98‑1 (MQ=255)
                                        gTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACg                                                        >  1:254149/1‑100 (MQ=255)
                                        gTTCCATGTCGGGTACTGTCTACCAACACGGAGGCGCTAACAAGTGACACCGACCCCTTGAAGTGCTTTATTGACTTACACCCTGATTACCGCTATCACg                                                        >  2:1803190/1‑100 (MQ=255)
                                          tCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGcc                                                      <  2:1277233/100‑1 (MQ=255)
                                           ccGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCa                                                     <  2:2071183/100‑1 (MQ=255)
                                            cGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAg                                                    <  2:2189591/100‑1 (MQ=255)
                                             gtgtCTGGTACTGTCTACCAAAACAGATGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAg                                                    <  2:1077987/99‑1 (MQ=255)
                                                 cGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCg                                               >  2:4130442/1‑100 (MQ=255)
                                                  gggTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGa                                              >  1:3612011/1‑100 (MQ=255)
                                                  gggTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGa                                              >  2:1462550/1‑100 (MQ=255)
                                                     tACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACa                                           <  2:2287990/100‑1 (MQ=255)
                                                      aCTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACaa                                          <  1:2223922/100‑1 (MQ=255)
                                                         gTCTACCAAAACAGTGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTAGCGCTATGACGCCAGGACCGAAGAatat                                       <  2:1627334/100‑1 (MQ=255)
                                                         gTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAatat                                       >  1:959224/1‑100 (MQ=255)
                                                         gTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAatat                                       >  1:361723/1‑100 (MQ=255)
                                                         gTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAatat                                       <  1:3908983/100‑1 (MQ=255)
                                                         gTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAatat                                       <  2:1537849/100‑1 (MQ=255)
                                                         gTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGCGCTTTTTGGACATACACCCTGATTACCGCTATGACGCCAGGAACGAACAatat                                       >  2:1181909/1‑100 (MQ=255)
                                                          tCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATAtt                                      >  1:219608/1‑100 (MQ=255)
                                                           cTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtc                                     <  1:805210/100‑1 (MQ=255)
                                                           cTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATAtt                                      >  2:1072101/1‑99 (MQ=255)
                                                           cTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtc                                     >  1:3071416/1‑100 (MQ=255)
                                                           cTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtc                                     >  1:1067368/1‑100 (MQ=255)
                                                            tACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtct                                    >  1:1802352/1‑100 (MQ=255)
                                                             aCCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtctc                                   >  1:1965922/1‑100 (MQ=255)
                                                             aCCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtcn                                    <  2:3598833/99‑2 (MQ=255)
                                                             aCCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATtctc                                   >  2:3988620/1‑100 (MQ=255)
                                                              ccAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACACTATTCTCg                                  >  2:456442/1‑100 (MQ=255)
                                                              ccAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCg                                  <  2:2262089/100‑1 (MQ=255)
                                                              ccAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCg                                  >  2:642924/1‑100 (MQ=255)
                                                              ccAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCg                                  >  1:3642449/1‑100 (MQ=255)
                                                              ccAAAACAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCg                                  >  2:1353754/1‑100 (MQ=255)
                                                                aaaaCAGAGGAGATAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGcc                                <  2:3134278/100‑1 (MQ=255)
                                                                        ggagATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGctc                        >  2:240238/1‑100 (MQ=255)
                                                                        ggagATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGctc                        <  2:1224869/100‑1 (MQ=255)
                                                                         gagaTAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGctcn                       <  2:453342/100‑2 (MQ=255)
                                                                          agaTAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTg                      <  2:2986382/100‑1 (MQ=255)
                                                                          agaTAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTg                      >  1:3851365/1‑100 (MQ=255)
                                                                           gaTAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGc                     >  2:1082976/1‑100 (MQ=255)
                                                                            aTAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCt                    >  1:1351728/1‑100 (MQ=255)
                                                                             tAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTa                   <  1:395662/100‑1 (MQ=255)
                                                                             tAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTa                   <  2:3691889/100‑1 (MQ=255)
                                                                             tAACAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTa                   <  1:2023548/100‑1 (MQ=255)
                                                                              aaCAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTAc                  <  1:1329800/100‑1 (MQ=255)
                                                                                cAAGTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGt                >  1:2622658/1‑100 (MQ=255)
                                                                                  aGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCg              <  1:250655/100‑1 (MQ=255)
                                                                                   gTGACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGc             >  2:240077/1‑100 (MQ=255)
                                                                                     gACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCAt           >  2:1340120/1‑100 (MQ=255)
                                                                                     gACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCAt           <  1:2770337/100‑1 (MQ=255)
                                                                                     gACACCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCAt           >  1:694406/1‑100 (MQ=255)
                                                                                        acCGACCCTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGa        <  1:935668/100‑1 (MQ=255)
                                                                                         ccGACCCTTTCAAGTGCGGTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGAt       >  2:3708999/1‑100 (MQ=255)
                                                                                           ctcccTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGAttt     <  1:3010025/98‑1 (MQ=255)
                                                                                             cccTTTTAAGTGCTTTTTTTTTTTTTTTTCTGATTACCGCTATGACGCCAGTACTCAACTATATTCTCGCCCTTAGCTCTGCTACGTCGCATTGCttt     >  2:290739/1‑98 (MQ=255)
                                                                                             cccTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCg   >  1:4074513/1‑100 (MQ=255)
                                                                                             cccTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCg   >  1:3765258/1‑100 (MQ=255)
                                                                                             cccTTTCCAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCg   <  2:3242723/100‑1 (MQ=255)
                                                                                             cccTTTCAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCg   <  1:2058818/100‑1 (MQ=255)
                                                                                                   |                                                                                              
ATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGTGCTTTTTGGACTTACACCCTGATTACCGCTATGACGCCAGGACCGAACAATATTCTCGCCCTTAGCTCTGCTACGTCGCATGGATTTCGT  >  NC_000913/2715340‑2715533

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: