Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,862,763 T→C 6.5% H143H (CAT→CAC)  ygbK → FliA‑regulated DUF1537 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,862,7630TC 6.5% 183.6 / 2.7 61H143H (CAT→CAC) ygbKFliA‑regulated DUF1537 family protein
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base C (2/2);  total (29/32)
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

ATATTGGCCCGGTTACCGATGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAG  >  NC_000913/2862615‑2862907
                                                                                                                                                    |                                                                                                                                                
atatTGGCCCGGTTACCGATGCCTTAATGGATGCTCTCGACCCGCCGATTACGGTCGTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGTGGTATTTGTTCGTAATGAATCAACTACTGGCCGAATCCGGGATGCGCCAt                                                                                                                                                  >  6:22942/1‑149 (MQ=255)
       cccGGTTACCGATGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCACCACCCgg                                                                                                                                           >  6:117355/1‑149 (MQ=255)
       cccGGTTACCGATGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCACCACCCgg                                                                                                                                           >  6:113862/1‑149 (MQ=255)
         cGGTTACCGATGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGt                                                                                                                                          >  7:120189/1‑148 (MQ=255)
                 gATGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCAt                                                                                                                                 <  8:151309/149‑1 (MQ=255)
                  aTGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCAt                                                                                                                                 >  3:17283/1‑148 (MQ=255)
                   tGCCTTACTGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGTCCCGTCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGa                                                                                                                               <  2:196067/149‑1 (MQ=255)
                      cTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGAcc                                                                                                                             >  3:74849/1‑148 (MQ=255)
                      cTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGAcc                                                                                                                             >  8:167680/1‑148 (MQ=255)
                      cTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGAcc                                                                                                                             >  8:159033/1‑148 (MQ=255)
                            ggATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGc                                                                                                                       >  7:15094/1‑148 (MQ=255)
                               tGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGCCGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATc                                                                                                                   >  1:184151/1‑149 (MQ=255)
                                  tctcGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTc                                                                                                                >  3:106288/1‑149 (MQ=255)
                                    tcGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTccc                                                                                                              >  5:64578/1‑149 (MQ=255)
                                        cacGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGGTTATCAGGGGGATTTGTTCGTAATGAATCAACTGCTGGCCGAAAAAGGGATGCGCCATCCACCGGTAAAGCCCAGGTCCGCCCGAGATCTTCCCCGTc                                                                                                          >  2:184130/1‑149 (MQ=255)
                                          cGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTg                                                                                                        <  8:15094/149‑1 (MQ=255)
                                           gCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTgg                                                                                                       <  1:1020/149‑1 (MQ=255)
                                           gCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTgg                                                                                                       <  3:24554/149‑1 (MQ=255)
                                           gCCGATTACGGTCTTCTCTCCGGCCCTGCCGCTCAACGGACGTACGGTTTATCAGGGGTATTTGGACGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACATCTATCATCCCCGTCTgg                                                                                                       >  2:148409/1‑149 (MQ=255)
                                            ccGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGt                                                                                                      <  6:64578/149‑1 (MQ=255)
                                            ccGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCACCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGt                                                                                                      <  5:117355/149‑1 (MQ=255)
                                            ccGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCACCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGt                                                                                                      <  5:113862/149‑1 (MQ=255)
                                             cGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGtt                                                                                                     >  6:118464/1‑149 (MQ=255)
                                                 tACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGaa                                                                                                  <  5:61452/148‑1 (MQ=255)
                                                 tACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGaa                                                                                                  <  7:94752/148‑1 (MQ=255)
                                                           tctcCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCca                                                                                        >  7:76787/1‑148 (MQ=255)
                                                             tcCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCGTCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGCAGCGCAATCCACAg                                                                                     >  3:161427/1‑149 (MQ=255)
                                                             tcCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAg                                                                                     >  2:109116/1‑149 (MQ=255)
                                                               cGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCcaca                                                                                      <  8:230609/146‑1 (MQ=255)
                                                               cGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCcaca                                                                                      <  4:74849/146‑1 (MQ=255)
                                                               cGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAggg                                                                                   <  7:159033/149‑1 (MQ=255)
                                                               cGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAg                                                                                     <  8:97994/147‑1 (MQ=255)
                                                               cGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGGAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAgg                                                                                    <  2:108240/148‑1 (MQ=255)
                                                                      gCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTgcg                                                                            >  4:38684/1‑149 (MQ=255)
                                                                       ccGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTgc                                                                             <  1:225631/147‑1 (MQ=255)
                                                                          gTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGcgt                                                                        >  3:201344/1‑149 (MQ=255)
                                                                              aCGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGCATCAACTGCTGGCCGAATCCGGGCTGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGtt                                                                    <  2:41623/149‑1 (MQ=255)
                                                                              aCGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGCCAGCTATCTTCCCCGGCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGtt                                                                    >  2:150068/1‑149 (MQ=255)
                                                                                ggACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGAc                                                                                                                          <  7:96248/93‑1 (MQ=255)
                                                                                ggACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGAc                                                                                                                          >  8:96248/1‑93 (MQ=255)
                                                                                ggACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTc                                                                  <  5:118464/149‑1 (MQ=255)
                                                                                    gTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGGTTCGGca                                                              >  2:239540/1‑149 (MQ=255)
                                                                                             aTCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGtttt                                                      <  1:212146/148‑1 (MQ=255)
                                                                                              tCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGGCGTTTCGGCACATGTTTTc                                                     <  7:69575/148‑1 (MQ=255)
                                                                                                       tttGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgtg                                           <  2:184151/149‑1 (MQ=255)
                                                                                                         tGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgtgt                                          <  8:76787/148‑1 (MQ=255)
                                                                                                          gTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGtgtg                                         <  4:218389/148‑1 (MQ=255)
                                                                                                           ttCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTgg                                        <  1:109116/148‑1 (MQ=255)
                                                                                                           ttCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGa                                       >  8:136684/1‑149 (MQ=255)
                                                                                                               tAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGcc                                   >  5:105484/1‑149 (MQ=255)
                                                                                                                aaTGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCg                                  >  4:216964/1‑149 (MQ=255)
                                                                                                                 aTGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGt                                 <  3:38684/149‑1 (MQ=255)
                                                                                                                    actccaCTGCTGGCCGACTCCGCGATGCGCCATCACCCGTTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTGGAAGCGCAATCCACCGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGGGTGGATGCCGTTCg                              <  3:28397/144‑1 (MQ=255)
                                                                                                                      tCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGt                             >  1:19184/1‑148 (MQ=255)
                                                                                                                        aaCTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCaa                          <  7:167680/149‑1 (MQ=255)
                                                                                                                        aaCTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCaa                          <  3:138471/149‑1 (MQ=255)
                                                                                                                        aaCTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGGCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCa                           <  3:160307/148‑1 (MQ=255)
                                                                                                                          ctgctgGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAaga                        <  8:120189/149‑1 (MQ=255)
                                                                                                                              tgGCCGAATCCGGGATGCGCCATCACCCGGTAAATCACATGACCGACCGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGTGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTg                    >  3:134453/1‑149 (MQ=255)
                                                                                                                                  cGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCa                           >  1:76045/1‑138 (MQ=255)
                                                                                                                                          ggctGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAg        <  5:22942/146‑1 (MQ=255)
                                                                                                                                                gcCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAagag  >  7:148212/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATATTGGCCCGGTTACCGATGCCTTAATGGATGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAG  >  NC_000913/2862615‑2862907

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