Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,550,199 A→G intergenic (‑93/‑519) malP ← / → malT maltodextrin phosphorylase/MalT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,550,1990AG85.2% 632.6 / 98.5 250intergenic (‑93/‑519)malP/malTmaltodextrin phosphorylase/MalT
Reads supporting (aligned to +/- strand):  ref base A (14/23);  new base G (97/116);  total (111/139)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

Only 100 of 258 total aligned reads displayed.

GGAGTTCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCAATACA  >  NC_000913/3550108‑3550291
                                                                                           |                                                                                            
ggAGTTCCAGTTTTCTTAGATTTCCAACACAACGTTATCGCTAGTTTGCCAGGCTCTATGTTGTCCTTCCTCATCCTGCGGGGGATTAGGCAg                                                                                             <  1:3726553/93‑1 (MQ=255)
 gAGTTCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCCCATCCTGCGGGGGATTAGGCggg                                                                                            <  1:6127550/93‑1 (MQ=255)
  aGTTCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCgggg                                                                                           >  2:2470553/1‑93 (MQ=255)
    ttCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGggag                                                                                         >  2:5503094/1‑93 (MQ=255)
    ttCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCGCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGggag                                                                                         <  2:4672485/93‑1 (MQ=255)
     tCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGggagg                                                                                        <  1:5874719/93‑1 (MQ=255)
     tCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGggagg                                                                                        <  1:1503235/93‑1 (MQ=255)
       cACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGgaggag                                                                                      >  2:1040287/1‑93 (MQ=255)
        aCTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGt                                                                                     <  2:4920652/93‑1 (MQ=255)
        aCTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGt                                                                                     <  2:4190974/93‑1 (MQ=255)
        aCTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGt                                                                                     <  2:4415565/93‑1 (MQ=255)
         cTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGtt                                                                                    >  1:182475/1‑93 (MQ=255)
          ttttCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTg                                                                                   >  2:2622315/1‑93 (MQ=255)
          ttttCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTg                                                                                   <  1:2813560/93‑1 (MQ=255)
           tttCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGc                                                                                  >  2:6351744/1‑93 (MQ=255)
           tttCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGc                                                                                  <  2:2083891/93‑1 (MQ=255)
            ttCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCg                                                                                 <  2:5696488/93‑1 (MQ=255)
             tCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCgg                                                                                <  2:362758/93‑1 (MQ=255)
             tCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCgg                                                                                <  2:3583689/93‑1 (MQ=255)
              cTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCggg                                                                               >  2:3846631/1‑93 (MQ=255)
               ttAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCgggg                                                                              >  1:1656543/1‑93 (MQ=255)
                tAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGa                                                                             <  1:3204072/93‑1 (MQ=255)
                tAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGa                                                                             <  1:893710/93‑1 (MQ=255)
                  gATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATg                                                                           >  1:1019305/1‑93 (MQ=255)
                    ttttCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAg                                                                         <  2:1809870/93‑1 (MQ=255)
                    ttttCAACACAACGTTATCGCTAGGTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAg                                                                         <  2:1855920/93‑1 (MQ=255)
                     tttCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGc                                                                        <  2:5032622/93‑1 (MQ=255)
                       tCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCaa                                                                      <  1:4281093/93‑1 (MQ=255)
                        cAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAg                                                                     <  1:5345873/93‑1 (MQ=255)
                         aaCACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAgg                                                                    <  2:5125011/93‑1 (MQ=255)
                         aaCACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAgg                                                                    >  1:2070714/1‑93 (MQ=255)
                          acacAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGa                                                                   >  1:4131018/1‑93 (MQ=255)
                            agaaCGTAATCGCTAGTTTGCCAGGCTCGATGTTAACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGaaa                                                                 <  2:1013542/91‑1 (MQ=255)
                            acaACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGaaa                                                                 >  2:2536122/1‑93 (MQ=255)
                            acaACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGaaa                                                                 <  1:1424410/93‑1 (MQ=255)
                            acaACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGaaa                                                                 >  1:4168090/1‑93 (MQ=255)
                             caACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAAt                                                                <  2:4583772/93‑1 (MQ=255)
                             caACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAAt                                                                >  2:3241465/1‑93 (MQ=255)
                             caACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAAt                                                                >  2:5626701/1‑93 (MQ=255)
                              aaCGTTATCGCTAGTTTGCCAGGCTCGATGTTTACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAAtg                                                               <  1:854046/93‑1 (MQ=255)
                              aaCGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGGTGAGCAAGGAAAtg                                                               >  2:3503068/1‑93 (MQ=255)
                               aCGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAAtgt                                                              >  1:1621063/1‑93 (MQ=255)
                                   tATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGAtc                                                          <  2:3767882/93‑1 (MQ=255)
                                       gCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCaa                                                      >  1:5762619/1‑93 (MQ=255)
                                        cTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAAc                                                     >  2:5816696/1‑93 (MQ=255)
                                        cTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAAc                                                     <  1:1174676/93‑1 (MQ=255)
                                         tAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAAcc                                                    <  2:114332/93‑1 (MQ=255)
                                          aGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCa                                                   <  1:3472379/93‑1 (MQ=255)
                                          aGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCa                                                   >  1:3948461/1‑93 (MQ=255)
                                           gTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCAc                                                  >  2:5749391/1‑93 (MQ=255)
                                            tttGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACt                                                 <  2:2315698/93‑1 (MQ=255)
                                             ttGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACtt                                                <  1:6090959/93‑1 (MQ=255)
                                               gCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTaa                                              >  1:3378248/1‑93 (MQ=255)
                                               gCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTaa                                              <  2:6183843/93‑1 (MQ=255)
                                                  aGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGc                                           <  1:1351261/93‑1 (MQ=255)
                                                  aGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGc                                           >  2:3367936/1‑93 (MQ=255)
                                                   ggCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCt                                          <  2:6010035/93‑1 (MQ=255)
                                                   ggCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCt                                          >  2:214575/1‑93 (MQ=255)
                                                     cTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAg                                        <  1:4646017/93‑1 (MQ=255)
                                                      tCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGt                                       >  1:1991969/1‑93 (MQ=255)
                                                      tCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGt                                       <  1:4705425/93‑1 (MQ=255)
                                                        gATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGc                                     >  1:1467611/1‑93 (MQ=255)
                                                               aCCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACcac                              <  2:2864189/93‑1 (MQ=255)
                                                               aCCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACcac                              <  1:3834295/93‑1 (MQ=255)
                                                               aCCTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACcac                              >  1:4462916/1‑93 (MQ=255)
                                                                ccTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACcaca                             <  1:6584841/93‑1 (MQ=255)
                                                                 cTTCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAg                            <  2:3090009/93‑1 (MQ=255)
                                                                 cTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGACATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAg                            <  1:6641273/93‑1 (MQ=255)
                                                                  ttCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAgg                           >  1:5648944/1‑93 (MQ=255)
                                                                   tCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGa                          <  1:2483661/93‑1 (MQ=255)
                                                                   tCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGa                          <  1:361792/93‑1 (MQ=255)
                                                                   tCCTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGa                          <  1:3232366/93‑1 (MQ=255)
                                                                    ccTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGAt                         >  1:2202341/1‑93 (MQ=255)
                                                                     cTCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGAtt                        <  1:5911462/93‑1 (MQ=255)
                                                                      tCATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTa                       >  1:5515592/1‑93 (MQ=255)
                                                                       cATCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAg                      <  2:6406160/93‑1 (MQ=255)
                                                                        aTCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGc                     <  1:5319888/93‑1 (MQ=255)
                                                                        aTCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGc                     <  1:2619126/93‑1 (MQ=255)
                                                                         tCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCa                    >  1:3031305/1‑93 (MQ=255)
                                                                         tCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCa                    >  1:5015953/1‑93 (MQ=255)
                                                                         tCCTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAAAGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAAGATTAGCa                    >  1:5120117/1‑93 (MQ=255)
                                                                           cTGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatc                  >  2:13582/1‑93 (MQ=255)
                                                                            tGCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatca                 >  1:4435191/1‑93 (MQ=255)
                                                                             gCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatcaa                <  2:884941/93‑1 (MQ=255)
                                                                             gCGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatcaa                <  2:3904844/93‑1 (MQ=255)
                                                                              cGGGGGATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatcaaa               <  2:4027870/93‑1 (MQ=255)
                                                                               gggggATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatcaaat              <  2:711517/93‑1 (MQ=255)
                                                                                ggggATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatcaaatc             <  2:472462/93‑1 (MQ=255)
                                                                                  ggATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCatcaaatcaa           <  2:5309783/93‑1 (MQ=255)
                                                                                   gATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAAt          >  1:6149455/1‑93 (MQ=255)
                                                                                   gATTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAAt          <  2:774421/93‑1 (MQ=255)
                                                                                    aTTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATg         <  2:4045945/93‑1 (MQ=255)
                                                                                    aTTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATg         <  1:5459122/93‑1 (MQ=255)
                                                                                    aTTAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATg         <  2:6440911/93‑1 (MQ=255)
                                                                                     ttAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGc        <  1:1094366/93‑1 (MQ=255)
                                                                                      tAGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCa       <  1:969560/93‑1 (MQ=255)
                                                                                       aGGCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCaa      <  1:2432015/93‑1 (MQ=255)
                                                                                        ggCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCAAt     >  2:5013467/1‑93 (MQ=255)
                                                                                        ggCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCAAt     <  1:677105/93‑1 (MQ=255)
                                                                                        ggCGGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCAAt     >  1:212478/1‑93 (MQ=255)
                                                                                           |                                                                                            
GGAGTTCCACTTTTCTTAGATTTTCAACACAACGTTATCGCTAGTTTGCCAGGCTCGATGTTGACCTTCCTCATCCTGCGGGGGATTAGGCAGGGAGGAGTTGCGGGGATGAGCAAGGAAATGTGATCTCAACCACTTAAAGCTAGTGCAAACCACAGGATTAGCATCAAATCAATGCAATACA  >  NC_000913/3550108‑3550291

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