Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,437,025 T→C Q595R (CAG→CGG)  ydbK ← putative oxidoreductase Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,437,0250TC86.4% 589.0 / 79.8 243Q595R (CAG→CGG) ydbKputative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (19/14);  new base C (105/105);  total (124/119)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

CTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGC  >  NC_000913/1436938‑1437116
                                                                                       |                                                                                           
 tgtgCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCtt                                                                                       >  1:6396387/1‑93 (MQ=255)
  gtgCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCttt                                                                                      <  2:6942094/93‑1 (MQ=255)
    gCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGc                                                                                    <  1:3951181/93‑1 (MQ=255)
      ggATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTa                                                                                  <  1:6274827/93‑1 (MQ=255)
      ggATTTACCGGTTGCAACGGAACTTCTTCTACGGACTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTa                                                                                  >  1:4621409/1‑93 (MQ=255)
       gATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTAc                                                                                 >  1:5251723/1‑93 (MQ=255)
        aTTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACt                                                                                >  1:4397853/1‑93 (MQ=255)
        aTTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACt                                                                                <  1:6636462/93‑1 (MQ=255)
         tttACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTg                                                                               >  1:1284794/1‑93 (MQ=255)
         tttACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTg                                                                               >  1:7215791/1‑93 (MQ=255)
         tttACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTg                                                                               <  1:4564337/93‑1 (MQ=255)
          ttACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGt                                                                              <  1:276951/93‑1 (MQ=255)
           tACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTa                                                                             >  2:3238954/1‑93 (MQ=255)
           tACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTa                                                                             >  1:1510252/1‑93 (MQ=255)
             ccGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAAc                                                                           >  1:1090034/1‑93 (MQ=255)
              cGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACt                                                                          >  2:4818071/1‑93 (MQ=255)
               ggTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACtt                                                                         >  2:7478060/1‑93 (MQ=255)
                 ttGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACtttt                                                                       >  1:2127812/1‑93 (MQ=255)
                 ttGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACtttt                                                                       >  1:5108652/1‑93 (MQ=255)
                  tGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTg                                                                      <  2:2258548/93‑1 (MQ=255)
                  tGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTg                                                                      <  1:3137143/93‑1 (MQ=255)
                      aCGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCaa                                                                  <  2:4476419/93‑1 (MQ=255)
                       cGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAAt                                                                 >  2:3390413/1‑93 (MQ=255)
                       cGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAAt                                                                 >  2:4159917/1‑93 (MQ=255)
                       cGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAAt                                                                 <  2:6792178/93‑1 (MQ=255)
                        ggAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATc                                                                >  1:3699166/1‑93 (MQ=255)
                         gAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCg                                                               >  2:561990/1‑93 (MQ=255)
                         gAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCg                                                               <  2:1569298/93‑1 (MQ=255)
                          aaCTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGc                                                              <  1:3824941/93‑1 (MQ=255)
                           aCTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCa                                                             >  2:3704798/1‑93 (MQ=255)
                               cttctACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCt                                                         <  2:2958644/93‑1 (MQ=255)
                               cttctACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCt                                                         >  2:5495938/1‑93 (MQ=255)
                               cttctACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCt                                                         <  1:3788072/93‑1 (MQ=255)
                               cttctACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCt                                                         <  1:5559601/93‑1 (MQ=255)
                               cttctACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCt                                                         >  1:2051488/1‑93 (MQ=255)
                                   tACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCaa                                                     <  2:5041306/93‑1 (MQ=255)
                                   tACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCaa                                                     <  1:2365076/93‑1 (MQ=255)
                                   tACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCaa                                                     <  1:4054469/93‑1 (MQ=255)
                                    aCGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAAt                                                    >  1:6042945/1‑93 (MQ=255)
                                    aCGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAAt                                                    <  1:5896372/93‑1 (MQ=255)
                                     cGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAAtt                                                   >  2:5341300/1‑93 (MQ=255)
                                      ggATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGTGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTc                                                  >  2:1098188/1‑93 (MQ=255)
                                       gATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCt                                                 >  1:7100632/1‑93 (MQ=255)
                                       gATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCt                                                 <  2:3521088/93‑1 (MQ=255)
                                       gATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCt                                                 >  1:1794066/1‑93 (MQ=255)
                                       gATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCt                                                 <  2:3906435/93‑1 (MQ=255)
                                           cACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGa                                             <  1:679612/93‑1 (MQ=255)
                                           cACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGa                                             >  2:6595065/1‑93 (MQ=255)
                                            aCGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAg                                            >  2:3285124/1‑93 (MQ=255)
                                             cgcgCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAgg                                           <  2:3158384/93‑1 (MQ=255)
                                              gcgccagcgccaGAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAggg                                          <  2:5815235/93‑1 (MQ=255)
                                               cgccagcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGgc                                         >  2:6389003/1‑93 (MQ=255)
                                                 ccagcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGgcgc                                       >  1:429232/1‑93 (MQ=255)
                                                    gcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTAt                                    <  2:7043166/93‑1 (MQ=255)
                                                    gcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTAt                                    >  2:2335249/1‑93 (MQ=255)
                                                       ccagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGc                                 >  2:6199191/1‑93 (MQ=255)
                                                       ccagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGc                                 >  2:7100337/1‑93 (MQ=255)
                                                        cagAGCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGcc                                <  2:5388368/93‑1 (MQ=255)
                                                         agagCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCa                               <  1:2075817/93‑1 (MQ=255)
                                                         agagCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCa                               <  2:1306825/93‑1 (MQ=255)
                                                         agagCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCa                               >  2:4610623/1‑93 (MQ=255)
                                                          gagCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAg                              <  1:742408/93‑1 (MQ=255)
                                                          gagCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAg                              <  2:1907257/93‑1 (MQ=255)
                                                           agCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAgg                             <  1:3463558/93‑1 (MQ=255)
                                                           agCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAgg                             >  2:5977900/1‑93 (MQ=255)
                                                            gCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGc                            <  2:4580262/93‑1 (MQ=255)
                                                            gCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGc                            >  1:3528899/1‑93 (MQ=255)
                                                            gCCTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGc                            <  2:4354998/93‑1 (MQ=255)
                                                              cTGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAg                          >  1:491713/1‑93 (MQ=255)
                                                               tGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGa                         >  1:2951587/1‑93 (MQ=255)
                                                               tGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGa                         >  2:881458/1‑93 (MQ=255)
                                                               tGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGa                         >  1:336885/1‑93 (MQ=255)
                                                               tGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGa                         >  1:1941484/1‑93 (MQ=255)
                                                               tGCCAGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGa                         >  2:1351779/1‑93 (MQ=255)
                                                                   aGTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAAttt                     >  1:77984/1‑93 (MQ=255)
                                                                    gTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTg                    >  2:4934165/1‑93 (MQ=255)
                                                                    gTTGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTg                    >  1:615220/1‑93 (MQ=255)
                                                                     ttGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGc                   >  2:702991/1‑93 (MQ=255)
                                                                     ttGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGc                   <  2:4975167/93‑1 (MQ=255)
                                                                     ttGCGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGc                   <  2:1178122/93‑1 (MQ=255)
                                                                      tGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCg                  <  2:39619/93‑1 (MQ=255)
                                                                        cGTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTc                <  2:7505685/93‑1 (MQ=255)
                                                                         gTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCa               <  1:343028/93‑1 (MQ=255)
                                                                         gTTCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCa               <  1:701947/93‑1 (MQ=255)
                                                                          ttCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAg              >  1:4481130/1‑93 (MQ=255)
                                                                           tCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGa             <  2:7211586/93‑1 (MQ=255)
                                                                           tCCACCAGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGa             <  1:4269595/93‑1 (MQ=255)
                                                                            ccaccaGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGAt            >  2:6643537/1‑93 (MQ=255)
                                                                            ccaccaGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGAt            <  2:259309/93‑1 (MQ=255)
                                                                             caccaGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATg           <  2:166963/93‑1 (MQ=255)
                                                                             caccaGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATg           <  1:3055440/93‑1 (MQ=255)
                                                                              accaGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATgg          <  2:528326/93‑1 (MQ=255)
                                                                               ccaGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGa         <  1:3903851/93‑1 (MQ=255)
                                                                                caGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaa        >  1:654189/1‑93 (MQ=255)
                                                                                 aGATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaa       >  1:5972736/1‑93 (MQ=255)
                                                                                  gATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaa      >  1:1434244/1‑93 (MQ=255)
                                                                                  gATCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaa      >  2:2982895/1‑93 (MQ=255)
                                                                                   aTCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaaa     <  2:3769758/93‑1 (MQ=255)
                                                                                   aTCCCGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaaa     <  2:4717554/93‑1 (MQ=255)
                                                                                   aTCCCGGCCTTTGCTACAGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaaa     <  1:760790/93‑1 (MQ=255)
                                                                                       |                                                                                           
CTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGC  >  NC_000913/1436938‑1437116

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