Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,291,463 A→G K262E (AAG→GAG)  galU → glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,291,4630AG99.2% 368.6 / ‑4.4 118K262E (AAG→GAG) galUglucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (52/65);  minor base T (1/0);  total (53/65)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 125 total aligned reads displayed.

CTCCGGGAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGAAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAA  >  NC_000913/1291374‑1291552
                                                                                         |                                                                                         
cTCCGGGAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGgaga                                                                                        >  2:4063578/1‑93 (MQ=255)
cTCCGGGAGCTGGTAATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGgaga                                                                                        <  2:4256845/93‑1 (MQ=255)
 tCCGGGAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGagag                                                                                       <  1:619839/93‑1 (MQ=255)
    gggAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCa                                                                                    <  1:1230246/93‑1 (MQ=255)
    gggAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCa                                                                                    <  1:1512134/93‑1 (MQ=255)
    gggAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCa                                                                                    >  2:1661362/1‑93 (MQ=255)
    gggAGCTGGTGATGAAAGTCAGCTTACCGACGCAATGGATATGCTATTCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCa                                                                                    <  2:1426372/93‑1 (MQ=255)
       aGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGa                                                                                 >  2:1420831/1‑93 (MQ=255)
        gCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGAc                                                                                >  2:3039986/1‑93 (MQ=255)
         cTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACt                                                                               >  1:1287066/1‑93 (MQ=255)
         cTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACt                                                                               >  1:1750400/1‑93 (MQ=255)
           ggTAATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGc                                                                             <  2:3745308/93‑1 (MQ=255)
            gTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCg                                                                            <  1:1802044/93‑1 (MQ=255)
              gatgaAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGt                                                                          >  1:2981185/1‑93 (MQ=255)
               ctgaAATTCATCTCGCCGACGCAATTGATATGCTGAGCGAAAAAGAACCGGTGGAAGCCTATCATATGAAAGGGGGGAGCCATGACTGCGGta                                                                         <  2:2380283/92‑1 (MQ=255)
                tgaAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGtaa                                                                        <  1:2794992/93‑1 (MQ=255)
                tgaAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGtaa                                                                        >  2:1942558/1‑93 (MQ=255)
                    aTTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATaaa                                                                    >  2:2293761/1‑93 (MQ=255)
                    aTTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATaaa                                                                    >  1:3057358/1‑93 (MQ=255)
                     ttCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAAt                                                                   <  2:2049645/93‑1 (MQ=255)
                     ttCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAAt                                                                   <  2:353037/93‑1 (MQ=255)
                     ttCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAAt                                                                   >  1:4599758/1‑93 (MQ=255)
                       cAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTa                                                                 <  1:3881432/93‑1 (MQ=255)
                       cAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTa                                                                 <  2:3608001/93‑1 (MQ=255)
                        aGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAg                                                                <  2:681897/93‑1 (MQ=255)
                         gCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAgg                                                               >  2:613951/1‑93 (MQ=255)
                         gCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAgg                                                               >  1:3433929/1‑93 (MQ=255)
                         gCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAgg                                                               <  2:467138/93‑1 (MQ=255)
                          cTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGt                                                              <  1:2354531/93‑1 (MQ=255)
                           tCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGtt                                                             <  1:2491739/93‑1 (MQ=255)
                           tCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGtt                                                             <  2:4331453/93‑1 (MQ=255)
                             aCCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTAc                                                           <  2:3726352/93‑1 (MQ=255)
                             aCCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTAc                                                           <  2:3712795/93‑1 (MQ=255)
                              ccGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACa                                                          <  1:3566130/93‑1 (MQ=255)
                                  cGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCa                                                      >  1:1938716/1‑93 (MQ=255)
                                   gCAACTGATATGCTGATCGAAAAAGAACCGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAg                                                     <  2:2308425/93‑1 (MQ=255)
                                    cAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAgg                                                    >  1:2310141/1‑93 (MQ=255)
                                      aTTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGcc                                                  >  2:2683179/1‑93 (MQ=255)
                                      aTTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGcc                                                  >  1:2048493/1‑93 (MQ=255)
                                      aTTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGcc                                                  >  1:1391479/1‑93 (MQ=255)
                                      aTTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGcc                                                  <  2:3250748/93‑1 (MQ=255)
                                       ttGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCt                                                 >  1:48339/1‑93 (MQ=255)
                                         gATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTc                                               >  2:2710439/1‑93 (MQ=255)
                                          atatGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCg                                              >  1:1481141/1‑93 (MQ=255)
                                           tatGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGt                                             >  2:1049098/1‑93 (MQ=255)
                                            atGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGtt                                            <  2:2609933/93‑1 (MQ=255)
                                              gCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGa                                          >  2:2195198/1‑93 (MQ=255)
                                              gCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGa                                          >  1:3804216/1‑93 (MQ=255)
                                              gCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGa                                          <  1:3730708/93‑1 (MQ=255)
                                               ctaaTCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGaa                                         <  2:1927515/90‑1 (MQ=255)
                                               cTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGaa                                         <  2:2224756/93‑1 (MQ=255)
                                               cTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGaa                                         <  1:3924180/93‑1 (MQ=255)
                                                tGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAAt                                        >  2:1000726/1‑93 (MQ=255)
                                                tGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAAt                                        >  2:2307692/1‑93 (MQ=255)
                                                 gATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATa                                       <  2:2129497/93‑1 (MQ=255)
                                                 gATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATa                                       >  2:3705947/1‑93 (MQ=255)
                                                 gATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGCCCTTCGTTGAATa                                       <  1:1017650/93‑1 (MQ=255)
                                                  aTCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGACTTCGTTGAATAc                                      >  1:4310595/1‑93 (MQ=255)
                                                   acGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACg                                     <  1:3508830/92‑1 (MQ=255)
                                                    cGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACgg                                    <  2:276727/93‑1 (MQ=255)
                                                    cGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACgg                                    <  2:4059401/93‑1 (MQ=255)
                                                     gAAAAAGAAACGGTGGAAGCGTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGt                                   <  2:1846015/93‑1 (MQ=255)
                                                     gAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGt                                   <  1:1688676/93‑1 (MQ=255)
                                                     gAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGt                                   >  2:694162/1‑93 (MQ=255)
                                                      aaaaaGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTa                                  <  1:472532/93‑1 (MQ=255)
                                                       aaaaGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTAt                                 >  1:3474503/1‑93 (MQ=255)
                                                        aaaGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTAtt                                <  2:3320610/93‑1 (MQ=255)
                                                        aaaGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTAtt                                <  1:3715715/93‑1 (MQ=255)
                                                        aaaGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTAtt                                >  1:635144/1‑93 (MQ=255)
                                                         aaGAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTc                               >  1:3697288/1‑93 (MQ=255)
                                                           gAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGt                             <  1:3120246/93‑1 (MQ=255)
                                                           gAAACGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGt                             <  2:3504302/93‑1 (MQ=255)
                                                            aaaCGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTc                            <  1:4117590/93‑1 (MQ=255)
                                                            aaaCGGAGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTc                            <  1:3057978/93‑1 (MQ=255)
                                                             aaCGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCa                           >  2:4636405/1‑93 (MQ=255)
                                                             aaCGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCa                           <  1:4338891/93‑1 (MQ=255)
                                                              aCGGTGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCAt                          >  1:172350/1‑93 (MQ=255)
                                                                  tGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAaca                      <  2:515156/93‑1 (MQ=255)
                                                                  tGGAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAaca                      <  1:4572605/93‑1 (MQ=255)
                                                                   ggAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAacac                     >  2:861810/1‑93 (MQ=255)
                                                                   ggAAGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAacac                     <  1:1616892/93‑1 (MQ=255)
                                                                     aaGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACAccc                   <  2:1869398/93‑1 (MQ=255)
                                                                      aGCCTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCt                  <  2:2244155/93‑1 (MQ=255)
                                                                        ccTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTg                <  1:1359801/93‑1 (MQ=255)
                                                                         cTATCATATGAAAGGGGAGCGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTgg               >  1:4573930/1‑93 (MQ=255)
                                                                         cTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTgg               <  1:634383/93‑1 (MQ=255)
                                                                         cTATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTgg               >  1:1368079/1‑93 (MQ=255)
                                                                          tATCATATGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGc              <  1:2643015/93‑1 (MQ=255)
                                                                              atatGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACgg          >  1:481110/1‑93 (MQ=255)
                                                                               tatGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGa         >  2:1889795/1‑93 (MQ=255)
                                                                               tatGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGa         >  1:1340154/1‑93 (MQ=255)
                                                                                atGAAAGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGaa        <  2:157629/93‑1 (MQ=255)
                                                                                   aaaGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAAttt     >  1:3533659/1‑93 (MQ=255)
                                                                                   aaaGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAAttt     >  1:1988589/1‑93 (MQ=255)
                                                                                   aaaGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAAttt     <  2:8893/93‑1 (MQ=255)
                                                                                    aaGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTa    <  1:605028/93‑1 (MQ=255)
                                                                                     aGGGGAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTaa   >  1:1271630/1‑93 (MQ=255)
                                                                                      ggggAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTaaa  >  2:2805060/1‑93 (MQ=255)
                                                                                      ggggAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTaaa  <  2:2675315/93‑1 (MQ=255)
                                                                                      ggggAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGGATTCGTCATAACACCCTTGGCACGGAATTTaaa  >  1:3475956/1‑93 (MQ=255)
                                                                                         |                                                                                         
CTCCGGGAGCTGGTGATGAAATTCAGCTCACCGACGCAATTGATATGCTGATCGAAAAAGAAACGGTGGAAGCCTATCATATGAAAGGGAAGAGCCATGACTGCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAA  >  NC_000913/1291374‑1291552

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