Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,992 C→T 6.4% E191E (GAG→GAA)  fbaB ← fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,9920CT 6.4% 197.3 / 2.6 62E191E (GAG→GAA) fbaBfructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/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 = 6.87e-01

TTTTGTTTGACGATATCTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGCGTTTGATCGTAGGTATTCG  >  NC_000913/2177849‑2178140
                                                                                                                                               |                                                                                                                                                    
ttttGTTTGACGATATCTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGt                                                                                                                                                  <  2:303642/148‑1 (MQ=255)
  ttgtttgACGATATCTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTgc                                                                                                                                                <  1:116563/148‑1 (MQ=255)
     tttgACGATATCTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCAc                                                                                                                                            >  5:240248/1‑149 (MQ=255)
      ttgACGATATCTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCCCAGCACTGTCACCATACCCAGCTCGTGCGCACg                                                                                                                                           >  6:33207/1‑149 (MQ=255)
                cTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGcc                                                                                                                                 >  1:155049/1‑149 (MQ=255)
                   cACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGcc                                                                                                                                 <  5:161712/146‑1 (MQ=255)
                   cACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCg                                                                                                                                <  8:88558/147‑1 (MQ=255)
                   cACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCg                                                                                                                                <  1:228087/147‑1 (MQ=255)
                            ttGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATttct                                                                                                                     <  3:57982/149‑1 (MQ=255)
                            ttGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATttct                                                                                                                     <  6:203165/149‑1 (MQ=255)
                                    aGATTGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTg                                                                                                             >  7:117650/1‑149 (MQ=255)
                                    aGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTg                                                                                                             >  7:116823/1‑149 (MQ=255)
                                    aGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTg                                                                                                             <  5:33207/149‑1 (MQ=255)
                                     gATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTgg                                                                                                            >  6:111540/1‑149 (MQ=255)
                                      aTGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGc                                                                                                           >  2:7291/1‑149 (MQ=255)
                                       tGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCg                                                                                                          >  3:4907/1‑149 (MQ=255)
                                          tttGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCACCGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCGCACCGCCCTGTCACCATACCCAGCTCGTGCGCACGTTCAAACGCCGCACAAATATCTTCAATCTGGGGACg                                                                                                       >  4:85279/1‑149 (MQ=255)
                                            tGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGt                                                                                                      <  1:197720/148‑1 (MQ=255)
                                            tGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGt                                                                                                      <  5:111540/148‑1 (MQ=255)
                                             gCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTg                                                                                                     <  4:4907/148‑1 (MQ=255)
                                               cTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAAttt                                                                                                                       <  3:82786/128‑1 (MQ=255)
                                               cTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAAttt                                                                                                                       >  4:82786/1‑128 (MQ=255)
                                                   ccGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAATGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTc                                                                                              <  4:312254/149‑1 (MQ=255)
                                                      gTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGa                                                                                           <  6:100142/149‑1 (MQ=255)
                                                               gcggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCaaaa                                                                                    >  5:232586/1‑147 (MQ=255)
                                                                 ggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGTTCGTGCGCACGTTCAAAAGCCGCAGAAAttt                                                                                                                       >  3:221190/1‑110 (MQ=255)
                                                                 ggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGTTCGTGCGCACGTTCAAAAGCCGCAGAAAttt                                                                                                                       <  2:329855/110‑1 (MQ=255)
                                                                 ggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGTTCGTGCGCACGTTCAAAAGCCGCAGAAAttt                                                                                                                       <  4:221190/110‑1 (MQ=255)
                                                                 ggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGTTCGTGCGCACGTTCAAAAGCCGCAGAAAttt                                                                                                                       >  1:329855/1‑110 (MQ=255)
                                                                 ggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAataga                                                                                <  3:15378/149‑1 (MQ=255)
                                                                 ggAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAataga                                                                                <  6:232586/149‑1 (MQ=255)
                                                                   aaaCATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAatagata                                                                              <  1:7291/149‑1 (MQ=255)
                                                                    aaCATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAatagatag                                                                             <  4:258912/149‑1 (MQ=255)
                                                                         ggTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCa                                                                        >  7:89130/1‑149 (MQ=255)
                                                                          gTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCAc                                                                       >  2:330179/1‑149 (MQ=255)
                                                                             aTCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACcaa                                                                    >  6:329811/1‑149 (MQ=255)
                                                                                aaCGCAATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTAGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACcaacc                                                                  <  8:183346/148‑1 (MQ=255)
                                                                                       tctttcttAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCaaccgcaaccgc                                                          >  2:145602/1‑149 (MQ=255)
                                                                                          ttcttAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGcc                                                       >  6:28435/1‑149 (MQ=255)
                                                                                            cttAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGccc                                                      <  6:240248/148‑1 (MQ=255)
                                                                                             ttAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCAt                                                    >  6:35138/1‑149 (MQ=255)
                                                                                               aaaGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGt                                                  >  5:312591/1‑149 (MQ=255)
                                                                                                   gCGGAGTTACTCATTTAGGCCCTCAGCACTGTCTACATACCCAGCTCGTGCGCACGTTCAAAAGCCTGAGATTTTTCTTCAATCTGGCTACTTGACTATTGCTAGCCAAAATAGATAGTCGCACCAACCGAAACCGATACCATGTTGaa                                              <  8:117650/149‑1 (MQ=255)
                                                                                                   gCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGaa                                              >  1:160420/1‑149 (MQ=255)
                                                                                                   gCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGaa                                              <  6:66383/149‑1 (MQ=255)
                                                                                                    cGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAAc                                             <  2:59402/149‑1 (MQ=255)
                                                                                                      gAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGc                                           <  8:116823/149‑1 (MQ=255)
                                                                                                         ttACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTg                                        >  3:258030/1‑149 (MQ=255)
                                                                                                             gCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTcc                                    <  4:116365/149‑1 (MQ=255)
                                                                                                             gCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTc                                     >  8:224639/1‑148 (MQ=255)
                                                                                                                   aGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAatag                                                                                 <  3:11318/98‑1 (MQ=255)
                                                                                                                   aGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAatag                                                                                 >  4:11318/1‑98 (MQ=255)
                                                                                                                     gCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTgg                             >  1:186450/1‑148 (MQ=255)
                                                                                                                              aCTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTGCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGCGt                   <  1:228659/149‑1 (MQ=255)
                                                                                                                                tGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGc                     <  7:35804/145‑1 (MQ=255)
                                                                                                                                tGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGc                     <  7:30969/145‑1 (MQ=255)
                                                                                                                                tGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGc                     >  8:30969/1‑145 (MQ=255)
                                                                                                                                tGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGc                     >  8:35804/1‑145 (MQ=255)
                                                                                                                                  tCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACCTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGCGTTTGa               >  4:81928/1‑149 (MQ=255)
                                                                                                                                        tACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCATCAATCTGGCGACGGGACTCTTCCGAGCCAACATAGATAGTCGCACCAACCGGAACCGCGCCCATGTTGCAAGCCTGCTCCACGCTGGCATACAGCGTTTGATCGTAg         <  5:103515/149‑1 (MQ=255)
                                                                                                                                           ccAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGCGTTTGATCGTAGGTa      >  7:79615/1‑149 (MQ=255)
                                                                                                                                           ccAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCCCAACCGCGCCCATGTTGAACGCCTGCTCCAAGCTGGCATACAGCGTTAGATCGTAGGTa      >  5:88418/1‑149 (MQ=255)
                                                                                                                                               cTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGCGTTTGATCGTAGGTATTCg  <  2:155049/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTTTGTTTGACGATATCTGCACCGATGGTTGCCGCCAGATGGTTTGCCTGACCGGTCAGGTCGGCGGAAACATGGTAATCAACGCCATCTTTCTTAAAGGCGGAGTTACGCAAATAGGCCCACAGCACTGTCACCATACCCAGCTCGTGCGCACGTTCAAAAGCCGCAGAAATTTCTTCAATCTGGCGACGTGACTCTTCCGAGCCAAAATAGATAGTCGCACCAACCGCAACCGCGCCCATGTTGAACGCCTGCTCCACGCTGGCATACAGCGTTTGATCGTAGGTATTCG  >  NC_000913/2177849‑2178140

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