Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,576 A→G 6.6% D154D (GAT→GAC)  cheZ ← chemotaxis regulator, protein phosphatase for CheY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,5760AG 6.6% 188.6 / 3.3 61D154D (GAT→GAC) cheZchemotaxis regulator, protein phosphatase for CheY
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACGGGCGTC  >  NC_000913/1966445‑1966722
                                                                                                                                   |                                                                                                                                                  
acacCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTaa                                                                                                                                    >  8:109965/1‑148 (MQ=255)
  acCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATg                                                                                                                                  <  7:102158/148‑1 (MQ=255)
    cGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGAc                                                                                                                                <  7:109965/148‑1 (MQ=255)
    cGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGAc                                                                                                                                <  8:49259/148‑1 (MQ=255)
    cGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGAc                                                                                                                                <  6:312047/148‑1 (MQ=255)
     ggCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGAcc                                                                                                                               >  7:239844/1‑148 (MQ=255)
           gCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGccc                                                                                                                          <  3:294758/147‑1 (MQ=255)
                  tCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGt                                                                                                                       >  5:209493/1‑143 (MQ=255)
                  tCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGAt                                                                                                                  <  5:144608/148‑1 (MQ=255)
                     aCCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCACCAACTGGCGATCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCt                                                                                                               >  8:131033/1‑148 (MQ=255)
                      ccTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGAGCCt                                                                                                               >  6:207565/1‑147 (MQ=255)
                             tCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATcc                                                                                                       >  5:152935/1‑148 (MQ=255)
                                aTTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCt                                                                                                      >  1:158292/1‑146 (MQ=255)
                                aTTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCt                                                                                                      <  2:158292/146‑1 (MQ=255)
                                aTTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCAACATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTgc                                                                                                    >  5:129347/1‑148 (MQ=255)
                                  tAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTgcgc                                                                                                  >  1:263434/1‑148 (MQ=255)
                                         aCTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACGTCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCcatcat                                                                                            >  3:37512/1‑147 (MQ=255)
                                         aCTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACGTCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCcatcat                                                                                            >  5:5000/1‑147 (MQ=255)
                                                 tttCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCcag                                                                                   >  1:116003/1‑148 (MQ=255)
                                                    cACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCcagca                                                                                 <  6:209493/147‑1 (MQ=255)
                                                    cACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCcagca                                                                                 <  6:152935/147‑1 (MQ=255)
                                                      cGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACGGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACCTCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGAGCCATCATGATTTCCAGCAGtt                                                                              >  5:288799/1‑148 (MQ=255)
                                                             gACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACGATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGtt                                                                       <  2:263434/148‑1 (MQ=255)
                                                              aCGCGACTCCTGTTCCGGGATGTTTTCCAACAGAACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTa                                                                      >  6:126790/1‑148 (MQ=255)
                                                                 cgACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGt                                                                    >  7:7440/1‑147 (MQ=255)
                                                                    cTCCTGTTCCGGTATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAg                                                                <  2:28882/148‑1 (MQ=255)
                                                                      ccTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCt                                                              >  5:75898/1‑148 (MQ=255)
                                                                          ttCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACGTCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTgg                                                            <  4:37512/146‑1 (MQ=255)
                                                                          ttCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACGTCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTgg                                                            <  6:5000/146‑1 (MQ=255)
                                                                          ttCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGTAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTgg                                                            <  1:325580/146‑1 (MQ=255)
                                                                           tCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTAt                                                         >  7:171898/1‑148 (MQ=255)
                                                                               ggATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATgcgc                                                     <  8:171898/148‑1 (MQ=255)
                                                                                   gTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCggg                                                  >  8:254549/1‑147 (MQ=255)
                                                                                     tttCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTa                                                >  6:54278/1‑147 (MQ=255)
                                                                                       tCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGa                                                                                          >  4:148288/1‑103 (MQ=255)
                                                                                       tCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGa                                                                                          <  3:148288/103‑1 (MQ=255)
                                                                                          aaCAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTg                                          >  4:153468/1‑148 (MQ=255)
                                                                                           aCAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGAATGCGCGGGTACATCTGc                                         >  6:141733/1‑148 (MQ=255)
                                                                                            cAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGcc                                        <  4:265657/148‑1 (MQ=255)
                                                                                                   aTCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGc                                         >  7:168112/1‑140 (MQ=255)
                                                                                                   aTCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGc                                         <  8:168112/140‑1 (MQ=255)
                                                                                                   aTCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAAtt                                 >  3:256429/1‑148 (MQ=255)
                                                                                                     cAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGtt                                                                       >  6:290166/1‑108 (MQ=255)
                                                                                                     cAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGtt                                                                       <  5:290166/108‑1 (MQ=255)
                                                                                                      aGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGt                               <  8:239844/147‑1 (MQ=255)
                                                                                                             ggCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGc                                                               <  8:40126/108‑1 (MQ=255)
                                                                                                             ggCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGc                                                               >  7:40126/1‑108 (MQ=255)
                                                                                                              gCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTAt                        <  8:7440/146‑1 (MQ=255)
                                                                                                                  tCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATc                       <  4:233353/143‑1 (MQ=255)
                                                                                                                  tCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATc                       >  3:233353/1‑143 (MQ=255)
                                                                                                                   cGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCGAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGttc                  <  5:126790/147‑2 (MQ=255)
                                                                                                                     aTCTCCTGAATGAAATCCATTATACGCTTTATGACATGCCCGGTGAGATACTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTAcc                <  8:126590/147‑1 (MQ=255)
                                                                                                                      actcCTGAATGACATCCATCATCCGCTTAATGACATGCCAGGTGAGATCCTCAAAATCCTGCGCCATCATGATTTCCAGCATATGCGCGTTACTAAAGCTGGTATGCGCGGGTACATCTGCAAGAAATAGTCGtgtaactgtaac                 >  7:126590/2‑145 (MQ=255)
                                                                                                                       ctcCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGt             >  5:304864/1‑148 (MQ=255)
                                                                                                                         ccTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGtt                   <  2:216232/140‑1 (MQ=255)
                                                                                                                         ccTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGtt                   <  2:211609/140‑1 (MQ=255)
                                                                                                                         ccTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGtt                   >  1:216232/1‑140 (MQ=255)
                                                                                                                         ccTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGtt                   >  1:211609/1‑140 (MQ=255)
                                                                                                                            gAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACg        >  5:280306/1‑148 (MQ=255)
                                                                                                                             aaTGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACg        <  6:129347/147‑1 (MQ=255)
                                                                                                                             aaTGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGAATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACg        <  5:141733/147‑1 (MQ=255)
                                                                                                                                  catccatcATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACGGGCGTc  <  7:131033/148‑1 (MQ=255)
                                                                                                                                  catccatcATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACGGGCGTc  <  2:116003/148‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                  
ACACCGGCTTTGCTGGTATCGACCTGAGGTCCATTAAGCAAACTCTGGTTTTCACGTTTTGGACGCGACTCCTGTTCCGGGATGTTTTCCAACAGCACCATCAGCAACTGGCGTTCGATCTCCTGAATGACATCCATCATCCGCTTAATGACCTGCCCGGTGAGATCCTGAAAATCCTGCGCCATCATGATTTCCAGCAGTTGCGCGTTAGTAAAGCTGGTATGCGCGGGTACATCTGCCAGAAATTGTCGTGTATCTGTTACCAGTTCACGGGCGTC  >  NC_000913/1966445‑1966722

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