Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 561,688 G→A 6.7% T664T (ACG→ACA)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6880GA 6.7% 161.7 / 2.6 61T664T (ACG→ACA) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base A (2/2);  total (27/34)
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

CAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTG  >  NC_000913/561541‑561836
                                                                                                                                                   |                                                                                                                                                    
cAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTAGCACGTAAGCCAGGGTGCTACCGCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCGTACGGTACGc                                                                                                                                                     >  1:1042720/1‑149 (MQ=255)
   cAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgg                                                                                                                                                  >  1:1455981/1‑149 (MQ=255)
    aaCGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                                                                                 <  2:398653/149‑1 (MQ=255)
             aaTAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGta                                                                                                                                        >  2:1093378/1‑149 (MQ=255)
              aTAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGAAAGCCCGGGTGATACCTCGAATAATGGGTCCACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGtat                                                                                                                                       <  1:51127/149‑1 (MQ=255)
                  ctggctggCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACt                                                                                                                                   >  1:1095169/1‑149 (MQ=255)
                   tggctggcAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAATCCAGGGTGATACCGCGAATAATGGGTCCACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTTCGGTGCGCGGGGGGGCGGGTATTACTa                                                                                                                                  >  1:1029591/1‑149 (MQ=255)
                    ggcttgcAGGTGGAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATAACTAGAATAATGGGTACACGGGCACCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTAc                                                                                                                                 <  1:528243/149‑1 (MQ=255)
                    ggctgtcAGGTGTAGTTGGGACATTACTGGAAGGCCACTACCTGAGATATCACGTAAGCAAGGGTGATATCTCGATTAATGGGTACACGGGCAGCGCCACGGCAAACTGGCATGCCGCTTACGGTACGATGGGTGGCGGGTATAACTAc                                                                                                                                 <  2:687205/149‑1 (MQ=255)
                    ggctggcAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTAc                                                                                                                                 <  2:1115696/149‑1 (MQ=255)
                      ctggcAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGa                                                                                                                               >  2:1467443/1‑149 (MQ=255)
                              tgaagctGGGACCTAACTGCAAGTCCAAAACCTTAGTTATCACGTATGACAGGGTGATAACTCGAATAATGGGTACACTGGCGGCGCAACGGCAAAATGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATc                                                                                                                       <  1:1034942/143‑1 (MQ=255)
                                tAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCaa                                                                                                                     <  2:1550342/149‑1 (MQ=255)
                                gagggtggACCTTACTGGAAGGCCACATCCCGAGTTATCAAGTAAGCCAGGGTGATACCTCGTATAATGGGTAAAAGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCaa                                                                                                                     <  2:107727/143‑1 (MQ=255)
                                   gTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACACTGGGGGGCGGGTATAACTACGACCGCGATCAAc                                                                                                                    >  2:1460576/1‑147 (MQ=255)
                                   gTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACACTGGGGGGCGGGTATAACTACGACCGCGATCAAc                                                                                                                    >  2:1412256/1‑147 (MQ=255)
                                   gTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACACTGGGGGGCGGGTATAACTACGACCGCGATCAAc                                                                                                                    <  1:1412256/147‑1 (MQ=255)
                                   gTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACACTGGGGGGCGGGTATAACTACGACCGCGATCAAc                                                                                                                    <  1:1460576/147‑1 (MQ=255)
                                               ttgAAGGCCACTACCTGAGTGAGCACGTTCTCCAGGGCGATACCTCGATTTATGGGTACACGGGCAGCGCCAGGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGGATAACTACGACCGCGATCAACATGACGTTAACTgg                                                                                                      <  2:288054/147‑1 (MQ=255)
                                                 gAAGGCCTCAACCTGAGTTATCACGTAAGCCAGGGAGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                                                                                    <  1:801405/149‑1 (MQ=255)
                                                 gAAGGCCACAACCTGAGTTATCGCGTAAGCCAGGGTGATACATCGTATAATGGGTACACGGGCCGCGCCCCGGCAAACTGGCGGGCCGGTTTCGGTACGCTGGGGGGCGGGTATAACTACGATCGCGATCAACATGACGTTATCTGGcg                                                                                                    >  2:725756/1‑148 (MQ=255)
                                                 gAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTATACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACTTTAACTGGCa                                                                                                    >  2:1216478/1‑149 (MQ=255)
                                                 gAAGGCCACAACCTGAGTTAACACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCCGCGCCACGGCAAACAGGCAGTCCGCAGACCGTCCCCTGGGGTGCGGGTATTACTACGCCCGCGATCACCATGTCGGTAAcggccc                                                                                                    >  2:871629/1‑144 (MQ=255)
                                                       cacaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                                                                              >  2:1311357/1‑149 (MQ=255)
                                                          aacaggAGGTATCTCGTAGTCCAGGGAGATACCTCGAATAACGGGTACCCGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGCGCGGGTACAACTACTACCGCGATCAACATGACGTTAACTGGCAGCTGTCTgg                                                                                           <  2:1052822/144‑1 (MQ=255)
                                                          aaCCTGAGGTATCACGCAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAACCTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGGATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTgg                                                                                           <  1:1539341/149‑1 (MQ=255)
                                                              tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgt                                                                                       <  2:23317/149‑1 (MQ=255)
                                                              tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgt                                                                                       >  2:539979/1‑149 (MQ=255)
                                                                     gcACGAACGCCTGGCTGATACCACGAATCATGGGTCCCCGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                                                                                <  2:650844/148‑1 (MQ=255)
                                                                        cGTAAGCCAGGGTGATACATCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGGTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                                                                             >  1:437669/1‑149 (MQ=255)
                                                                             gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                                                                        <  1:1467443/149‑1 (MQ=255)
                                                                                  gttGATACCGCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGGCGGGCATGAAAATgg                                                                   <  2:744129/147‑1 (MQ=255)
                                                                                           tCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACAGCGATCAACATGACGTAAACTGGCCGCTGTATGGCGGGGTGAGCGGGCAGGACAATGGAATAACGCt                                                          >  2:1046133/1‑149 (MQ=255)
                                                                                             gAATAATGGGTACTCGGTCAGCGCCACGGGAATCTGGCAGGCCGCTTACGGTACGCGGGGGGGCGGGTATAACTACGTCCGCGATCAACATGACGTTAACTGGGATCTGTCTGGCGGTGTGGTCGGGCATGAAATTGGCATAACGCTGa                                                        <  1:418010/149‑1 (MQ=255)
                                                                                             gAATAATGGGGACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCGTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGa                                                        <  1:1169908/149‑1 (MQ=255)
                                                                                               ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGCTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc                                                      >  2:919743/1‑149 (MQ=255)
                                                                                               ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc                                                      <  2:436414/149‑1 (MQ=255)
                                                                                                 aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcca                                                    <  2:810448/149‑1 (MQ=255)
                                                                                                  atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGGACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagccag                                                   >  2:1375963/1‑149 (MQ=255)
                                                                                                   tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCTTAACGCTGagccagc                                                  >  1:108062/1‑149 (MQ=255)
                                                                                                    gggTCCACGGGCAGCGCCACGGCAAACGGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTTTACCTACGAGCGCGATCAAGATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagccagcc                                                 <  2:1219438/149‑1 (MQ=255)
                                                                                                    gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGCCCGCGATCAACATGCCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagccagcc                                                 >  2:736595/1‑149 (MQ=255)
                                                                                                     ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCTCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCt                                                <  1:428539/149‑1 (MQ=255)
                                                                                                     ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCAGGAAAATGGCATACCGCTGAGCCAGCCt                                                >  1:80673/1‑149 (MQ=255)
                                                                                                       tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCttt                                              >  1:787354/1‑149 (MQ=255)
                                                                                                         cccGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCTGGTATAACAACGACCGCGATCAACATGACCGTAACTCGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAg                                            <  2:479928/147‑1 (MQ=255)
                                                                                                         cacGGGCAGCGCCACGGCACACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAg                                            <  2:1455981/149‑1 (MQ=255)
                                                                                                              gCAGCGCCACGGCAAACTGGCATGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACCTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGAt                                       <  1:313990/149‑1 (MQ=255)
                                                                                                              gCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGAt                                       >  1:378754/1‑149 (MQ=255)
                                                                                                                aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAAGTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGGGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATAc                                     <  1:346509/149‑1 (MQ=255)
                                                                                                                   gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCaa                                  >  2:537112/1‑149 (MQ=255)
                                                                                                                    ccACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCCGCATTTAGGGGATACCAAt                                 >  2:1346380/1‑149 (MQ=255)
                                                                                                                     cACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATg                                >  1:562674/1‑149 (MQ=255)
                                                                                                                        ggCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGGATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGttt                             <  1:1368889/149‑1 (MQ=255)
                                                                                                                          cAAACTGGCAGGCCGCTTACGGTACGTTGTGTGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTg                           <  2:942718/149‑1 (MQ=255)
                                                                                                                          cAAACAGGCAGGCCGCTTACGGTACGCCGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTg                           <  1:400247/149‑1 (MQ=255)
                                                                                                                           aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGa                          >  1:1449647/1‑149 (MQ=255)
                                                                                                                             aCTGGCCGGCCGCTTACGGTACGCAGGGGGGCGGGTATAACTACGCCCGCGATCACCATGACGGTCCCTGGCAGCTGACTGGAGGGGGGGGCGGGCATGAAAAGGGCATAACGCTGAGCCAGCCTTTCGGGGATACCAATGTTTTGAtt                        <  2:1042720/149‑1 (MQ=255)
                                                                                                                             aCTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGAtt                        <  2:1537897/149‑1 (MQ=255)
                                                                                                                                     gCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGGGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACGAATGTTTTGATTAAAGCGcc                >  2:337600/1‑149 (MQ=255)
                                                                                                                                         cTTACGGTACGCGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGGCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGgc            <  2:567180/149‑1 (MQ=255)
                                                                                                                                                   gCTGGGGGGCGGTTATAACTACGACCGCGATCAACATGACTTTAACTGGCAGCTGTGTGTCGGTGTGGTCGGGCATGAAAATGGAATAACGATGAGCATGCATTTAGGGGATACAAATGTTTTGATTAATGCGCCTGGAGCAGGCGgtg  >  2:878280/1‑149 (MQ=255)
                                                                                                                                                   gCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTATGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAAAGTTTTGATTAAAGCGCCTGGCGCAGGCGgtg  >  2:322529/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAGCAACGGGCAAAATAGCTGGCTGGCAGGTGTAGGTGGGACCTTACTGGAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTG  >  NC_000913/561541‑561836

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