Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,054 (C)7→6 coding (1279/2526 nt) ecpC ← ECP production outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,0480C.97.5% 355.0 / 1.9 81coding (1285/2526 nt)ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (36/43);  total (37/44)
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.90e-01

CAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGG  >  NC_000913/307950‑308143
                                                                                                  |                                                                                               
caGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcc                                                                                                >  2:1016503/1‑100 (MQ=255)
 aGCGAACTAAAGCCTCCCGGTAGAGTGTCGCTGATGCTGCCGAGGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccc                                                                                               <  1:2721844/100‑1 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  2:2079427/1‑100 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  2:3199199/100‑1 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  1:1387154/1‑100 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              <  1:2872564/100‑1 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  2:3136837/1‑100 (MQ=255)
  gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccc                                                                                              >  1:307065/1‑100 (MQ=255)
    gAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccccc                                                                                            >  2:722563/1‑100 (MQ=255)
     aaCTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccccca                                                                                           <  1:45063/100‑2 (MQ=255)
       cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaag                                                                                         >  1:1152497/1‑97 (MQ=255)
       cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaag                                                                                         >  2:223715/1‑97 (MQ=255)
       cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcccccca                                                                                           >  2:1594392/1‑97 (MQ=255)
        tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaagc                                                                                        <  1:2624614/100‑5 (MQ=255)
        tAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccaagc                                                                                        >  1:727062/1‑96 (MQ=255)
         aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCg                                                                                      >  1:94379/1‑100 (MQ=255)
           aGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGCGCTGTCACTGGCCAGCATATTTTGCAGGATAACGTTGATCG‑CCCCCCAAGCGGc                                                                                    >  2:849412/1‑100 (MQ=255)
             ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAg                                                                                  <  2:2433854/100‑1 (MQ=255)
             ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATACTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAg                                                                                  >  1:361087/1‑100 (MQ=255)
              cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGc                                                                                 >  1:3158292/1‑100 (MQ=255)
              cTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGc                                                                                 >  1:2256440/1‑100 (MQ=255)
                nccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGt                                                                               >  1:2081523/2‑100 (MQ=255)
                 ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTc                                                                              >  1:2776022/1‑100 (MQ=255)
                 ccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTc                                                                              <  1:2389689/100‑1 (MQ=255)
                     tAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGAc                                                                          <  2:470028/100‑1 (MQ=255)
                      agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGAc                                                                          <  1:3022842/100‑1 (MQ=255)
                        ngTGGCGCTGATTCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCTGCAGCGTCAGACggg                                                                       >  2:936893/2‑100 (MQ=255)
                              gcTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTcac                                                                 <  2:1136127/100‑1 (MQ=255)
                               cTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTcacc                                                                <  1:2023492/100‑1 (MQ=255)
                                tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccc                                                               <  2:557467/100‑1 (MQ=255)
                                 gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccca                                                              >  1:1846689/1‑100 (MQ=255)
                                     ctgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGc                                                          <  2:2084984/100‑1 (MQ=255)
                                      tgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                         <  1:129905/100‑1 (MQ=255)
                                      tgcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                         >  1:3119822/1‑100 (MQ=255)
                                       gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCt                                                        <  1:1622156/100‑1 (MQ=255)
                                        ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                       <  2:651435/100‑1 (MQ=255)
                                        ccGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                       <  2:1846667/100‑1 (MQ=255)
                                         cGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGa                                                      >  1:693388/1‑100 (MQ=255)
                                            tgctgcTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTa                                                   >  2:3266835/1‑100 (MQ=255)
                                             gctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                  >  1:1337692/1‑100 (MQ=255)
                                              ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                 <  1:592235/100‑1 (MQ=255)
                                              ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATc                                                 >  1:128525/1‑100 (MQ=255)
                                               tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                <  2:2642646/100‑1 (MQ=255)
                                               tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCg                                                >  1:1861611/1‑100 (MQ=255)
                                                gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                               >  1:3211936/1‑100 (MQ=255)
                                                 ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTa                                              <  1:1009540/100‑1 (MQ=255)
                                                 ctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTa                                              <  1:2454513/100‑1 (MQ=255)
                                                   ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  1:491829/100‑1 (MQ=255)
                                                   ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  2:2535312/100‑1 (MQ=255)
                                                   ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATc                                            <  1:3045137/100‑1 (MQ=255)
                                                     aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:345802/100‑1 (MQ=255)
                                                     aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:1635938/100‑1 (MQ=255)
                                                     aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCa                                          <  1:217941/100‑1 (MQ=255)
                                                       gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        <  2:2372731/100‑1 (MQ=255)
                                                       gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        <  2:2497221/100‑1 (MQ=255)
                                                       gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCata                                        >  1:2318515/1‑100 (MQ=255)
                                                          cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  2:306372/100‑1 (MQ=255)
                                                          cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  1:1016515/100‑1 (MQ=255)
                                                          cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                     <  2:3251078/100‑1 (MQ=255)
                                                           tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                    <  2:2074855/100‑1 (MQ=255)
                                                           tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                    >  2:1074206/1‑100 (MQ=255)
                                                            gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                   >  1:733817/1‑100 (MQ=255)
                                                             tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGt                                  >  2:269496/1‑100 (MQ=255)
                                                              cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCCTATCCCTATCCCg                                   >  2:1626623/1‑99 (MQ=255)
                                                              cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                 <  2:1637243/100‑1 (MQ=255)
                                                              cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                 >  1:3219029/1‑100 (MQ=255)
                                                               aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                <  1:276624/100‑1 (MQ=255)
                                                               aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTg                                >  1:1369782/1‑100 (MQ=255)
                                                                cTGGCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATCATCGTATCCATATCCCGTTgc                               <  1:2407540/100‑1 (MQ=255)
                                                                atcggcAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgc                               <  2:1220548/95‑1 (MQ=255)
                                                                   gCCAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgc                            >  1:1675421/1‑100 (MQ=255)
                                                                    ccAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgcc                           <  1:1594409/100‑1 (MQ=255)
                                                                    ccAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccg                             >  2:233771/1‑98 (MQ=255)
                                                                     cAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          <  2:2256414/100‑1 (MQ=255)
                                                                     cAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          >  2:407979/1‑100 (MQ=255)
                                                                     cAGCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                          <  1:134316/100‑1 (MQ=255)
                                                                       gCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAg                        >  2:1491277/1‑100 (MQ=255)
                                                                           atTTTGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                    >  1:1740457/1‑100 (MQ=255)
                                                                            ttttGCAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAg                   <  2:2882011/100‑1 (MQ=255)
                                                                             cttccgatctAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                  <  1:943751/91‑1 (MQ=255)
                                                                                 cAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc              >  1:86954/1‑100 (MQ=255)
                                                                                 cAGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc              >  2:1888520/1‑100 (MQ=255)
                                                                                  aGGTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACt             <  1:177441/100‑1 (MQ=255)
                                                                                   ggTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTc            >  2:1034942/1‑100 (MQ=255)
                                                                                    gTTAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCa           <  2:1391734/100‑1 (MQ=255)
                                                                                      tAACGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAAt         <  2:588745/100‑1 (MQ=255)
                                                                                       aaCGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGTTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTTAGCGTACTCa           >  2:1385957/1‑97 (MQ=255)
                                                                                       aaCGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGcagatcggaag        >  2:943739/1‑90 (MQ=255)
                                                                                       aaCGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATg        <  1:3136874/100‑1 (MQ=255)
                                                                                        aCGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa       >  2:1129614/1‑100 (MQ=255)
                                                                                         cGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa       <  2:1547487/99‑1 (MQ=255)
                                                                                         cGTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAg      >  1:3170514/1‑100 (MQ=255)
                                                                                          gTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc     >  1:2924300/1‑100 (MQ=255)
                                                                                          gTTGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc     <  1:425718/100‑1 (MQ=255)
                                                                                            tGATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCg   >  1:2102095/1‑100 (MQ=255)
                                                                                             gATCG‑CCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg  <  2:727054/100‑1 (MQ=255)
                                                                                                  |                                                                                               
CAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGG  >  NC_000913/307950‑308143

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