Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,282,191 C→T D127D (GAC→GAT agaB → PTS system galactosamine‑specific IIB component 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,282,1910CT88.0% 889.3 / 72.4 299D127D (GAC→GATagaBPTS system galactosamine‑specific IIB component 1
Reads supporting (aligned to +/- strand):  ref base C (14/22);  new base T (126/137);  total (140/159)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 300 total aligned reads displayed.

AGGTGGTATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAATCCC  >  NC_000913/3282101‑3282280
                                                                                          |                                                                                          
aGGTGGTATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCa                                                                                          <  1:4133148/93‑1 (MQ=255)
aGGTGGTATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCa                                                                                          >  2:2126971/1‑93 (MQ=255)
  gtggtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAgg                                                                                        <  2:4924080/93‑1 (MQ=255)
  gtggtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAgg                                                                                        <  2:5101323/93‑1 (MQ=255)
   tggtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGa                                                                                       <  1:5269377/93‑1 (MQ=255)
   tggtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGa                                                                                       <  2:1078306/93‑1 (MQ=255)
   tggtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGa                                                                                       <  1:1990521/93‑1 (MQ=255)
     gtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGAtc                                                                                     >  1:1664644/1‑93 (MQ=255)
     gtATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGAtc                                                                                     >  1:3545736/1‑93 (MQ=255)
      tATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGAtct                                                                                    <  2:4539354/93‑1 (MQ=255)
       aTTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGAtctc                                                                                   >  1:5320579/1‑93 (MQ=255)
       aTTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGAtctc                                                                                   <  2:4191305/93‑1 (MQ=255)
        ttGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCa                                                                                  <  1:1265201/93‑1 (MQ=255)
        ttGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCa                                                                                  <  2:7296453/93‑1 (MQ=255)
          gACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACg                                                                                >  1:7098750/1‑93 (MQ=255)
             cTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGAc                                                                             >  2:3970565/1‑93 (MQ=255)
             cTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGAc                                                                             >  2:2501279/1‑93 (MQ=255)
             cTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGAc                                                                             >  1:204462/1‑93 (MQ=255)
              tGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACt                                                                            >  1:3400680/1‑93 (MQ=255)
              tGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACt                                                                            >  2:1803549/1‑93 (MQ=255)
                aaaGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTa                                                                          <  2:6919764/93‑1 (MQ=255)
                aaaGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTa                                                                          <  2:4743759/93‑1 (MQ=255)
                  aGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACg                                                                        <  1:3295546/93‑1 (MQ=255)
                  aGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACg                                                                        >  1:6229257/1‑93 (MQ=255)
                  aGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACg                                                                        >  1:4563454/1‑93 (MQ=255)
                   gATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGt                                                                       >  2:581674/1‑93 (MQ=255)
                    aTGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGtt                                                                      <  2:5575892/93‑1 (MQ=255)
                    aTGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGtt                                                                      <  2:3218366/93‑1 (MQ=255)
                     tGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGttt                                                                     <  1:1763782/93‑1 (MQ=255)
                      gTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGtttt                                                                    >  1:7498611/1‑93 (MQ=255)
                         aaCGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTttatt                                                                 <  1:836675/93‑1 (MQ=255)
                          aCGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTttatta                                                                >  1:7443478/1‑93 (MQ=255)
                          aCGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTttatta                                                                >  1:1160791/1‑93 (MQ=255)
                           cGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTaa                                                               >  2:157186/1‑93 (MQ=255)
                            gTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTaaa                                                              <  2:3856098/93‑1 (MQ=255)
                             tCGGCAATATGCATTTCTCGGAAGGGAAATAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAaac                                                             <  2:2255627/93‑1 (MQ=255)
                             tCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAaac                                                             >  2:343781/1‑93 (MQ=255)
                               ggCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAaacaa                                                           >  1:1086689/1‑93 (MQ=255)
                               ggCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAaacaa                                                           <  1:7240973/93‑1 (MQ=255)
                               ggCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAaacaa                                                           >  1:3727629/1‑93 (MQ=255)
                               ggCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAaacaa                                                           <  2:573199/93‑1 (MQ=255)
                                  aaTATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGt                                                        >  1:5034769/1‑93 (MQ=255)
                                  aaTATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGt                                                        >  2:6777049/1‑93 (MQ=255)
                                  aaTATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGt                                                        >  2:648988/1‑93 (MQ=255)
                                   atatGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTg                                                       <  1:2354081/93‑1 (MQ=255)
                                    tatGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTgg                                                      >  2:72325/1‑93 (MQ=255)
                                    tatGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTgg                                                      >  2:7039903/1‑93 (MQ=255)
                                    tatGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTgg                                                      <  2:5314589/93‑1 (MQ=255)
                                    tatGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTgg                                                      >  2:5969280/1‑93 (MQ=255)
                                    tatGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTgg                                                      >  2:331510/1‑93 (MQ=255)
                                     atGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGc                                                     >  2:167782/1‑93 (MQ=255)
                                         aTTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGa                                                >  2:3797225/1‑93 (MQ=255)
                                           ttCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAAt                                              >  1:3141948/1‑93 (MQ=255)
                                            tctcGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATg                                             >  2:5852137/1‑93 (MQ=255)
                                            tctcGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATg                                             >  1:6016518/1‑93 (MQ=255)
                                              tcGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGtt                                           >  2:2320839/1‑93 (MQ=255)
                                               cGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGttt                                          <  2:6607207/93‑1 (MQ=255)
                                               cGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGttt                                          >  1:1865282/1‑93 (MQ=255)
                                               cGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGttt                                          <  1:4953963/93‑1 (MQ=255)
                                                ggAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGtttt                                         >  1:2587371/1‑93 (MQ=255)
                                                 gAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttc                                        <  1:4613002/93‑1 (MQ=255)
                                                 gAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttc                                        >  1:6119156/1‑93 (MQ=255)
                                                 gAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttc                                        <  2:4348044/93‑1 (MQ=255)
                                                 gAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttc                                        <  1:717185/93‑1 (MQ=255)
                                                  aaGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttca                                       >  2:789068/1‑93 (MQ=255)
                                                      gAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttcattca                                   >  2:3578081/1‑93 (MQ=255)
                                                      gAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTttcattca                                   <  1:1525590/93‑1 (MQ=255)
                                                       aaaaaGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAg                                  <  2:522692/93‑1 (MQ=255)
                                                         aaaGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGa                                <  1:1827536/93‑1 (MQ=255)
                                                         aaaGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGa                                <  2:3087460/93‑1 (MQ=255)
                                                         aaaGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGa                                >  1:5548673/1‑93 (MQ=255)
                                                         aaaGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGa                                <  1:6680991/93‑1 (MQ=255)
                                                          aaGCAAATCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGAc                               >  1:1471611/1‑93 (MQ=255)
                                                              aaaTCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTcc                           <  2:2181322/93‑1 (MQ=255)
                                                               aaTCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTccc                          >  2:4102999/1‑93 (MQ=255)
                                                               aaTCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTccc                          <  2:4427793/93‑1 (MQ=255)
                                                               aaTCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTccc                          >  1:2951542/1‑93 (MQ=255)
                                                                aTCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCt                         <  1:4377138/93‑1 (MQ=255)
                                                                aTCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCt                         >  1:5917587/1‑93 (MQ=255)
                                                                 tCAGCAGTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTg                        <  1:23982/93‑1 (MQ=255)
                                                                     cagTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGa                    <  2:2023093/93‑1 (MQ=255)
                                                                      agTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTCGAATGTTTTCATTCAGGACGTCCCTGGCGa                    >  1:2368840/1‑93 (MQ=255)
                                                                      agTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGAt                   <  2:2440083/93‑1 (MQ=255)
                                                                       gTAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATc                  >  2:2793137/1‑93 (MQ=255)
                                                                        tAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCa                 <  2:1258215/93‑1 (MQ=255)
                                                                        tAAAGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCa                 <  1:5771832/93‑1 (MQ=255)
                                                                         aaaGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCaa                >  2:1366252/1‑93 (MQ=255)
                                                                          aaGTTTATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCaaa               >  1:6454239/1‑93 (MQ=255)
                                                                             tttATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAg            >  1:1426661/1‑93 (MQ=255)
                                                                             tttATGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAg            <  2:2978820/93‑1 (MQ=255)
                                                                                 tGTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACa        <  2:5001172/93‑1 (MQ=255)
                                                                                 tGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACa        <  2:6336379/93‑1 (MQ=255)
                                                                                  gTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACaa       >  1:404160/1‑93 (MQ=255)
                                                                                  gTCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACaa       >  1:2689817/1‑93 (MQ=255)
                                                                                   tCGATGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACaaa      <  1:4768968/93‑1 (MQ=255)
                                                                                    cGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAAt     <  1:6952105/93‑1 (MQ=255)
                                                                                    cGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAAt     <  1:6342639/93‑1 (MQ=255)
                                                                                       tGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAATccc  >  2:6630962/1‑93 (MQ=255)
                                                                                       tGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAATccc  <  2:2694495/93‑1 (MQ=255)
                                                                                       tGATCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAATccc  <  1:5579175/93‑1 (MQ=255)
                                                                                          |                                                                                          
AGGTGGTATTGACCTGAAAGATGTCAACGTCGGCAATATGCATTTCTCGGAAGGGAAAAAGCAAATCAGCAGTAAAGTTTATGTCGATGACCAGGATCTCACGGACTTACGTTTTATTAAACAACGTGGCGTGAATGTTTTCATTCAGGACGTCCCTGGCGATCAAAAAGAACAAATCCC  >  NC_000913/3282101‑3282280

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