Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,287,833 T→C I460T (ATC→ACC)  yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,287,8330TC99.6% 717.5 / NA 223I460T (ATC→ACC) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (106/115);  total (108/115)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

TGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGCAATT  >  NC_000913/3287743‑3287922
                                                                                          |                                                                                         
tGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATAccc                                                                                         <  2:269000/93‑1 (MQ=255)
 gTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATAcccc                                                                                        <  2:1703185/93‑1 (MQ=255)
   cGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCgg                                                                                      >  2:4072589/1‑93 (MQ=255)
    gATGGATAGCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGa                                                                                     >  2:6125639/1‑93 (MQ=255)
     aTGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGaa                                                                                    <  2:3961143/93‑1 (MQ=255)
      tGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGaaa                                                                                   >  1:4484364/1‑93 (MQ=255)
       ggATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAAc                                                                                  >  1:3051823/1‑93 (MQ=255)
         atatCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCa                                                                                >  2:186877/1‑93 (MQ=255)
         atatCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCa                                                                                >  2:4192808/1‑93 (MQ=255)
          tatCAACAGAGTTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCaa                                                                               <  2:3708510/93‑1 (MQ=255)
          tatCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCaa                                                                               >  1:6039072/1‑93 (MQ=255)
            tCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAaca                                                                             >  2:3054676/1‑93 (MQ=255)
               aCAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCa                                                                          >  1:4571879/1‑93 (MQ=255)
               aCAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCa                                                                          >  1:5886489/1‑93 (MQ=255)
                 agagCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAat                                                                        >  1:5028546/1‑93 (MQ=255)
                   accTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAatat                                                                      <  1:4653939/91‑1 (MQ=255)
                    gCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATc                                                                     <  1:6216636/93‑1 (MQ=255)
                    gCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATc                                                                     <  1:429068/93‑1 (MQ=255)
                      tGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGc                                                                   <  2:3302692/93‑1 (MQ=255)
                         cacCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGt                                                                >  2:4145472/1‑93 (MQ=255)
                          acCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGCCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTc                                                               <  2:1928863/93‑1 (MQ=255)
                           cccAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCa                                                              <  1:3901779/93‑1 (MQ=255)
                             cAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGTTGTCAGc                                                            <  2:5984691/93‑1 (MQ=255)
                             cAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGc                                                            <  2:342253/93‑1 (MQ=255)
                             cAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGc                                                            <  2:5328669/93‑1 (MQ=255)
                                  cGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGctacta                                                       <  1:1083494/93‑1 (MQ=255)
                                  cGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGctacta                                                       <  2:1451507/93‑1 (MQ=255)
                                  cGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGctacta                                                       <  1:3661967/93‑1 (MQ=255)
                                   gATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTAt                                                      <  1:5865867/93‑1 (MQ=255)
                                   gATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTAt                                                      >  1:5291158/1‑93 (MQ=255)
                                    aTAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATc                                                     <  2:1464881/93‑1 (MQ=255)
                                     tAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCg                                                    >  1:4210440/1‑93 (MQ=255)
                                        ccAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCta                                                 <  1:5374217/93‑1 (MQ=255)
                                        ccAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCta                                                 >  1:5575448/1‑93 (MQ=255)
                                        ccAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCta                                                 <  1:3372637/93‑1 (MQ=255)
                                        ccAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCta                                                 <  2:5485345/93‑1 (MQ=255)
                                         cAGCACTCTTTTCATGGCTATTAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtat                                                >  1:3993541/1‑93 (MQ=255)
                                         cAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtat                                                <  2:263400/93‑1 (MQ=255)
                                         cAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtat                                                >  1:2619179/1‑93 (MQ=255)
                                         cAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtat                                                >  1:879964/1‑93 (MQ=255)
                                          atcgctctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtata                                               <  1:1514222/89‑1 (MQ=255)
                                          aGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtatc                                               >  2:1377004/1‑92 (MQ=255)
                                            cACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATAcc                                             >  1:244337/1‑93 (MQ=255)
                                             aCTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCa                                            >  2:1685907/1‑93 (MQ=255)
                                             aCTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCa                                            >  1:1491025/1‑93 (MQ=255)
                                              ctctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCaa                                           >  1:811563/1‑93 (MQ=255)
                                              ctctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCaa                                           <  2:864838/93‑1 (MQ=255)
                                              ctctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCaa                                           >  1:436863/1‑93 (MQ=255)
                                                 ttttCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGa                                        >  2:3458023/1‑93 (MQ=255)
                                                  tttCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGAt                                       <  2:58500/93‑1 (MQ=255)
                                                     cATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGc                                    >  2:526364/1‑93 (MQ=255)
                                                     cATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGc                                    >  2:2005952/1‑93 (MQ=255)
                                                      aTGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCt                                   <  2:5441129/93‑1 (MQ=255)
                                                      aTGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCt                                   >  1:2822829/1‑93 (MQ=255)
                                                      aTGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCt                                   >  2:1345781/1‑93 (MQ=255)
                                                       tGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTa                                  <  1:2283514/93‑1 (MQ=255)
                                                        ggCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAt                                 >  2:5805528/1‑93 (MQ=255)
                                                         gCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtt                                >  2:2399215/1‑93 (MQ=255)
                                                         gCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtt                                >  1:3510694/1‑93 (MQ=255)
                                                           tATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtttt                              >  1:5400110/1‑93 (MQ=255)
                                                           tATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtttt                              >  1:3834369/1‑93 (MQ=255)
                                                            aTCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTa                             >  1:4818219/1‑93 (MQ=255)
                                                            aTCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTa                             <  1:1539664/93‑1 (MQ=255)
                                                             tCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAg                            <  2:5064151/93‑1 (MQ=255)
                                                               aaTGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtt                          <  2:151438/93‑1 (MQ=255)
                                                               aaTGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtt                          <  1:1717222/93‑1 (MQ=255)
                                                                aTGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGttt                         <  2:431867/93‑1 (MQ=255)
                                                                aTGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGttt                         >  1:1926638/1‑93 (MQ=255)
                                                                 tGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtttt                        >  2:5753821/1‑93 (MQ=255)
                                                                 tGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtttt                        >  2:720173/1‑93 (MQ=255)
                                                                   gCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTaa                      <  1:4621491/93‑1 (MQ=255)
                                                                    cGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAAt                     <  1:3154523/93‑1 (MQ=255)
                                                                    cGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAAt                     <  2:56488/93‑1 (MQ=255)
                                                                    cGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAAt                     >  1:2150425/1‑93 (MQ=255)
                                                                     gTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATg                    <  1:5960049/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                   <  1:1335667/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                   <  2:3970326/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                   <  1:2079234/93‑1 (MQ=255)
                                                                      tACGAAGTACATCAAAGACACCCAGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                   <  2:2479888/93‑1 (MQ=255)
                                                                       aCGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGaa                  >  1:1821187/1‑93 (MQ=255)
                                                                       aCGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGaa                  <  2:4706289/93‑1 (MQ=255)
                                                                       aCGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGaa                  <  2:2466164/93‑1 (MQ=255)
                                                                         gCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGc                >  2:4598563/1‑93 (MQ=255)
                                                                          cAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGcc               <  2:1988607/93‑1 (MQ=255)
                                                                          cAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGcc               >  2:3534519/1‑93 (MQ=255)
                                                                          cAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGcc               >  2:4793537/1‑93 (MQ=255)
                                                                          cAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGcc               >  1:6166659/1‑93 (MQ=255)
                                                                          cAGTACATAAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGcc               >  2:6157216/1‑93 (MQ=255)
                                                                           aGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCa              >  1:5726154/1‑93 (MQ=255)
                                                                           aGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCa              >  2:1066458/1‑93 (MQ=255)
                                                                            gcaCATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCaa             >  2:2766231/3‑93 (MQ=255)
                                                                            gTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCaa             >  1:959251/1‑93 (MQ=255)
                                                                              aCATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATa           <  1:1136840/93‑1 (MQ=255)
                                                                              aCATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATa           >  2:382949/1‑93 (MQ=255)
                                                                              aCATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATa           <  1:754967/93‑1 (MQ=255)
                                                                               cATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAc          >  1:5654533/1‑93 (MQ=255)
                                                                                aTCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAcc         <  2:1279604/93‑1 (MQ=255)
                                                                                aTCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAcc         <  1:2201421/93‑1 (MQ=255)
                                                                                 tCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAccc        >  2:4534203/1‑93 (MQ=255)
                                                                                 tCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAccc        >  2:4256484/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGCAATT  >  NC_000913/3287743‑3287922

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