Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,692,558 A→G L513P (CTG→CCG)  uidA ← beta‑D‑glucuronidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,692,5580AG99.0% 649.0 / ‑3.4 207L513P (CTG→CCG) uidAbeta‑D‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (112/93);  total (112/95)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 218 total aligned reads displayed.

ACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACC  >  NC_000913/1692468‑1692647
                                                                                          |                                                                                         
acgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGc                                                                                         >  2:9801110/1‑93 (MQ=255)
acgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGc                                                                                         >  2:511241/1‑93 (MQ=255)
  ggCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGccc                                                                                       <  2:6242025/93‑1 (MQ=255)
   gcgcTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCg                                                                                      >  2:7524643/1‑93 (MQ=255)
   gcgcTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCg                                                                                      >  1:201676/1‑93 (MQ=255)
   gcgcTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCg                                                                                      >  1:2743274/1‑93 (MQ=255)
      cTGACGCGCCCAAAGACGCGGTGATGCATATCCAGCCATGGACACTGATACTCTTCACTCCACATGTCGGTGAACAATGATTGCGGCCCGGCt                                                                                   >  2:4166167/1‑93 (MQ=255)
      cTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCt                                                                                   >  1:9219847/1‑93 (MQ=255)
        gACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTaa                                                                                 <  1:9801407/93‑1 (MQ=255)
           gcgATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGt                                                                              <  1:5294231/93‑1 (MQ=255)
            cgATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTa                                                                             <  2:4630389/93‑1 (MQ=255)
            cgATCAAAGACGCGCTGATACATATCCAGCCATGCCCGCTGATGCTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTa                                                                             <  1:2070833/93‑1 (MQ=255)
             gATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTAt                                                                            >  2:5883337/1‑93 (MQ=255)
             gATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTAt                                                                            >  2:211836/1‑93 (MQ=255)
              aTCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGGTAACGGATc                                                                           >  2:8429979/1‑93 (MQ=255)
              aTCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATc                                                                           >  2:3266241/1‑93 (MQ=255)
              aTCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATc                                                                           >  2:9056059/1‑93 (MQ=255)
                cAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCa                                                                         <  2:5643226/93‑1 (MQ=255)
                cAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCa                                                                         >  2:4095878/1‑93 (MQ=255)
                 aaaGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCAc                                                                        >  1:9970484/1‑93 (MQ=255)
                  aaGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACg                                                                       >  1:1049867/1‑93 (MQ=255)
                  aaGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACg                                                                       >  1:4843770/1‑93 (MQ=255)
                    gACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGcc                                                                     <  2:3892103/93‑1 (MQ=255)
                    gACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGcc                                                                     >  2:1616683/1‑93 (MQ=255)
                    gACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGcc                                                                     >  1:3228571/1‑93 (MQ=255)
                       gcgGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTa                                                                  <  2:4004403/93‑1 (MQ=255)
                         ggTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTAtt                                                                >  1:1283224/1‑93 (MQ=255)
                         ggTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTAtt                                                                >  2:289839/1‑93 (MQ=255)
                           tGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTCTTCg                                                              >  1:6228576/1‑93 (MQ=255)
                           tGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCg                                                              <  2:724388/93‑1 (MQ=255)
                           tGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCg                                                              >  1:9475998/1‑93 (MQ=255)
                           tGATACATATCCAGCCACGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCg                                                              >  1:3810471/1‑93 (MQ=255)
                            gATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCgg                                                             <  1:9058897/93‑1 (MQ=255)
                             aTACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGt                                                            <  2:2342102/93‑1 (MQ=255)
                              tACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGtg                                                           >  1:5725723/1‑93 (MQ=255)
                               aCATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGtga                                                          <  2:8614440/93‑1 (MQ=255)
                                cATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGtgat                                                         <  2:4860587/93‑1 (MQ=255)
                                 atatCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGtgatg                                                        >  2:1275673/1‑93 (MQ=255)
                                  tatCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGtgatga                                                       >  2:4094447/1‑93 (MQ=255)
                                  tatCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGtgatga                                                       >  1:2403332/1‑93 (MQ=255)
                                    tCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATa                                                     <  2:8847788/93‑1 (MQ=255)
                                     ccAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATaa                                                    >  2:5330063/1‑93 (MQ=255)
                                      cAGCCATGCACACTGATCCTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAAt                                                   <  1:1944/93‑1 (MQ=255)
                                        gCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCg                                                 >  2:5949176/1‑93 (MQ=255)
                                         ccATGCACACTGATACTCTTCACTCCTCATGGCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCgg                                                <  1:3568698/93‑1 (MQ=255)
                                          cATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGc                                               >  1:1872916/1‑93 (MQ=255)
                                          cATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGc                                               <  1:2818062/93‑1 (MQ=255)
                                          cATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGc                                               <  2:1473079/93‑1 (MQ=255)
                                            tGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTg                                             >  2:8348909/1‑93 (MQ=255)
                                            tGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTg                                             >  1:3026058/1‑93 (MQ=255)
                                             gCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGa                                            <  2:8142970/93‑1 (MQ=255)
                                             gCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGa                                            <  2:2697830/93‑1 (MQ=255)
                                              cacaCTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGAt                                           <  2:3499759/93‑1 (MQ=255)
                                                 acTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCa                                        <  2:5684051/93‑1 (MQ=255)
                                                 acTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCa                                        <  1:3631091/93‑1 (MQ=255)
                                                 acTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCa                                        <  2:295652/93‑1 (MQ=255)
                                                  cTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAg                                       >  1:7793796/1‑93 (MQ=255)
                                                   tGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGt                                      >  1:3822844/1‑93 (MQ=255)
                                                   tGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGt                                      >  1:6464471/1‑93 (MQ=255)
                                                    gATTCTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCGGATGCAGtt                                     >  2:9252138/1‑93 (MQ=255)
                                                    gATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAggt                                     >  2:775743/1‑91 (MQ=255)
                                                    gATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGtt                                     <  2:6563620/93‑1 (MQ=255)
                                                     aTACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGttt                                    <  1:6925189/93‑1 (MQ=255)
                                                     aTACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGttt                                    <  1:7309552/93‑1 (MQ=255)
                                                      tACTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtc                                   >  2:4706471/1‑93 (MQ=255)
                                                       aCTCTTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtct                                  >  1:239058/1‑93 (MQ=255)
                                                        ctctTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtctc                                 <  1:5609855/93‑1 (MQ=255)
                                                        ctctTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtctc                                 <  1:7431379/93‑1 (MQ=255)
                                                        cgctTCACTCCACATGTCGTTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTTGGTGATGATAATCGGCTGATGCAGTTtctc                                 <  1:1571088/91‑1 (MQ=255)
                                                         tctTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTcc                                <  2:538064/93‑1 (MQ=255)
                                                          ctTCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCt                               <  1:1732459/93‑1 (MQ=255)
                                                           ttCACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTg                              >  2:2515418/1‑93 (MQ=255)
                                                             cACTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgcc                            <  1:3955991/93‑1 (MQ=255)
                                                              ccTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgcca                           <  1:8373905/92‑1 (MQ=255)
                                                              aCTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCAGGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgcca                           <  2:4529200/93‑1 (MQ=255)
                                                              aCTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgcca                           <  1:5024223/93‑1 (MQ=255)
                                                              aCTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgcca                           >  1:6846608/1‑93 (MQ=255)
                                                               cTCCACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccag                          <  1:6363155/93‑1 (MQ=255)
                                                                 ccACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggc                        >  2:6253200/1‑93 (MQ=255)
                                                                 ccACATGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggc                        >  2:2293320/1‑93 (MQ=255)
                                                                  cacaTGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggcc                       >  2:1528665/1‑93 (MQ=255)
                                                                   acaTGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggcca                      >  2:5192893/1‑93 (MQ=255)
                                                                   acaTGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggcca                      <  2:9857507/93‑1 (MQ=255)
                                                                    caTGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggccag                     >  1:8582864/1‑93 (MQ=255)
                                                                     aTGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGa                    >  2:4836824/1‑93 (MQ=255)
                                                                      tGTCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGaa                   <  2:9719760/93‑1 (MQ=255)
                                                                        tCGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGt                 <  2:9833466/93‑1 (MQ=255)
                                                                         cGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGtt                <  2:3624054/93‑1 (MQ=255)
                                                                         cGGTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGtt                >  2:5521221/1‑93 (MQ=255)
                                                                          ggTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTc               >  1:7632904/1‑93 (MQ=255)
                                                                          ggTGTACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTc               <  2:9608027/93‑1 (MQ=255)
                                                                          ggTGTACATTGAGTGCGGCCCGAATAACGTATCCACGCCGTCTTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTc               <  1:8938720/93‑1 (MQ=255)
                                                                             gtACATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCttt            >  1:4117232/1‑93 (MQ=255)
                                                                                cATTGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTc         >  2:8274170/1‑93 (MQ=255)
                                                                                  ttGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCa       >  1:376738/1‑93 (MQ=255)
                                                                                  ttGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCa       <  1:2922409/93‑1 (MQ=255)
                                                                                  ttGAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCa       <  1:9319929/93‑1 (MQ=255)
                                                                                  ttGAGCGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCa       <  1:7042817/93‑1 (MQ=255)
                                                                                    gAGTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGt     <  1:6025153/93‑1 (MQ=255)
                                                                                       tGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGAGGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTAcc  <  1:5801247/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACC  >  NC_000913/1692468‑1692647

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