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,8330TC94.4% 1295.6 / 45.2 448I460T (ATC→ACC) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (14/11);  new base C (235/188);  total (249/199)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 471 total aligned reads displayed.

TGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGCAATTGGG  >  NC_000913/3287743‑3287925
                                                                                          |                                                                                            
tgtgGATGGATATCAACAGAGCTGACTCCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATAccc                                                                                            <  1:657932/93‑1 (MQ=255)
  tCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCg                                                                                          <  2:1459677/93‑1 (MQ=255)
  tCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCg                                                                                          >  2:1831304/1‑93 (MQ=255)
    gATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGa                                                                                        >  2:937690/1‑93 (MQ=255)
     aTGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCACTGGCGTACGCAGTACATCACAGATACACCGGaa                                                                                       >  1:9712921/1‑93 (MQ=255)
     aTGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGaa                                                                                       >  1:5137426/1‑93 (MQ=255)
     aTGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGaa                                                                                       >  1:7212248/1‑93 (MQ=255)
      tGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGaaa                                                                                      <  1:241321/93‑1 (MQ=255)
      tGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGACTATCAATAGCGTACGCAGTACATCAAAGATACACCAGaaa                                                                                      >  2:4247382/1‑93 (MQ=255)
       ggATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAAc                                                                                     >  2:7173289/1‑93 (MQ=255)
        gATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAAcc                                                                                    >  1:1285288/1‑93 (MQ=255)
        gATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAAcc                                                                                    >  1:1931415/1‑93 (MQ=255)
         atatCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCa                                                                                   >  2:9707906/1‑93 (MQ=255)
           atCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAac                                                                                 >  1:7464618/1‑93 (MQ=255)
           atCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAac                                                                                 <  1:8595915/93‑1 (MQ=255)
            tCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAaca                                                                                >  1:9161020/1‑93 (MQ=255)
              aaCAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACAcc                                                                              >  2:9487553/1‑93 (MQ=255)
                  gagCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAata                                                                          <  1:6272908/93‑1 (MQ=255)
                   agCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAatat                                                                         >  1:638948/1‑93 (MQ=255)
                   agCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAatat                                                                         <  1:5025033/93‑1 (MQ=255)
                   agCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAatat                                                                         <  1:5782822/93‑1 (MQ=255)
                    gCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATc                                                                        >  2:6077115/1‑93 (MQ=255)
                     cTGACACCCAATTCGATAACCAGCACTCTTTTCATTGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCg                                                                       >  1:8942444/1‑93 (MQ=255)
                       gACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCt                                                                     >  1:1171260/1‑93 (MQ=255)
                        acacCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTg                                                                    >  2:4096961/1‑93 (MQ=255)
                        acacCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTg                                                                    <  1:4339967/93‑1 (MQ=255)
                         cacCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGt                                                                   <  1:1510468/93‑1 (MQ=255)
                          acCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTc                                                                  >  2:8232855/1‑93 (MQ=255)
                          acCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTc                                                                  >  1:8194100/1‑93 (MQ=255)
                           cccAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCa                                                                 <  1:5951278/93‑1 (MQ=255)
                           cccAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCa                                                                 >  1:8047656/1‑93 (MQ=255)
                            ccAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAg                                                                <  1:2614801/93‑1 (MQ=255)
                               aTTCGATAACCAGCACTCTTTTCATGGCTATCAAAGGCGTACGCAGTGCATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGcta                                                             >  1:5732339/1‑93 (MQ=255)
                                ttCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGctat                                                            >  2:9566269/1‑92 (MQ=255)
                                atCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGctac                                                            <  1:1117281/92‑1 (MQ=255)
                                  cGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGctacta                                                          <  2:3341207/93‑1 (MQ=255)
                                   gATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTAt                                                         >  2:3182430/1‑93 (MQ=255)
                                   gATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTAt                                                         >  2:4454501/1‑93 (MQ=255)
                                   gATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTAt                                                         <  1:1703898/93‑1 (MQ=255)
                                     tAACCAGCTCTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCg                                                       <  1:4079140/93‑1 (MQ=255)
                                     tAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCg                                                       >  2:7865339/1‑93 (MQ=255)
                                     tAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCg                                                       >  2:9861492/1‑93 (MQ=255)
                                        ccAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCta                                                    <  2:3523001/93‑1 (MQ=255)
                                        ccAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCta                                                    <  1:3184154/93‑1 (MQ=255)
                                          aGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtata                                                  >  1:8243087/1‑93 (MQ=255)
                                          aGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtata                                                  <  2:927086/93‑1 (MQ=255)
                                           gcgctctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATAc                                                 >  2:6629389/4‑93 (MQ=255)
                                             aCTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCa                                               >  2:186416/1‑93 (MQ=255)
                                              ctctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCaa                                              >  1:7313576/1‑93 (MQ=255)
                                               tctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAAc                                             >  2:3797104/1‑93 (MQ=255)
                                               tctTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAAc                                             >  2:9103517/1‑93 (MQ=255)
                                                 ttttCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGACACCAACACCAATATCGCTGTAAGCTACTATCGCTATACCAACGa                                           >  1:2954271/1‑93 (MQ=255)
                                                  tttCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGAt                                          <  1:6868950/93‑1 (MQ=255)
                                                   ttCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCCATACCAACGATg                                         >  1:7688176/1‑93 (MQ=255)
                                                    tCATGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATgg                                        >  2:2300687/1‑93 (MQ=255)
                                                      aTGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCt                                      >  2:4312326/1‑93 (MQ=255)
                                                       tGGCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTa                                     <  1:791836/93‑1 (MQ=255)
                                                         gCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtt                                   >  1:3810965/1‑93 (MQ=255)
                                                         gCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtt                                   >  1:2137934/1‑93 (MQ=255)
                                                         gCTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtt                                   <  1:9993606/93‑1 (MQ=255)
                                                          cTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAttt                                  >  1:9434309/1‑93 (MQ=255)
                                                          cTATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAttt                                  <  1:482059/93‑1 (MQ=255)
                                                           tATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtttt                                 >  2:532682/1‑93 (MQ=255)
                                                           tATCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAtttt                                 >  2:1206199/1‑93 (MQ=255)
                                                            aTCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTa                                >  2:7002673/1‑93 (MQ=255)
                                                            aTCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTa                                >  2:6279771/1‑93 (MQ=255)
                                                             tCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAg                               <  1:2546120/93‑1 (MQ=255)
                                                             tCAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAg                               >  2:3684697/1‑93 (MQ=255)
                                                              cAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGt                              <  2:1718493/93‑1 (MQ=255)
                                                              cAATGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGt                              >  2:2432569/1‑93 (MQ=255)
                                                               aaTGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtt                             >  2:8588909/1‑93 (MQ=255)
                                                               aaTGGCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtt                             <  2:6425588/93‑1 (MQ=255)
                                                                   gCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTaa                         <  1:4674148/93‑1 (MQ=255)
                                                                   gCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTaa                         <  2:2996291/93‑1 (MQ=255)
                                                                   gCGTACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTaa                         <  1:3683616/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      >  1:8350110/1‑93 (MQ=255)
                                                                      tACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      <  1:6187892/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      >  2:8399360/1‑93 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      >  2:6697989/1‑93 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      <  2:7682031/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      >  1:3338217/1‑93 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      <  1:1018906/93‑1 (MQ=255)
                                                                      tACGCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGa                      >  1:5635028/1‑93 (MQ=255)
                                                                         gCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGc                   >  1:9607765/1‑93 (MQ=255)
                                                                         gCAGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGc                   >  1:6508111/1‑93 (MQ=255)
                                                                           aGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATCCCAACGCTGGCTATTTTCGTTTTAATGaaccac                 >  1:1344015/1‑89 (MQ=255)
                                                                           aGTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCa                 >  2:2556789/1‑93 (MQ=255)
                                                                            gTACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCaa                <  2:8681000/93‑1 (MQ=255)
                                                                             tACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAAt               >  2:2126055/1‑93 (MQ=255)
                                                                             tACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAAt               <  2:4869951/93‑1 (MQ=255)
                                                                             tACATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAAt               <  2:281233/93‑1 (MQ=255)
                                                                               cATCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAc             >  1:2859350/1‑93 (MQ=255)
                                                                                aTCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAcc            >  1:5288302/1‑93 (MQ=255)
                                                                                aTCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAcc            >  2:2675688/1‑93 (MQ=255)
                                                                                aTCAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATAcc            >  2:7500351/1‑93 (MQ=255)
                                                                                  cAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCg          >  1:3009679/1‑93 (MQ=255)
                                                                                  cAAAGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTCTACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCg          >  1:3572028/1‑93 (MQ=255)
                                                                                   aaaGATACCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGc         >  2:4711243/1‑93 (MQ=255)
                                                                                     aGATACCCCGGAAACCAACACCAGTATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGCaa       <  1:6226127/93‑1 (MQ=255)
                                                                                      gATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGCAAt      >  2:586108/1‑93 (MQ=255)
                                                                                          |                                                                                            
TGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAATGAAGCCAATACCCGCAATTGGG  >  NC_000913/3287743‑3287925

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