Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,099,661 (A)5→6 intergenic (‑19/+164) yggN ← / ← b2959 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,099,6561.A95.9% 371.6 / 8.6 122intergenic (‑14/+169)yggN/b2959hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base A (53/63);  total (54/68)
Fisher's exact test for biased strand distribution p-value = 6.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 155 total aligned reads displayed.

TCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCG  >  NC_000913/3099581‑3099726
                                                                            |                                                                         
tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAg                                                                            <  1:1548583/76‑1 (MQ=255)
tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAg                                                                            <  1:4620735/76‑1 (MQ=255)
 cGGCGTGAGCGGTCATTGCCGTCACTGAACGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                          <  1:3038599/76‑1 (MQ=255)
 cGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                          <  2:3682691/76‑1 (MQ=255)
 cGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                          >  1:3391577/1‑76 (MQ=255)
  ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                         >  2:1143465/1‑76 (MQ=255)
  ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                         >  2:1105180/1‑76 (MQ=255)
  ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                         >  1:164163/1‑76 (MQ=255)
  ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                         <  1:2005030/76‑1 (MQ=255)
  ggCGTGAGCGGTCATTCCCGTCACTGAAAGTGCTGCCGCGAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                         <  2:3460143/76‑1 (MQ=255)
    cgggAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaa                                                                      <  1:4746925/73‑1 (MQ=255)
     gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                    <  2:917949/76‑1 (MQ=255)
     gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                    >  2:2384312/1‑76 (MQ=255)
     gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                    >  1:585220/1‑76 (MQ=255)
     gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                    <  1:2690731/76‑1 (MQ=255)
       gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaat                                                                  <  1:2008712/76‑3 (MQ=255)
       gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTGGCGCATCATAGTCTTCCCTCAAGaaaaaat                                                                  <  2:2422537/76‑3 (MQ=255)
        aGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatc                                                                 >  1:153501/1‑73 (MQ=255)
             tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                            <  2:1368157/76‑1 (MQ=255)
              cATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                           <  1:1721424/76‑1 (MQ=255)
              cATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                           <  1:1873752/76‑1 (MQ=255)
               aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                                          >  2:4090838/1‑76 (MQ=255)
               aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                                          <  2:4644367/76‑1 (MQ=255)
                ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGc                                                         <  2:2788220/76‑1 (MQ=255)
                ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGc                                                         <  1:1567158/76‑1 (MQ=255)
                 tGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                        <  2:2296842/76‑1 (MQ=255)
                 tGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                        >  1:3378816/1‑76 (MQ=255)
                 ggCCGTCACTGAAAGTGCTGCCGCCAGCAGCTTTTTGCGGATCATAGTCTTCCGTCAAGAAAGGATCGGGCCAGCa                                                        <  2:3124083/75‑1 (MQ=255)
                  gCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAt                                                       <  2:1938393/76‑1 (MQ=255)
                   ccGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                      >  2:1501322/1‑76 (MQ=255)
                    cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                     >  1:676783/1‑76 (MQ=255)
                    cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                     >  1:4868769/1‑76 (MQ=255)
                     gTCACTGCAATTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGTCCATCATTTg                                                    <  1:4522989/76‑1 (MQ=255)
                        aCTGAAAGTGCTGCCGCCAGCAGCATTTTGGGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                 <  2:2181319/76‑1 (MQ=255)
                        aCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                 <  1:94108/76‑1 (MQ=255)
                          tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGc                                               <  1:3063872/76‑1 (MQ=255)
                            aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                                             <  2:2184031/76‑1 (MQ=255)
                              aGTGCTGCCGCCGGCGGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                           <  2:4333585/76‑1 (MQ=255)
                              aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                           <  1:1300417/76‑1 (MQ=255)
                              aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAGATCAGGCCAGCATTTGCTGGCCCCa                                           <  2:1191990/76‑1 (MQ=255)
                               gTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTGAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAg                                          >  1:2321525/1‑76 (MQ=255)
                                tgctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCCAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                         >  2:3296895/1‑76 (MQ=255)
                                tgctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                         >  2:4860722/1‑76 (MQ=255)
                                tgctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                         >  2:3268685/1‑76 (MQ=255)
                                 gctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                            <  2:3057824/72‑1 (MQ=255)
                                 gctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                        <  1:2343926/76‑1 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      <  2:2959440/76‑1 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      <  2:3320409/76‑1 (MQ=255)
                                    gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                     <  1:1532157/76‑1 (MQ=255)
                                    gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                     >  1:3741380/1‑76 (MQ=255)
                                     ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                    <  2:2430997/76‑1 (MQ=255)
                                     ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                    >  1:3240601/1‑76 (MQ=255)
                                     ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGAt                                    >  1:3971948/1‑76 (MQ=255)
                                      cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                   <  1:1429788/76‑1 (MQ=255)
                                      cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                   <  2:2468205/76‑1 (MQ=255)
                                      cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                   <  1:4415406/76‑1 (MQ=255)
                                        ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAAc                                 >  1:1733247/1‑76 (MQ=255)
                                        ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAAc                                 >  1:3926900/1‑76 (MQ=255)
                                         cagcagCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATGGATAACa                                >  2:1426722/1‑76 (MQ=255)
                                           gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                              <  1:525575/76‑1 (MQ=255)
                                            cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                             >  1:289780/1‑76 (MQ=255)
                                            cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                             <  2:3741380/76‑1 (MQ=255)
                                            cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                             <  1:577666/76‑1 (MQ=255)
                                             agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                            <  2:2479270/76‑1 (MQ=255)
                                             agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                            >  2:2859382/1‑76 (MQ=255)
                                                aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                         >  1:249364/1‑76 (MQ=255)
                                                aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                         >  2:400801/1‑76 (MQ=255)
                                                aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                         >  2:3691632/1‑76 (MQ=255)
                                                aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                         >  2:2875591/1‑76 (MQ=255)
                                                aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                         >  1:4382244/1‑76 (MQ=255)
                                                 ttgtGCGCATCATAGTCTTCCCTCAAGAAAAAAACAGGCCCGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgc                        <  2:3261367/73‑1 (MQ=255)
                                                 ttatGCGCATCATAGACTTCCCTCAAGAAAAAAACAGGCCAGCATTTGCTGGCGCCAGATTGATAACAAAGTAcgc                        <  2:4533329/73‑1 (MQ=255)
                                                  tttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                       <  2:2587909/76‑1 (MQ=255)
                                                  tttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                       <  1:4030458/76‑1 (MQ=255)
                                                  tttGCGCAGCATAGTCTTCCCTCAAGAAAAAAGCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                       <  2:1361469/76‑1 (MQ=255)
                                                   ttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                      <  2:4163779/76‑1 (MQ=255)
                                                    tGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcgct                     >  2:1809307/1‑74 (MQ=255)
                                                    tGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt                     >  1:3366020/1‑76 (MQ=255)
                                                     gcgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                    <  1:360447/76‑1 (MQ=255)
                                                        catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                 <  1:103065/76‑1 (MQ=255)
                                                         atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                >  1:2189537/1‑76 (MQ=255)
                                                         atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                >  1:1328891/1‑76 (MQ=255)
                                                         atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGCTTTCCa                >  2:2679666/1‑76 (MQ=255)
                                                          tcatAGTCTCCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCAt               <  1:405397/76‑1 (MQ=255)
                                                           catTGTCTTCCCTCAAGAAAAGCTGAGGCTGGCGTTTGCTGGCCCCAGATTGATACCAAAGTACGCTTTTTGCATg              <  2:887787/76‑1 (MQ=255)
                                                                tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa         >  2:798955/1‑76 (MQ=255)
                                                                tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa         >  1:1720456/1‑76 (MQ=255)
                                                                tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTGCCATGCCGGa         >  2:3808211/1‑76 (MQ=255)
                                                                 ctacccTCAAGAAAAGAGCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt        <  1:3664864/73‑1 (MQ=255)
                                                                 cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt        >  2:2968799/1‑76 (MQ=255)
                                                                 cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt        <  1:642955/76‑1 (MQ=255)
                                                                  ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCGGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg       <  1:3335911/76‑1 (MQ=255)
                                                                  ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTCCCATGCCGGCTg       >  2:16920/1‑76 (MQ=255)
                                                                   tCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgc      >  1:912075/1‑76 (MQ=255)
                                                                     ccgcAAGAAAAAGTGAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcgg    <  2:2112961/73‑1 (MQ=38)
                                                                     ccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcgg    <  2:4433165/76‑1 (MQ=255)
                                                                      cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCCTGCCGGATgcggc   >  2:2042425/1‑76 (MQ=38)
                                                                      cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggc   >  1:1302170/1‑76 (MQ=255)
                                                                      cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggc   <  1:119861/76‑1 (MQ=255)
                                                                       tCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcg  <  2:4372289/76‑1 (MQ=255)
                                                                            |                                                                         
TCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCG  >  NC_000913/3099581‑3099726

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