Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 999,709 (C)6→7 coding (494/2601 nt) elfC → putative outer membrane fimbrial subunit export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913999,7031.C94.1% 271.2 / 12.7 85coding (488/2601 nt)elfCputative outer membrane fimbrial subunit export usher protein
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (42/38);  total (45/40)
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.64e-01

AGATTATTCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑CC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGTGAT  >  NC_000913/999604‑999797
                                                                                                    |                                                                                               
agaTTATGCTGGATGCCAGCGTAGCGTTTAAGGTGAATAAACGCCGTCTGGAAACTTCAGTACCGCAAATCGCCATCAAAAGTATTGCCCGTGGTTATGt                                                                                                  <  1:2300612/100‑1 (MQ=255)
 gaTTATTCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑c                                                                                                >  1:547073/1‑100 (MQ=255)
  attattCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc                                                                                               >  1:1776882/1‑100 (MQ=255)
   ttattCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑c                                                                                             <  1:1254102/100‑1 (MQ=255)
    tattCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑cc                                                                                            >  2:1064670/1‑100 (MQ=255)
     attCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑ccc                                                                                           <  2:867181/100‑1 (MQ=255)
       tCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑ccccc                                                                                         >  1:2074612/1‑99 (MQ=255)
       tCCGGATGCCAGCGTAACGTTTAATGTCAATAAACTCCGTTGGGTAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGGGGTTATGT‑cc‑ccccc                                                                                         <  2:2477094/100‑2 (MQ=255)
        ccGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑ccccct                                                                                        <  1:2594867/100‑3 (MQ=255)
        ccGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑ccccct                                                                                        >  1:1407166/1‑98 (MQ=255)
        ccGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑cc‑ccccct                                                                                        <  1:744049/100‑3 (MQ=255)
           gATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATGGCCGTCAAAAGTAATGCCCGTGGTTATGT‑cc‑ccccctgaa                                                                                     <  2:2435371/100‑6 (MQ=255)
             tGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCGGTGGTTATGTCCC‑CCCCTGAACg                                                                                    >  2:2485817/1‑100 (MQ=255)
             tGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACg                                                                                    <  1:1662305/100‑1 (MQ=255)
               ccAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGtt                                                                                  >  2:61042/1‑100 (MQ=255)
                cAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTg                                                                                 >  1:820187/1‑100 (MQ=255)
                 aGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTgg                                                                                <  2:1753002/100‑1 (MQ=255)
                 aGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTgg                                                                                <  1:480426/100‑1 (MQ=255)
                 aGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTgg                                                                                <  1:534540/100‑1 (MQ=255)
                   cGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGa                                                                              <  2:2125237/100‑1 (MQ=255)
                   cGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGa                                                                              <  2:2125562/100‑1 (MQ=255)
                    gTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGAt                                                                             >  1:332560/1‑100 (MQ=255)
                    gTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGAt                                                                             >  1:2124910/1‑100 (MQ=255)
                    gTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAACGCCCGTGGTTCTGTCCC‑CCCCTGAACGTTGCGAt                                                                             >  1:1011131/1‑100 (MQ=255)
                     tAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATg                                                                            >  1:1826204/1‑100 (MQ=255)
                        cGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAg                                                                         <  1:1250898/100‑1 (MQ=255)
                         gTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAgg                                                                        >  1:1693282/1‑100 (MQ=255)
                         gTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAgg                                                                        >  1:2643939/1‑100 (MQ=255)
                            tAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCATAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGAt                                                                     >  1:2560214/1‑100 (MQ=255)
                             aaTGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATc                                                                    <  1:304870/100‑1 (MQ=255)
                               tgtgAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCaa                                                                  >  2:979831/1‑100 (MQ=255)
                                gtgAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAAc                                                                 >  2:2182982/1‑100 (MQ=255)
                                gtgAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑CC‑CCCCTGAACGTTGGGATGAAGGGATCAAcg                                                                >  1:1337366/1‑100 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             <  2:2069838/100‑1 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             <  2:2070072/100‑1 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             <  1:53249/100‑1 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             >  1:577351/1‑100 (MQ=255)
                                     tAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCt                                                            <  1:2557991/100‑1 (MQ=255)
                                     tAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCt                                                            <  1:76943/100‑1 (MQ=255)
                                       aaCTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTAt                                                          >  2:14898/1‑100 (MQ=255)
                                        aCTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTAtt                                                         <  2:1670299/100‑1 (MQ=255)
                                         cgccGTCGGGAAAGTTCAGGACCGCAAATCGCCATCAAAAGTAATTCCCGTGGTTATGT‑CC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTAc                                                       <  1:1314419/98‑1 (MQ=255)
                                         cTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTa                                                        <  1:2453036/100‑1 (MQ=255)
                                         cTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGT‑GTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTAc                                                       >  2:1704649/1‑100 (MQ=255)
                                           ccGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACt                                                      <  2:1801840/100‑1 (MQ=255)
                                            cGTCTGGAAAGTTCAGTACCGCAAATCGCCATCAAAAGTAATGCGCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGTAGGGATCAACGCGCTATTACTg                                                     <  2:1002291/100‑1 (MQ=255)
                                             gTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTgg                                                    <  2:2178262/100‑1 (MQ=255)
                                             gTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTgg                                                    <  1:2123704/100‑1 (MQ=255)
                                              tCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTggg                                                   <  2:677485/100‑1 (MQ=255)
                                                tGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGat                                                 >  1:1562999/1‑100 (MQ=255)
                                                tGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGat                                                 >  2:123729/1‑100 (MQ=255)
                                                tGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGat                                                 >  1:30428/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTCCTGGGata                                                >  2:274554/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGata                                                >  2:2338135/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGata                                                >  1:547949/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGata                                                >  2:2026931/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAACCGCCATCAAAAGTAATGCCCGTTGTTATGTCCC‑CCCCTGAACGATGGAATGAAGGGATCAACGCGCTATTACTGGGata                                                >  2:655187/1‑100 (MQ=255)
                                                  gAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGatat                                               >  2:1768706/1‑100 (MQ=255)
                                                   aaaTTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATAtt                                              >  2:1170598/1‑100 (MQ=255)
                                                   aaaTTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATAtt                                              <  2:1497899/100‑1 (MQ=255)
                                                    aaTTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTc                                             >  2:1198588/1‑100 (MQ=255)
                                                      tttCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGt                                           >  2:1776212/1‑100 (MQ=255)
                                                      tttCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGt                                           >  1:266014/1‑100 (MQ=255)
                                                       ttCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGtt                                          >  2:50759/1‑100 (MQ=255)
                                                          aGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑CC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGc                                      <  1:351278/100‑1 (MQ=255)
                                                           gTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑CC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCg                                     >  1:1744777/1‑100 (MQ=255)
                                                             aCCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGTTATTACTGGGATATTCGTTTAGCgg                                    >  1:771139/1‑100 (MQ=255)
                                                              ccGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGACGCAGCGAACAACGCGCTATTACTGGGATATTCGTTTaaaggc                                   >  2:757083/1‑95 (MQ=255)
                                                                 cAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCt                                >  2:243190/1‑100 (MQ=255)
                                                                  aaaTCGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTa                               <  2:922487/100‑1 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           >  2:1036777/1‑100 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           <  1:1597764/100‑1 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           >  2:2275586/1‑100 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           <  1:380181/100‑1 (MQ=255)
                                                                       gCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAgn                          <  2:645845/100‑2 (MQ=255)
                                                                       gCCATCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGt                          <  1:2263891/100‑1 (MQ=255)
                                                                        ccATCAAAAGTAATGCCCGTGGTTATGT‑CCACCACTGAAAGTTGGGATGAAGGGATCAACGCGGTATTACTGGGATATTCGTTTAGCGGGGCTAAAAGTa                         >  2:450758/1‑100 (MQ=255)
                                                                          aTCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTAtt                       >  2:1333784/1‑100 (MQ=255)
                                                                           tCAAAAGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTc                      <  2:56273/100‑1 (MQ=255)
                                                                             aaaaGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCAt                    <  2:322433/100‑1 (MQ=255)
                                                                                aGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagc                 <  2:805723/100‑1 (MQ=255)
                                                                                aGTAATGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagc                 <  1:1979013/100‑1 (MQ=255)
                                                                                   aaTGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagcagc              <  2:925275/100‑1 (MQ=255)
                                                                                   aaTGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagcagc              <  1:2330049/100‑1 (MQ=255)
                                                                                    aTGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAgcg             <  2:87359/100‑1 (MQ=255)
                                                                                     tGCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAgcgc            <  1:1201249/100‑1 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  1:1086039/1‑100 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  2:804912/1‑100 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  1:1992235/1‑100 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  2:563354/1‑100 (MQ=255)
                                                                                       cccGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAg          <  2:1199003/100‑1 (MQ=255)
                                                                                       cccGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAg          <  1:2221160/100‑1 (MQ=255)
                                                                                        ccGTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGa         >  1:143129/1‑100 (MQ=255)
                                                                                          gTGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGAcn       <  1:867494/100‑2 (MQ=255)
                                                                                           tGGTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAg      >  1:1679577/1‑100 (MQ=255)
                                                                                            ggTTATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGt     <  1:1461553/100‑1 (MQ=255)
                                                                                               tATGTCCC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGTGAt  >  1:1993515/1‑100 (MQ=255)
                                                                                                    |                                                                                               
AGATTATTCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGT‑CC‑CCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGTGAT  >  NC_000913/999604‑999797

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