Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,238,265 T→C Y35C (TAC→TGC)  mglB ← galactose‑binding transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,238,2650TC99.6% 847.9 / NA 243Y35C (TAC→TGC) mglBgalactose‑binding transport protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (126/116);  total (127/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

TGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAA  >  NC_000913/2238175‑2238353
                                                                                          |                                                                                        
tGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCAc                                                                                        <  2:403208/93‑1 (MQ=255)
  agaATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACtt                                                                                      >  1:4164197/1‑93 (MQ=255)
  agaATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACtt                                                                                      <  1:3266936/93‑1 (MQ=255)
  agaATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACtt                                                                                      >  2:115943/1‑93 (MQ=255)
   gaATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACCGACATAAAGTTATCGTCGCACTta                                                                                     >  2:5549448/1‑93 (MQ=255)
    aaTCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTtat                                                                                    <  1:1639939/93‑1 (MQ=255)
    aaTCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTtat                                                                                    <  2:2233898/93‑1 (MQ=255)
    aaTCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTtat                                                                                    >  2:179757/1‑93 (MQ=255)
       cattcatCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGa                                                                                 >  1:5178796/1‑93 (MQ=255)
       cattcatCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGa                                                                                 <  2:4393349/93‑1 (MQ=255)
        attcatCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGAt                                                                                <  1:5401631/93‑1 (MQ=255)
        attcatCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGAt                                                                                >  1:3849777/1‑93 (MQ=255)
         ttcatCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGAtt                                                                               >  1:2519574/1‑93 (MQ=255)
          tcatcaGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTg                                                                              <  2:4156081/93‑1 (MQ=255)
            atcaGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGtt                                                                            <  1:6178188/93‑1 (MQ=255)
             tcaGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTa                                                                           >  1:5611576/1‑93 (MQ=255)
             tcaGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTa                                                                           <  1:3792245/93‑1 (MQ=255)
              cagcagCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTac                                                                          >  1:1504824/1‑93 (MQ=255)
               agcagcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTaca                                                                         <  1:939538/93‑1 (MQ=255)
               agcagcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTaca                                                                         >  2:4817650/1‑93 (MQ=255)
                gcagcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTacac                                                                        >  1:4566916/1‑93 (MQ=255)
                 cagcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACAcc                                                                       <  1:5537367/93‑1 (MQ=255)
                 cagcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACAcc                                                                       >  2:3454391/1‑93 (MQ=255)
                  agcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCa                                                                      >  2:1082364/1‑93 (MQ=255)
                  agcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCa                                                                      >  1:6094850/1‑93 (MQ=255)
                  agcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCa                                                                      >  1:4054297/1‑93 (MQ=255)
                   gcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCaa                                                                     <  2:6096091/93‑1 (MQ=255)
                    cTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAAt                                                                    >  1:4015414/1‑93 (MQ=255)
                     tGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATg                                                                   <  2:5827276/93‑1 (MQ=255)
                     tGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATg                                                                   <  2:1749872/93‑1 (MQ=255)
                     tGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATg                                                                   >  1:4884373/1‑93 (MQ=255)
                     tGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAGTTAACGGCGCACTTATAGATTGTTACAACAATGc                                                                  >  2:1457300/1‑93 (MQ=255)
                        aCATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGa                                                                <  2:3129646/93‑1 (MQ=255)
                          aTCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGt                                                              >  1:256916/1‑93 (MQ=255)
                           tCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTa                                                             >  1:440649/1‑93 (MQ=255)
                           tCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTa                                                             >  1:2646917/1‑93 (MQ=255)
                                 gcgGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATcagca                                                       >  2:13106/1‑93 (MQ=255)
                                 gcgGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATcagca                                                       <  2:917021/93‑1 (MQ=255)
                                   ggCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCagcagc                                                     <  1:245417/93‑1 (MQ=255)
                                    gCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCg                                                    <  1:4703629/93‑1 (MQ=255)
                                      tttCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCgtg                                                  >  1:2401436/1‑93 (MQ=255)
                                      tttCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCgtg                                                  <  2:505583/93‑1 (MQ=255)
                                       ttCGCATCTTGCTCAATAGCCGTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCgtgt                                                 <  1:1008624/93‑1 (MQ=255)
                                        tCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGtgtg                                                <  2:3256676/93‑1 (MQ=255)
                                          gCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCa                                              <  1:2494052/93‑1 (MQ=255)
                                          gCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCa                                              >  1:6448357/1‑93 (MQ=255)
                                          gCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCa                                              >  2:4859628/1‑93 (MQ=255)
                                           cATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAg                                             >  1:5563959/1‑93 (MQ=255)
                                            aTCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGc                                            >  2:4488071/1‑93 (MQ=255)
                                             tCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCg                                           <  1:695315/93‑1 (MQ=255)
                                             tCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCg                                           >  2:2759751/1‑93 (MQ=255)
                                             tCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCg                                           <  2:407504/93‑1 (MQ=255)
                                               ttGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGc                                         >  1:3749352/1‑93 (MQ=255)
                                                  cTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCAcc                                      <  1:2759433/93‑1 (MQ=255)
                                                  cTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCAcc                                      >  2:5179552/1‑93 (MQ=255)
                                                   tCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCg                                     >  2:1000915/1‑93 (MQ=255)
                                                   tCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCg                                     >  2:15600/1‑93 (MQ=255)
                                                    cAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGa                                    >  1:4493727/1‑93 (MQ=255)
                                                     aaTAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGaa                                   <  1:5296806/93‑1 (MQ=255)
                                                     aaTAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGaa                                   <  1:3385079/93‑1 (MQ=255)
                                                      aTAGCCTTGCGCACTACAGACATAAAGTTATCATCGCACTCATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAAt                                  >  2:742691/1‑93 (MQ=255)
                                                       tAGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATa                                 <  2:2468817/93‑1 (MQ=255)
                                                        aGCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATaa                                >  2:3631475/1‑93 (MQ=255)
                                                         gCCTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAAc                               >  1:1948467/1‑93 (MQ=255)
                                                           cTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACAt                             >  2:1523870/1‑93 (MQ=255)
                                                           cTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACAt                             <  2:5107571/93‑1 (MQ=255)
                                                           cTTGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACAt                             >  2:6459533/1‑93 (MQ=255)
                                                             tGCGCACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGc                           <  2:4303256/93‑1 (MQ=255)
                                                              gcgcACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCt                          <  2:5916342/93‑1 (MQ=255)
                                                               cgcACTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTg                         <  1:2274325/93‑1 (MQ=255)
                                                                  aCTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGcc                      <  1:3922938/93‑1 (MQ=255)
                                                                   cTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCa                     <  1:2715776/93‑1 (MQ=255)
                                                                   cTACAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCa                     <  1:2987539/93‑1 (MQ=255)
                                                                     aCAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATc                   >  2:1850220/1‑93 (MQ=255)
                                                                     aCAGACATAAAGTTATCGTCGCACTTATAGATTGTTAAACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACGGAATAAAATGCTGGCCATc                   >  1:2361995/1‑93 (MQ=255)
                                                                      cAGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATca                  >  2:6155711/1‑93 (MQ=255)
                                                                       aGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATcac                 >  1:6283984/1‑93 (MQ=255)
                                                                       aGACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATcac                 >  1:341701/1‑93 (MQ=255)
                                                                        gACATAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATcaca                <  2:1028322/93‑1 (MQ=255)
                                                                            tAAAGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCAcagcag            >  2:1892704/1‑93 (MQ=255)
                                                                             aaaGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGa           <  1:5908486/93‑1 (MQ=255)
                                                                             aaaGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGa           >  2:164041/1‑93 (MQ=255)
                                                                               aGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACa         <  1:1295070/93‑1 (MQ=255)
                                                                               aGTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACa         >  2:5483008/1‑93 (MQ=255)
                                                                                gTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAg        >  2:3726019/1‑93 (MQ=255)
                                                                                gTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAg        <  2:6247227/93‑1 (MQ=255)
                                                                                gTTATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAg        <  2:2950914/93‑1 (MQ=255)
                                                                                 ttATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAgg       <  1:1674490/93‑1 (MQ=255)
                                                                                  tATCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAggg      <  1:2106044/93‑1 (MQ=255)
                                                                                   aTCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGt     <  2:4990486/93‑1 (MQ=255)
                                                                                   aTCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGt     <  2:3029404/93‑1 (MQ=255)
                                                                                   aTCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGt     >  2:3559842/1‑93 (MQ=255)
                                                                                   aTCGTCGCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGt     <  2:6402529/93‑1 (MQ=255)
                                                                                    tcgtcgCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGtt    <  2:4649734/93‑1 (MQ=255)
                                                                                    tcgtcgCACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGACATCACAGCAGACAggggt    >  1:566738/1‑91 (MQ=255)
                                                                                     cgtcgcACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTa   <  1:610011/93‑1 (MQ=255)
                                                                                     cgtcgcACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTa   <  2:4757398/93‑1 (MQ=255)
                                                                                     cgtcgcACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTa   <  1:1136915/93‑1 (MQ=255)
                                                                                      gtcgcACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTaa  >  2:3907563/1‑93 (MQ=255)
                                                                                      gtcgcACTTATAGATTGTTACACCAATGCGAGGATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATGACAGCAGACAGGGTTaa  >  2:1695450/1‑93 (MQ=255)
                                                                                          |                                                                                        
TGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAA  >  NC_000913/2238175‑2238353

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