Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,612,094 T→A 6.2% T78S (ACC→TCC)  yneG ← DUF4186 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,612,0940TA 6.2% 199.5 / 3.8 65T78S (ACC→TCC) yneGDUF4186 family protein
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base A (2/2);  total (36/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 = 6.73e-01

TGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTTCTCCAGGCAA  >  NC_000913/1611948‑1612239
                                                                                                                                                  |                                                                                                                                                 
tGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTc                                                                                                                                                  <  1:170069/148‑1 (MQ=255)
 gTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTcgc                                                                                                                                                >  6:216568/1‑149 (MQ=255)
      tttCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTgg                                                                                                                                           >  3:292144/1‑149 (MQ=255)
          cGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGt                                                                                                                                        <  7:239380/148‑1 (MQ=255)
              aTTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGa                                                                                                                                   <  5:124692/149‑1 (MQ=255)
               ttCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAg                                                                                                                                  <  6:372242/149‑1 (MQ=255)
                       gCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAg                                                                                                                                  <  8:75522/141‑1 (MQ=255)
                       gCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAg                                                                                                                                  >  7:75522/1‑141 (MQ=255)
                       gCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATaaa                                                                                                                           <  6:150072/148‑1 (MQ=255)
                               aTCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGAt                                                                                                                   <  1:223717/148‑1 (MQ=255)
                                 cTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAAAACCCTGCGTAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGcc                                                                                                                >  3:198425/1‑149 (MQ=255)
                                    tACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTgg                                                                                                                                           <  6:10362/119‑1 (MQ=255)
                                    tACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTgg                                                                                                                                           >  5:10362/1‑119 (MQ=255)
                                    tACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCcgc                                                                                                             <  8:315228/149‑1 (MQ=255)
                                         ccagccagTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATg                                                                                                        >  2:209425/1‑149 (MQ=255)
                                              cagTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGt                                                                                                   <  5:216568/149‑1 (MQ=255)
                                                 tGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTg                                                                                                >  8:56484/1‑149 (MQ=255)
                                                   aTAGATTACCGCCCCGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCt                                                                                              >  1:80952/1‑149 (MQ=255)
                                                   aTAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCt                                                                                              >  5:366059/1‑149 (MQ=255)
                                                       aTTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTAcc                                                                                          >  6:134949/1‑149 (MQ=255)
                                                        ttACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCg                                                                                         >  2:361716/1‑149 (MQ=255)
                                                         tACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTg                                                                                                                                       <  4:18729/102‑1 (MQ=255)
                                                         tACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTg                                                                                                                                       >  3:18729/1‑102 (MQ=255)
                                                              ccACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGt                                                                                                                                        <  6:188138/96‑1 (MQ=255)
                                                              ccACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGt                                                                                                                                        >  5:188138/1‑96 (MQ=255)
                                                              ccACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCAt                                                                                     >  4:104832/1‑147 (MQ=255)
                                                                             gttgttCTTCTCTTAACGCAACACCCTGCGGAATGTTATGCAACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCTGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCgg                                                                     <  4:198425/148‑1 (MQ=255)
                                                                              ttgttCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTGGCTTACCGTCATTAGCCGGTAATGCCGGAg                                                                   >  4:72270/1‑149 (MQ=255)
                                                                                gttCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGAGAGCAGGTCGCGGTGGCGTGATGAGAGATAAACACCGGATGCCCGCGCATGGGTGTTTGCTTACCTTAATTAGACGTTAAATCCTGAGc                                                                  >  5:294486/1‑148 (MQ=255)
                                                                                gttCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGc                                                                  >  6:230714/1‑148 (MQ=255)
                                                                                gttCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGc                                                                  <  1:361716/148‑1 (MQ=255)
                                                                                    ttcGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGAc                                                             >  1:295639/1‑149 (MQ=255)
                                                                                     tcGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGACGCgg                                                                                                                                              >  4:352458/1‑67 (MQ=255)
                                                                                     tcGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGACGCgg                                                                                                                                              >  4:346685/1‑67 (MQ=255)
                                                                                     tcGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGACGCgg                                                                                                                                              <  3:346685/67‑1 (MQ=255)
                                                                                     tcGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGACGCgg                                                                                                                                              <  3:352458/67‑1 (MQ=255)
                                                                                               aaaCACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGAGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGAc                                                   >  5:174873/1‑148 (MQ=255)
                                                                                                   acCCTGCGGAATGTTATGCCACTTAGCCAGCCAACCGCGACAGCAGGTCGCTGTGGCGTGCTGCGCGATAAACACCGGCTGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGAAATGCCGGAGCAAGACGTTTAGCGACCAAAt                                              >  2:145318/1‑149 (MQ=255)
                                                                                                     ccTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATcc                                            >  8:269226/1‑149 (MQ=255)
                                                                                                     ccTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATcc                                            >  3:308842/1‑149 (MQ=255)
                                                                                                         cGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATccgc                                          <  5:134949/147‑1 (MQ=255)
                                                                                                          ggAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCcgccgc                                       <  3:104832/149‑1 (MQ=255)
                                                                                                           gAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGAGACAAAATCCGCCGCa                                      >  4:241780/1‑149 (MQ=255)
                                                                                                            aaTGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCAt                                     >  1:170016/1‑149 (MQ=255)
                                                                                                                 tATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCAtgttgt                                >  3:313753/1‑149 (MQ=255)
                                                                                                                  aTGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATgttgtt                               >  3:13499/1‑149 (MQ=255)
                                                                                                                   tGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGCCAACATCCGCCGCATGTTGTTc                              >  8:160840/1‑149 (MQ=255)
                                                                                                                   tGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATgttgtt                               <  7:56484/148‑1 (MQ=255)
                                                                                                                   tGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATgttgtt                               <  7:147398/148‑1 (MQ=255)
                                                                                                                   tGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTc                              >  7:384146/1‑149 (MQ=255)
                                                                                                                   tGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTCGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATgttgtt                               <  3:72270/148‑1 (MQ=255)
                                                                                                                    gCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTg                                                                                                <  7:232567/82‑1 (MQ=255)
                                                                                                                    gCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTg                                                                                                >  8:232567/1‑82 (MQ=255)
                                                                                                                            gCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGAAGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTgg                     <  4:292144/149‑1 (MQ=255)
                                                                                                                              cAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGgcg                   >  1:380570/1‑149 (MQ=255)
                                                                                                                               aGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGa                                                              <  6:34657/105‑1 (MQ=255)
                                                                                                                               aGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGa                                                              >  5:34657/1‑105 (MQ=255)
                                                                                                                               aGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGgcgc                  <  1:209425/149‑1 (MQ=255)
                                                                                                                                 aCAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGgcgcgc                <  5:230714/149‑1 (MQ=255)
                                                                                                                                           cagcagGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTTCTCCAg      <  4:13499/149‑1 (MQ=255)
                                                                                                                                           cagcagGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTTCTCCAg      >  1:230724/1‑149 (MQ=255)
                                                                                                                                              cagGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTtctc         >  4:170561/1‑143 (MQ=255)
                                                                                                                                              cagGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTtctc         <  3:170561/143‑1 (MQ=255)
                                                                                                                                              cagGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGCGGGTTTGCTAACCGTCATTAGCCGGTAAAGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTTCTCCAGGc    >  6:133145/1‑148 (MQ=255)
                                                                                                                                               agGTCGCGGTGGCGTGCTGAGCTATAAACACCGTCTGCCCGCGCATGGGGGTTGGCGTCCCGTCCTTATCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTTGCGCGCCTTTCTCCAGGCaa  <  1:145318/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGTCATTAGCCGGTAATGCCGGAGCAAGACGTTTAGCGACAAAATCCGCCGCATGTTGTTCGATGACTGGCGCGCCTTTCTCCAGGCAA  >  NC_000913/1611948‑1612239

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: ‑