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.A99.2% 420.2 / NA 133intergenic (‑14/+169)yggN/b2959hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base A (66/65);  minor base T (0/1);  total (67/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

TCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCG  >  NC_000913/3099581‑3099726
                                                                            |                                                                      
tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAg                                                                         >  1:1279116/1‑76 (MQ=255)
 cGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCCTTTTGCGCATCATAGTCTTCCCTCAAGa                                                                       >  2:2927675/1‑76 (MQ=255)
  ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                      <  2:4490214/76‑1 (MQ=255)
  ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaa                                                                      <  2:2959938/76‑1 (MQ=255)
    cGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaa                                                                    >  1:3179216/1‑76 (MQ=255)
     gTGGGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                   <  2:345799/76‑1 (MQ=255)
     gTGAGCGGTCTGTGCCTCCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                   <  1:871469/76‑1 (MQ=255)
     gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                   <  2:290199/76‑1 (MQ=255)
     gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaa                                                                   >  1:4480642/1‑76 (MQ=255)
      tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                  >  2:837812/1‑75 (MQ=255)
      tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                  <  2:5213884/76‑2 (MQ=255)
      tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                  <  2:3036535/76‑2 (MQ=255)
      tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                  >  1:1690470/1‑75 (MQ=255)
      tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                  >  1:1409793/1‑75 (MQ=255)
       gAGCGGTCATTGCCGTCACTGAGAGAGCTGGCGCCTGGAGCATTGGGCGAATCATAGTCTGCCATAAAGaaaaaat                                                                 >  2:2365918/1‑74 (MQ=255)
       gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaat                                                                 <  1:2576351/76‑3 (MQ=255)
        aGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCACCATAGTCTTCCCTCAAGaaaaaatc                                                                <  1:5519524/76‑4 (MQ=255)
           ggTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                              <  1:174746/76‑1 (MQ=255)
            gTCATGGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGc                                                             <  2:3308113/76‑1 (MQ=255)
             tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                            >  1:4846012/1‑76 (MQ=255)
             tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                            <  2:4337349/76‑1 (MQ=255)
             tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                            <  1:3776356/76‑1 (MQ=255)
             tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                            >  1:2052947/1‑76 (MQ=255)
              cATTGCCGTCACTGAGAGTTCTGCCGCTGGCATCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                           <  1:3703098/76‑1 (MQ=255)
              cATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                           <  2:1322270/76‑1 (MQ=255)
              cATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                           <  2:2032193/76‑1 (MQ=255)
              cATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTAAAGAAAAAATCAGGCCa                                                           <  1:1604904/76‑1 (MQ=255)
               aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                                          <  1:1357210/76‑1 (MQ=255)
                  gCCGTCACTTAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAt                                                       <  1:861261/76‑1 (MQ=255)
                  gCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAt                                                       >  1:4309963/1‑76 (MQ=255)
                  gCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAt                                                       >  2:3785560/1‑76 (MQ=255)
                  gCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCATCAAGAAAAAATCAGTCCAGcct                                                       >  2:1868461/1‑74 (MQ=255)
                    cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                     <  1:264739/76‑1 (MQ=255)
                    cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                     >  1:5448310/1‑76 (MQ=255)
                    cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                     >  1:4514946/1‑76 (MQ=255)
                    cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                     <  2:4307658/76‑1 (MQ=255)
                     gTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTg                                                    >  2:512760/1‑76 (MQ=255)
                       cACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                  <  1:5057506/76‑1 (MQ=255)
                       cACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                  <  1:3506482/76‑1 (MQ=255)
                        aCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                 <  2:1690470/76‑1 (MQ=255)
                         cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                                                <  2:4188897/76‑1 (MQ=255)
                          tGAAAGTGCTTCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCGCAAGAAAAAATCCGGCCGGCATTTTCTGGc                                               <  1:4997130/76‑1 (MQ=255)
                           gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcc                                              >  2:242345/1‑76 (MQ=255)
                           gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcc                                              >  2:1577734/1‑76 (MQ=255)
                           gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcc                                              >  1:5803845/1‑76 (MQ=255)
                            aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                                             >  1:4193599/1‑76 (MQ=255)
                              aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                           <  2:1305681/76‑1 (MQ=255)
                              aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                           >  1:2980964/1‑76 (MQ=255)
                               gTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAg                                          >  2:565239/1‑76 (MQ=255)
                               gTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAGTTGCTGGCCCCAg                                          >  1:2909862/1‑76 (MQ=255)
                                ggctgcCGCCACCGGCATTTTGCGCAACATATTCTGCCCTCAAGAAAAAAGCAGGCCAGCATTTGCTGGCCCCAGa                                         <  2:568845/75‑1 (MQ=255)
                                 gctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                        >  2:4416560/1‑76 (MQ=255)
                                  ctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                                       <  1:2571500/76‑1 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      >  2:3648904/1‑76 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      >  2:3315414/1‑76 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      >  1:5431996/1‑76 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      <  2:5501775/76‑1 (MQ=255)
                                   tgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTg                                      >  2:2724382/1‑76 (MQ=255)
                                   gcccgccAGCAGAATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAGCAGGACAGCATTTGCTGGCCCCAGATTg                                      >  2:4041208/3‑76 (MQ=255)
                                    gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                     <  2:1433613/76‑1 (MQ=255)
                                    gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                     >  1:4012343/1‑76 (MQ=255)
                                      cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGACACAGATTGATa                                   >  2:1501509/1‑76 (MQ=255)
                                       gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                  >  2:2519985/1‑76 (MQ=255)
                                       gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCAAGCATTTGCTGGCCCCAGATTGATaa                                  <  2:528393/76‑1 (MQ=255)
                                        ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAAc                                 >  1:587009/1‑76 (MQ=255)
                                        ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAAc                                 >  1:2958777/1‑76 (MQ=255)
                                         cagcagCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACa                                >  2:2468612/1‑76 (MQ=255)
                                           gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                              <  2:4766678/76‑1 (MQ=255)
                                             agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                            <  2:4865076/76‑1 (MQ=255)
                                             agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                            >  1:4251025/1‑76 (MQ=255)
                                              gcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                           >  1:3549362/1‑76 (MQ=255)
                                               caTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                          <  2:3007821/76‑1 (MQ=255)
                                                aTTTGGCGCATCATAGTCTTACATAAAGAAAAAATCAGGCAAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                         <  1:2136561/76‑1 (MQ=255)
                                                 ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgc                        >  2:2191951/1‑76 (MQ=255)
                                                  tttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                       <  2:3530064/76‑1 (MQ=255)
                                                   ttGCTCATCATCGTCTTGCGACAAGTAAGAATCAGGACAGGATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                      <  1:788451/76‑1 (MQ=255)
                                                   ttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                      >  2:1211093/1‑76 (MQ=255)
                                                   ttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                      <  2:3179216/76‑1 (MQ=255)
                                                    tGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt                     >  2:5172662/1‑76 (MQ=255)
                                                    tGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtt                     >  2:5023814/1‑76 (MQ=255)
                                                     gcgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                    >  2:383963/1‑76 (MQ=255)
                                                     gcgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                    >  1:1247839/1‑76 (MQ=255)
                                                     gcgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                    >  1:4006671/1‑76 (MQ=255)
                                                     gcgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                    <  1:3106815/76‑1 (MQ=255)
                                                      cgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtttt                   >  2:4738796/1‑76 (MQ=255)
                                                       gcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                  <  2:3875729/76‑1 (MQ=255)
                                                       gcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                  >  1:3891435/1‑76 (MQ=255)
                                                        catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                 >  1:5725201/1‑76 (MQ=255)
                                                          tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCAt               >  2:4293663/1‑76 (MQ=255)
                                                          tcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCAt               >  2:4350201/1‑76 (MQ=255)
                                                            ctAGTCTTCCCTCAAGAAAAAATCAGGCCAGCGTTTGCTGGCCCCAGATTGATAACAAAGTGCGGGTTTTCCATGc             <  2:3736604/75‑1 (MQ=255)
                                                             tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGcc            <  2:3610601/76‑1 (MQ=255)
                                                             tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGcc            <  1:774410/76‑1 (MQ=255)
                                                              aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCg           >  2:3674113/1‑76 (MQ=255)
                                                               gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCgg          <  2:3306877/76‑1 (MQ=255)
                                                                tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa         <  1:2415279/76‑1 (MQ=255)
                                                                tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa         >  2:4240535/1‑76 (MQ=255)
                                                                 cTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt        >  1:2358602/1‑76 (MQ=255)
                                                                    cccTCAAGAAAAAAGCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcg     <  2:4828634/76‑1 (MQ=38)
                                                                     ccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcgg    <  1:2013672/76‑1 (MQ=255)
                                                                            |                                                                      
TCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCG  >  NC_000913/3099581‑3099726

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