Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,610 G→A R450R (CGC→CGT yaaJ ← inner membrane transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,6100GA85.7% 831.7 / 84.3 307R450R (CGC→CGTyaaJinner membrane transport protein
Reads supporting (aligned to +/- strand):  ref base G (27/17);  new base A (132/131);  total (159/148)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 315 total aligned reads displayed.

GATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAG  >  NC_000913/6518‑6701
                                                                                            |                                                                                           
  tAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGacg                                                                                           <  1:6992272/93‑1 (MQ=255)
    gTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCa                                                                                         >  1:6858365/1‑93 (MQ=255)
    gTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCa                                                                                         >  2:743750/1‑93 (MQ=255)
     ttGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCAc                                                                                        >  1:2281015/1‑93 (MQ=255)
      tgtgATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACg                                                                                       >  2:5954583/1‑93 (MQ=255)
        tgATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGcc                                                                                     <  2:1864509/93‑1 (MQ=255)
          aTTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGa                                                                                   <  2:1452062/93‑1 (MQ=255)
          aTTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGa                                                                                   <  2:127777/93‑1 (MQ=255)
          aTTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGa                                                                                   >  1:2089816/1‑93 (MQ=255)
           ttACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAg                                                                                  >  2:6674400/1‑93 (MQ=255)
            tACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGt                                                                                 <  1:4204855/93‑1 (MQ=255)
            tACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGt                                                                                 <  1:463126/93‑1 (MQ=255)
             aCTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGtt                                                                                <  2:3674083/93‑1 (MQ=255)
              cTCCTGCGAAACATTATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGttt                                                                               >  2:832647/1‑93 (MQ=255)
               tCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTa                                                                              <  2:45125/93‑1 (MQ=255)
               tCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTa                                                                              <  2:3137087/93‑1 (MQ=255)
               tCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTa                                                                              <  1:7119867/93‑1 (MQ=255)
                  tGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGc                                                                           >  1:2920590/1‑93 (MQ=255)
                  tGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGc                                                                           <  1:6731049/93‑1 (MQ=255)
                  tGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGc                                                                           <  2:1521991/93‑1 (MQ=255)
                   gCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCt                                                                          <  2:6360879/93‑1 (MQ=255)
                   gCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCt                                                                          >  2:7061137/1‑93 (MQ=255)
                     gAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTgg                                                                        <  2:518885/93‑1 (MQ=255)
                     gAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTgg                                                                        >  1:1776632/1‑93 (MQ=255)
                     gAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTgg                                                                        >  1:4434635/1‑93 (MQ=255)
                     gAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGAGCGCACGCCGAGCTTACGCTgg                                                                        >  1:760027/1‑93 (MQ=255)
                      aaaCATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGc                                                                       <  1:7320383/93‑1 (MQ=255)
                       aaCATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCg                                                                      >  2:4078541/1‑93 (MQ=255)
                          atcatcCTACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAg                                                                   >  2:2024571/1‑93 (MQ=255)
                           tcatcCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGa                                                                  >  2:3788624/1‑93 (MQ=255)
                           tcatcCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGa                                                                  >  1:253612/1‑93 (MQ=255)
                           tcaccCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGa                                                                  <  1:2712194/93‑1 (MQ=255)
                              tcCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATaa                                                               <  1:4368465/93‑1 (MQ=255)
                                 cACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAAtta                                                            >  1:1670487/1‑91 (MQ=255)
                                 cACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCa                                                            <  2:5356618/93‑1 (MQ=255)
                                 cACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCa                                                            <  2:3917286/93‑1 (MQ=255)
                                 cACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCa                                                            >  2:423646/1‑93 (MQ=255)
                                 cACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCa                                                            >  1:569379/1‑93 (MQ=255)
                                  aCGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCAc                                                           <  1:275999/93‑1 (MQ=255)
                                     cgTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTgg                                                        >  1:3028164/1‑93 (MQ=255)
                                      gTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGc                                                       >  2:7498004/1‑93 (MQ=255)
                                       tCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCa                                                      >  1:5627229/1‑93 (MQ=255)
                                       tCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCa                                                      >  2:5643762/1‑93 (MQ=255)
                                        ccGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGAACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCaa                                                     >  2:4685492/1‑93 (MQ=255)
                                        ccGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCaa                                                     <  2:2500676/93‑1 (MQ=255)
                                         cGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAAt                                                    >  2:4341921/1‑93 (MQ=255)
                                              aaaGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTAt                                               >  1:1979768/1‑93 (MQ=255)
                                              aaaGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTAt                                               <  1:360007/93‑1 (MQ=255)
                                               aaGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATg                                              >  2:7240065/1‑93 (MQ=255)
                                                 gCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGaa                                            <  2:886352/93‑1 (MQ=255)
                                                 gCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGaa                                            <  1:2994255/93‑1 (MQ=255)
                                                  cTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAAc                                           <  1:7492598/93‑1 (MQ=255)
                                                   tGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAAcc                                          >  1:618448/1‑93 (MQ=255)
                                                     gCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcac                                        >  2:4635151/1‑93 (MQ=255)
                                                     gCGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcac                                        <  1:6069772/93‑1 (MQ=255)
                                                      cGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcaca                                       <  1:556194/93‑1 (MQ=255)
                                                      cGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcaca                                       <  1:1412233/93‑1 (MQ=255)
                                                      cGACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcaca                                       >  1:1478606/1‑93 (MQ=255)
                                                       gACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAg                                      >  1:6268713/1‑93 (MQ=255)
                                                       gACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAg                                      <  2:1663660/93‑1 (MQ=255)
                                                       gACCGATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAg                                      >  1:3428162/1‑93 (MQ=255)
                                                           gATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCga                                  >  1:4102175/1‑93 (MQ=255)
                                                           gATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCga                                  >  1:51584/1‑93 (MQ=255)
                                                           gATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCga                                  >  2:2253526/1‑93 (MQ=255)
                                                           gATATCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCga                                  >  2:2760253/1‑93 (MQ=255)
                                                            atatCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCgag                                 >  2:1993674/1‑93 (MQ=255)
                                                            atatCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCgag                                 >  1:529542/1‑93 (MQ=255)
                                                            atatCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCgag                                 <  2:4896156/93‑1 (MQ=255)
                                                            atatCCGGAAAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCgag                                 <  1:3772966/93‑1 (MQ=255)
                                                             tatCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCgaga                                <  2:4314302/93‑1 (MQ=255)
                                                              atCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGagag                               >  1:4283049/1‑93 (MQ=255)
                                                              atCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGagag                               >  1:3306356/1‑93 (MQ=255)
                                                               tCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGc                              <  2:3394992/93‑1 (MQ=255)
                                                               tCCGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGc                              <  2:4402526/93‑1 (MQ=255)
                                                                ccGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCa                             <  2:1374276/93‑1 (MQ=255)
                                                                 cGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAg                            >  2:3427890/1‑93 (MQ=255)
                                                                 cGGATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAg                            <  1:5973840/93‑1 (MQ=255)
                                                                  ggATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGt                           >  2:6895695/1‑93 (MQ=255)
                                                                  ggATAACGCAATGGATCAAACACCGGACGCCCGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGACAGCAGt                           >  1:4746130/1‑93 (MQ=255)
                                                                  ggATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGt                           >  2:6266007/1‑93 (MQ=255)
                                                                  ggATAACGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGt                           <  1:1281847/93‑1 (MQ=255)
                                                                   gATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTa                          >  1:3217931/1‑93 (MQ=255)
                                                                   gATAACGCAAGGGATCAAACCCCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTa                          <  1:48376/93‑1 (MQ=255)
                                                                      aaCGCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTaaaa                       <  1:1155016/93‑1 (MQ=255)
                                                                         gCAATGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATgg                    <  1:195839/93‑1 (MQ=255)
                                                                          cAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGc                   <  1:6424997/93‑1 (MQ=255)
                                                                           aaTGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCg                  >  2:58272/1‑93 (MQ=255)
                                                                           aaTGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCg                  >  2:2342145/1‑93 (MQ=255)
                                                                           aaTGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCg                  <  1:3952732/93‑1 (MQ=255)
                                                                            aTGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCgg                 <  2:3268515/93‑1 (MQ=255)
                                                                            aTGGATCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCgg                 >  1:3500889/1‑93 (MQ=255)
                                                                               gATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCa              <  1:4940251/93‑1 (MQ=255)
                                                                                 tCAAACACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCaaa            >  2:6558539/1‑93 (MQ=255)
                                                                                    aaCACCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTg         >  1:2204836/1‑93 (MQ=255)
                                                                                     acacCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTgg        >  1:1001797/1‑93 (MQ=255)
                                                                                       acCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTa      <  2:5565006/93‑1 (MQ=255)
                                                                                       acCGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTa      >  1:72093/1‑93 (MQ=255)
                                                                                         cGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAt    >  1:1374547/1‑93 (MQ=255)
                                                                                         cGGACGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAt    <  1:6682487/93‑1 (MQ=255)
                                                                                           ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGCATGAACCACAGGCGAGAGCAGTAAAATGGCGGGCAAATTGGTAATAg  >  1:953058/1‑93 (MQ=255)
                                                                                            |                                                                                           
GATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAG  >  NC_000913/6518‑6701

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