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

AGATTATTCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGTGAT  >  NC_000913/999604‑999797
                                                                                                    |                                                                                               
agaTTATGCTGGATGCCAGCGTAGCGTTTAAGGTGAATAAACGCCGTCTGGAAACTTCAGTACCGCAAATCGCCATCAAAAGTATTGCCCGTGGTTATGt                                                                                                  <  1:2300612/100‑1 (MQ=255)
 gaTTATTCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTc                                                                                                >  1:547073/1‑100 (MQ=255)
  attattCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTcc                                                                                               >  1:1776882/1‑100 (MQ=255)
   ttattCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTccc                                                                                             <  1:1254102/100‑1 (MQ=255)
    tattCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTcccc                                                                                            >  2:1064670/1‑100 (MQ=255)
     attCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTccccc                                                                                           <  2:867181/100‑1 (MQ=255)
       tCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTccccccc                                                                                         >  1:2074612/1‑99 (MQ=255)
       tCCGGATGCCAGCGTAACGTTTAATGTCAATAAACTCCGTTGGGTAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGGGGTTATGTccccccc                                                                                         <  2:2477094/100‑2 (MQ=255)
        ccGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTccccccct                                                                                        <  1:2594867/100‑3 (MQ=255)
        ccGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTccccccct                                                                                        >  1:1407166/1‑98 (MQ=255)
        ccGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTccccccct                                                                                        <  1:744049/100‑3 (MQ=255)
           gATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATGGCCGTCAAAAGTAATGCCCGTGGTTATGTccccccctgaa                                                                                     <  2:2435371/100‑6 (MQ=255)
             tGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCGGTGGTTATGTCCCCCCCTGAACg                                                                                    >  2:2485817/1‑100 (MQ=255)
             tGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACg                                                                                    <  1:1662305/100‑1 (MQ=255)
               ccAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGtt                                                                                  >  2:61042/1‑100 (MQ=255)
                cAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTg                                                                                 >  1:820187/1‑100 (MQ=255)
                 aGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTgg                                                                                <  2:1753002/100‑1 (MQ=255)
                 aGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTgg                                                                                <  1:480426/100‑1 (MQ=255)
                 aGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTgg                                                                                <  1:534540/100‑1 (MQ=255)
                   cGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGa                                                                              <  2:2125237/100‑1 (MQ=255)
                   cGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGa                                                                              <  2:2125562/100‑1 (MQ=255)
                    gTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGAt                                                                             >  1:332560/1‑100 (MQ=255)
                    gTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGAt                                                                             >  1:2124910/1‑100 (MQ=255)
                    gTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAACGCCCGTGGTTCTGTCCCCCCCTGAACGTTGCGAt                                                                             >  1:1011131/1‑100 (MQ=255)
                     tAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATg                                                                            >  1:1826204/1‑100 (MQ=255)
                        cGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAg                                                                         <  1:1250898/100‑1 (MQ=255)
                         gTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAgg                                                                        >  1:1693282/1‑100 (MQ=255)
                         gTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAgg                                                                        >  1:2643939/1‑100 (MQ=255)
                            tAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCATAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGAt                                                                     >  1:2560214/1‑100 (MQ=255)
                             aaTGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATc                                                                    <  1:304870/100‑1 (MQ=255)
                               tgtgAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCaa                                                                  >  2:979831/1‑100 (MQ=255)
                                gtgAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAAc                                                                 >  2:2182982/1‑100 (MQ=255)
                                gtgAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCTGAACGTTGGGATGAAGGGATCAAcg                                                                >  1:1337366/1‑100 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             <  2:2069838/100‑1 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             <  2:2070072/100‑1 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             <  1:53249/100‑1 (MQ=255)
                                    aTAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACgcgc                                                             >  1:577351/1‑100 (MQ=255)
                                     tAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCt                                                            <  1:2557991/100‑1 (MQ=255)
                                     tAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCt                                                            <  1:76943/100‑1 (MQ=255)
                                       aaCTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTAt                                                          >  2:14898/1‑100 (MQ=255)
                                        aCTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTAtt                                                         <  2:1670299/100‑1 (MQ=255)
                                         cgccGTCGGGAAAGTTCAGGACCGCAAATCGCCATCAAAAGTAATTCCCGTGGTTATGTCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTAc                                                       <  1:1314419/98‑1 (MQ=255)
                                         cTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTa                                                        <  1:2453036/100‑1 (MQ=255)
                                         cTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTAc                                                       >  2:1704649/1‑100 (MQ=255)
                                           ccGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACt                                                      <  2:1801840/100‑1 (MQ=255)
                                            cGTCTGGAAAGTTCAGTACCGCAAATCGCCATCAAAAGTAATGCGCGTGGTTATGTCCCCCCCTGAACGTTGGGATGTAGGGATCAACGCGCTATTACTg                                                     <  2:1002291/100‑1 (MQ=255)
                                             gTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTgg                                                    <  2:2178262/100‑1 (MQ=255)
                                             gTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTgg                                                    <  1:2123704/100‑1 (MQ=255)
                                              tCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTggg                                                   <  2:677485/100‑1 (MQ=255)
                                                tGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGat                                                 >  1:1562999/1‑100 (MQ=255)
                                                tGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGat                                                 >  2:123729/1‑100 (MQ=255)
                                                tGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGat                                                 >  1:30428/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTCCTGGGata                                                >  2:274554/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGata                                                >  2:2338135/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGata                                                >  1:547949/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGata                                                >  2:2026931/1‑100 (MQ=255)
                                                 ggAAATTTCAGTACCGCAAACCGCCATCAAAAGTAATGCCCGTTGTTATGTCCCCCCCTGAACGATGGAATGAAGGGATCAACGCGCTATTACTGGGata                                                >  2:655187/1‑100 (MQ=255)
                                                  gAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGatat                                               >  2:1768706/1‑100 (MQ=255)
                                                   aaaTTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATAtt                                              >  2:1170598/1‑100 (MQ=255)
                                                   aaaTTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATAtt                                              <  2:1497899/100‑1 (MQ=255)
                                                    aaTTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTc                                             >  2:1198588/1‑100 (MQ=255)
                                                      tttCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGt                                           >  2:1776212/1‑100 (MQ=255)
                                                      tttCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGt                                           >  1:266014/1‑100 (MQ=255)
                                                       ttCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGtt                                          >  2:50759/1‑100 (MQ=255)
                                                          aGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGc                                      <  1:351278/100‑1 (MQ=255)
                                                           gTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCg                                     >  1:1744777/1‑100 (MQ=255)
                                                             aCCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGTTATTACTGGGATATTCGTTTAGCgg                                    >  1:771139/1‑100 (MQ=255)
                                                              ccGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGACGCAGCGAACAACGCGCTATTACTGGGATATTCGTTTaaaggc                                   >  2:757083/1‑95 (MQ=255)
                                                                 cAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCt                                >  2:243190/1‑100 (MQ=255)
                                                                  aaaTCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTa                               <  2:922487/100‑1 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           >  2:1036777/1‑100 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           <  1:1597764/100‑1 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           >  2:2275586/1‑100 (MQ=255)
                                                                      cGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAg                           <  1:380181/100‑1 (MQ=255)
                                                                       gCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAgn                          <  2:645845/100‑2 (MQ=255)
                                                                       gCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGt                          <  1:2263891/100‑1 (MQ=255)
                                                                        ccATCAAAAGTAATGCCCGTGGTTATGTCCACCACTGAAAGTTGGGATGAAGGGATCAACGCGGTATTACTGGGATATTCGTTTAGCGGGGCTAAAAGTa                         >  2:450758/1‑100 (MQ=255)
                                                                          aTCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTAtt                       >  2:1333784/1‑100 (MQ=255)
                                                                           tCAAAAGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTc                      <  2:56273/100‑1 (MQ=255)
                                                                             aaaaGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCAt                    <  2:322433/100‑1 (MQ=255)
                                                                                aGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagc                 <  2:805723/100‑1 (MQ=255)
                                                                                aGTAATGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagc                 <  1:1979013/100‑1 (MQ=255)
                                                                                   aaTGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagcagc              <  2:925275/100‑1 (MQ=255)
                                                                                   aaTGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATagcagc              <  1:2330049/100‑1 (MQ=255)
                                                                                    aTGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAgcg             <  2:87359/100‑1 (MQ=255)
                                                                                     tGCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAgcgc            <  1:1201249/100‑1 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  1:1086039/1‑100 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  2:804912/1‑100 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  1:1992235/1‑100 (MQ=255)
                                                                                      gCCCGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCa           >  2:563354/1‑100 (MQ=255)
                                                                                       cccGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAg          <  2:1199003/100‑1 (MQ=255)
                                                                                       cccGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAg          <  1:2221160/100‑1 (MQ=255)
                                                                                        ccGTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGa         >  1:143129/1‑100 (MQ=255)
                                                                                          gTGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGAcn       <  1:867494/100‑2 (MQ=255)
                                                                                           tGGTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAg      >  1:1679577/1‑100 (MQ=255)
                                                                                            ggTTATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGt     <  1:1461553/100‑1 (MQ=255)
                                                                                               tATGTCCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGTGAt  >  1:1993515/1‑100 (MQ=255)
                                                                                                    |                                                                                               
AGATTATTCCGGATGCCAGCGTAACGTTTAATGTGAATAAACTCCGTCTGGAAATTTCAGTACCGCAAATCGCCATCAAAAGTAATGCCCGTGGTTATGTCCCCCCTGAACGTTGGGATGAAGGGATCAACGCGCTATTACTGGGATATTCGTTTAGCGGGGCTAACAGTATTCATAGCAGCGCAGACAGTGAT  >  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: