Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,622,927 T→C F387L (TTT→CTT)  yjjJ → toxin YjjJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,622,9270TC99.7% 1529.7 / ‑4.9 356F387L (TTT→CTT) yjjJtoxin YjjJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (167/188);  minor base A (1/0);  total (168/188)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

Only 100 of 361 total aligned reads displayed.

AACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGC  >  NC_000913/4622829‑4623024
                                                                                                  |                                                                                                 
aaCCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGctt                                                                                                 <  2:8600026/101‑1 (MQ=255)
aaCCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGctt                                                                                                 >  1:8600768/1‑101 (MQ=255)
  ccGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGa                                                                                               <  2:8600383/101‑1 (MQ=255)
   cGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGAt                                                                                              >  2:8600772/1‑101 (MQ=255)
    gCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATc                                                                                             >  2:8600777/1‑101 (MQ=255)
    gCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATc                                                                                             >  1:8600778/1‑101 (MQ=255)
     ccATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCt                                                                                            >  1:8600782/1‑101 (MQ=255)
     ccATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCt                                                                                            <  1:8600260/101‑1 (MQ=255)
       aTTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTa                                                                                          >  2:8600788/1‑101 (MQ=255)
       aTTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTa                                                                                          <  1:8600461/101‑1 (MQ=255)
        tttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTaa                                                                                         >  1:8600790/1‑101 (MQ=255)
         ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAAc                                                                                        <  1:8600201/101‑1 (MQ=255)
         ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAAc                                                                                        <  1:8600288/101‑1 (MQ=255)
            cgcTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTc                                                                                     >  2:8600796/1‑101 (MQ=255)
              cTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAg                                                                                   >  2:8600803/1‑101 (MQ=255)
              cTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAg                                                                                   >  1:8600801/1‑101 (MQ=255)
              cTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAg                                                                                   <  1:8600433/101‑1 (MQ=255)
               tGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGt                                                                                  <  2:8600464/101‑1 (MQ=255)
                gACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTa                                                                                 <  1:8599786/101‑1 (MQ=255)
                gACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTa                                                                                 <  2:8600249/101‑1 (MQ=255)
                gACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTa                                                                                 <  2:8600071/101‑1 (MQ=255)
                gACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTa                                                                                 >  2:8600807/1‑101 (MQ=255)
                 aCGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTaa                                                                                <  1:8600347/101‑1 (MQ=255)
                  cGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTaaa                                                                               >  1:8600809/1‑101 (MQ=255)
                  cGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTaaa                                                                               <  1:8600103/101‑1 (MQ=255)
                   gCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTaaaa                                                                              >  1:8600815/1‑101 (MQ=255)
                    cccGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAg                                                                             <  2:8600092/101‑1 (MQ=255)
                      cGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAgcg                                                                           <  1:8600296/101‑1 (MQ=255)
                      cGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAgcg                                                                           <  1:8600151/101‑1 (MQ=255)
                      cGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAgcg                                                                           >  2:8600818/1‑101 (MQ=255)
                      cGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAgcg                                                                           <  2:8600324/101‑1 (MQ=255)
                       gTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGt                                                                                  <  1:8600222/93‑1 (MQ=255)
                       gTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTaaaa                                                                              <  1:8600035/97‑1 (MQ=255)
                        tCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCt                                                                         >  1:8600821/1‑101 (MQ=255)
                         cTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCtt                                                                        >  1:8600825/1‑101 (MQ=255)
                          tACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTg                                                                       <  2:8599974/101‑1 (MQ=255)
                           aCGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTgg                                                                      <  1:8599951/101‑1 (MQ=255)
                            cGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGt                                                                     >  2:8600829/1‑101 (MQ=255)
                             gACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGtt                                                                    >  1:8600830/1‑101 (MQ=255)
                               cATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTaa                                                                  <  1:8600323/101‑1 (MQ=255)
                                aTGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAAc                                                                 <  2:8600409/101‑1 (MQ=255)
                                 tgctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAgc                                                                            >  1:8600837/1‑89 (MQ=255)
                                 tgctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACg                                                                <  2:8600302/101‑1 (MQ=255)
                                 tgctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACg                                                                <  1:8600242/101‑1 (MQ=255)
                                  gctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACgg                                                               <  2:8600346/101‑1 (MQ=255)
                                  gctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACgg                                                               <  1:8599888/101‑1 (MQ=255)
                                  gctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACgg                                                               >  1:8600840/1‑101 (MQ=255)
                                   ctgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                              <  1:8600138/101‑1 (MQ=255)
                                    tgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCg                                                             <  2:8600367/101‑1 (MQ=255)
                                    tgcCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCg                                                             >  2:8600841/1‑101 (MQ=255)
                                      ccGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGAt                                                           <  2:8600206/101‑1 (MQ=255)
                                        gATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATcc                                                         >  2:8600847/1‑101 (MQ=255)
                                         aTGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATccc                                                        <  2:8600334/101‑1 (MQ=255)
                                          tGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCg                                                       <  2:8600392/101‑1 (MQ=255)
                                            gTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCt                                                     <  2:8600285/101‑1 (MQ=255)
                                             tCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTg                                                    >  1:8600851/1‑101 (MQ=255)
                                               tATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgc                                                  <  2:8600216/101‑1 (MQ=255)
                                               tATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgc                                                  >  2:8600854/1‑101 (MQ=255)
                                               tATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgc                                                  <  2:8600559/101‑1 (MQ=255)
                                                aTGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgcg                                                 >  1:8600858/1‑101 (MQ=255)
                                                 tGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgcgc                                                <  1:8600273/101‑1 (MQ=255)
                                                  gccccAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACg                                                                <  2:8600171/81‑1 (MQ=255)
                                                    aCCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagc                                             >  2:8600861/1‑101 (MQ=255)
                                                    aCCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagc                                             <  1:8600379/101‑1 (MQ=255)
                                                     ccAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagca                                            <  2:8600317/101‑1 (MQ=255)
                                                      cAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGcagcag                                           >  1:8600865/1‑101 (MQ=255)
                                                        aaCAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGtt                                         >  2:8600868/1‑101 (MQ=255)
                                                            gcgcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTgg                                     <  1:8600122/101‑1 (MQ=255)
                                                             cgcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGa                                                            >  1:8600873/1‑77 (MQ=255)
                                                             cgcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGc                                    <  1:8600558/101‑1 (MQ=255)
                                                              gcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCa                                   >  2:8600876/1‑101 (MQ=255)
                                                               cTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCaa                                  >  1:8600880/1‑101 (MQ=255)
                                                                 ggAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAAc                                >  1:8600883/1‑101 (MQ=255)
                                                                  gAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACg                               >  2:8600885/1‑101 (MQ=255)
                                                                  gAATGCTGCGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACg                               <  2:8600532/101‑1 (MQ=255)
                                                                     tgctgcGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTc                            <  1:8600158/101‑1 (MQ=255)
                                                                     tgctgcGTGATGCTGCCATTGAGGTGAAGATTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTc                            >  2:8600888/1‑101 (MQ=255)
                                                                      gctgcGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCg                           <  2:8600270/101‑1 (MQ=255)
                                                                       ctgcGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGc                          >  1:8600892/1‑101 (MQ=255)
                                                                        tgcGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGcc                         <  2:8600450/101‑1 (MQ=255)
                                                                         gcGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCa                                   <  2:8600142/90‑1 (MQ=255)
                                                                          cGTGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCag                       >  2:8600896/1‑101 (MQ=255)
                                                                            tGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCagag                     >  2:8600901/1‑101 (MQ=255)
                                                                            tGATGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCagag                     >  1:8600903/1‑101 (MQ=255)
                                                                              aTGCTGCCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCagag                     >  2:8600906/1‑99 (MQ=255)
                                                                                gctgcCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTaa                                                                  >  1:8600910/1‑52 (MQ=255)
                                                                                  tgcCATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCcgc               <  1:8600446/101‑1 (MQ=255)
                                                                                    ccATTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGt             <  2:8600358/101‑1 (MQ=255)
                                                                                      aTTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCcgcg              >  1:8600918/1‑98 (MQ=255)
                                                                                      aTTGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTAt           >  2:8600917/1‑101 (MQ=255)
                                                                                        tGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCa         <  1:8600494/101‑1 (MQ=255)
                                                                                        tGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCa         <  2:8600338/101‑1 (MQ=255)
                                                                                        tGAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCa         <  2:8600292/101‑1 (MQ=255)
                                                                                         gAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAg        <  1:8600183/101‑1 (MQ=255)
                                                                                         gAGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAg        >  2:8600923/1‑101 (MQ=255)
                                                                                          aGGTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGc       >  1:8600925/1‑101 (MQ=255)
                                                                                           ggTGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCg      >  2:8600931/1‑101 (MQ=255)
                                                                                             tGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAg    <  2:8600238/101‑1 (MQ=255)
                                                                                             tGAAGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAg    >  1:8600933/1‑101 (MQ=255)
                                                                                               aaGCTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGc  >  1:8600936/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
AACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGC  >  NC_000913/4622829‑4623024

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