Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,054 (C)7→8 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,0471.C97.6% 288.4 / 1.2 83coding (1286/2526 nt)ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (42/39);  total (43/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 = 8.25e-01

CATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307898‑308192
                                                                                                                                                      |                                                                                                                                                 
cATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGc                                                                                                                                                  >  1:1493648/1‑151 (MQ=255)
 aTCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGcc                                                                                                                                                 >  2:1557804/1‑151 (MQ=255)
   cGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccc                                                                                                                                                <  2:1698375/150‑1 (MQ=255)
    gCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccc                                                                                                                                              <  1:763684/151‑1 (MQ=255)
      tACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccc                                                                                                                                            >  2:878017/1‑151 (MQ=255)
           cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccccaagc                                                                                                                                       >  2:1697323/1‑146 (MQ=255)
           cGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGccccccccaag                                                                                                                                        <  1:878530/150‑5 (MQ=255)
            gCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCg                                                                                                                                       <  2:1669554/151‑1 (MQ=255)
            gCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCg                                                                                                                                       >  2:592748/1‑151 (MQ=255)
             cAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCgg                                                                                                                                      >  1:226229/1‑151 (MQ=255)
              aaTTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCg                                                                                                                                       >  2:228183/1‑149 (MQ=255)
               attgattgCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGc                                                                                                                                     >  1:726005/1‑150 (MQ=255)
                      gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTc                                                                                                                               >  2:582447/1‑149 (MQ=255)
                      gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCa                                                                                                                              <  1:535209/150‑1 (MQ=255)
                      gCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAg                                                                                                                             <  2:1712978/151‑1 (MQ=255)
                       ccAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGa                                                                                                                            <  2:1539981/151‑1 (MQ=255)
                        cAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGAc                                                                                                                           >  1:1225448/1‑151 (MQ=255)
                         aaTGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACg                                                                                                                          >  2:1448016/1‑151 (MQ=255)
                         aaTGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACg                                                                                                                          >  2:25041/1‑151 (MQ=255)
                         aaTGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACg                                                                                                                          <  1:1040554/151‑1 (MQ=255)
                          aTGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACgg                                                                                                                         <  1:1683692/151‑1 (MQ=255)
                          aTGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACgg                                                                                                                         <  1:1241394/151‑1 (MQ=255)
                            gCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACggg                                                                                                                        >  2:1068000/1‑150 (MQ=255)
                               ggTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTc                                                                                                                    <  2:131163/151‑1 (MQ=255)
                                   ttttCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCCGCGTCAGACGGGTTTcaccc                                                                                                                >  2:445755/1‑151 (MQ=255)
                                      tCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGCTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacccac                                                                                                              <  1:50024/150‑1 (MQ=255)
                                       ccTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcacccacc                                                                                                             <  1:1364401/150‑1 (MQ=255)
                                          gATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGcc                                                                                                          <  1:236587/150‑1 (MQ=255)
                                          gATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCt                                                                                                         >  1:632182/1‑151 (MQ=255)
                                          gATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCt                                                                                                         >  2:790542/1‑151 (MQ=255)
                                          gATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCtt                                                                                                        <  1:798845/151‑2 (MQ=255)
                                           aTTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTg                                                                                                        <  2:724128/151‑1 (MQ=255)
                                             tAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAt                                                                                                      >  1:1287683/1‑151 (MQ=255)
                                             tAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGAt                                                                                                      >  2:1506880/1‑151 (MQ=255)
                                                cccACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTAt                                                                                                   >  1:315876/1‑151 (MQ=255)
                                                   acaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                                                                                >  2:651198/1‑151 (MQ=255)
                                                   acaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                                                                                >  1:1230123/1‑151 (MQ=255)
                                                   acaGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGt                                                                                                >  2:493745/1‑151 (MQ=255)
                                                     aGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTAt                                                                                              <  2:1366918/151‑1 (MQ=255)
                                                      gCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTAt                                                                                              <  1:547814/150‑1 (MQ=255)
                                                           cTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCatat                                                                                        >  1:1508860/1‑151 (MQ=255)
                                                             aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCGCCCACCGCCTGATTATCGTATCCATATc                                                                                       <  1:1389633/150‑1 (MQ=255)
                                                             aaaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTCGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATcc                                                                                      >  2:1389589/1‑151 (MQ=255)
                                                              aaGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATccc                                                                                     <  2:1572960/151‑1 (MQ=255)
                                                                gCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCg                                                                                    <  1:1596985/150‑1 (MQ=255)
                                                                 ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                                                                  <  1:680128/151‑1 (MQ=255)
                                                                 ccTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGACCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGtt                                                                                  <  2:601895/151‑1 (MQ=255)
                                                                    cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCCGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgcc                                                                               >  2:482084/1‑151 (MQ=255)
                                                                    cccGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgcc                                                                               <  2:49158/151‑1 (MQ=255)
                                                                       ggTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTgccgc                                                                             <  1:1525718/150‑1 (MQ=255)
                                                                        gTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGccc                                                                           <  1:326028/151‑1 (MQ=255)
                                                                          agagTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAg                                                                         <  2:1412321/151‑1 (MQ=255)
                                                                               ggCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                                                                     <  2:1008879/150‑1 (MQ=255)
                                                                                   cTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGt                                                                 >  2:1172061/1‑150 (MQ=255)
                                                                                    tGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc                                                               >  2:686189/1‑151 (MQ=255)
                                                                                     gATGCTGCCGATGCTGCTCCATTAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACt                                                              <  2:1485367/151‑1 (MQ=255)
                                                                                     gATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACt                                                              >  2:837941/1‑151 (MQ=255)
                                                                                           gcCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                        <  2:463230/151‑1 (MQ=255)
                                                                                             cGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc                                                      <  1:1631623/151‑1 (MQ=255)
                                                                                               aTGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCg                                                    >  2:1675599/1‑151 (MQ=255)
                                                                                                 gctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg                                                   <  2:77666/150‑1 (MQ=255)
                                                                                                  ctgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCgg                                                   >  2:970325/1‑149 (MQ=255)
                                                                                                   tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGa                                                  >  1:487335/1‑149 (MQ=255)
                                                                                                   tgctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAgg                                                >  2:1455755/1‑151 (MQ=255)
                                                                                                    gctCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAg                                                 <  1:25040/149‑1 (MQ=255)
                                                                                                      tCCATGAGCTATCACTGGCCAGCATATTTTGCAGGTTAACGGTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCTTATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc                                              >  2:965050/1‑150 (MQ=255)
                                                                                                       ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc                                              <  1:1067916/149‑1 (MQ=255)
                                                                                                       ccATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCg                                             >  1:165254/1‑150 (MQ=255)
                                                                                                          tGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGc                                         >  1:89155/1‑151 (MQ=255)
                                                                                                          tGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGc                                         >  1:868702/1‑151 (MQ=255)
                                                                                                              cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCt                                     >  2:511057/1‑151 (MQ=255)
                                                                                                                gTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCt                                     <  2:726062/149‑1 (MQ=255)
                                                                                                                  cACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCt                                  >  1:1648617/1‑150 (MQ=255)
                                                                                                                   aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAg                                <  2:1640648/151‑1 (MQ=255)
                                                                                                                      ggCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCa                             <  1:1447917/151‑1 (MQ=255)
                                                                                                                      ggCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCa                             <  1:971738/151‑1 (MQ=255)
                                                                                                                       gCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCa                             >  1:229559/1‑150 (MQ=255)
                                                                                                                        ccAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAAc                           >  2:1290742/1‑151 (MQ=255)
                                                                                                                         cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAAct                          >  2:842026/1‑151 (MQ=255)
                                                                                                                              tatTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCtt                     >  2:246190/1‑151 (MQ=255)
                                                                                                                               atTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCttt                    >  2:1056837/1‑150 (MQ=255)
                                                                                                                                  ttGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGc                 >  2:8453/1‑151 (MQ=255)
                                                                                                                                    gCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTg               >  1:1227312/1‑151 (MQ=255)
                                                                                                                                     cAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTgg              <  2:1225543/151‑1 (MQ=255)
                                                                                                                                       ggTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTc            <  2:1493749/151‑1 (MQ=255)
                                                                                                                                         ttAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTcgc          >  2:958727/1‑151 (MQ=255)
                                                                                                                                          tAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGt                                               <  1:151471/113‑1 (MQ=255)
                                                                                                                                          tAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGt                                               >  2:151480/1‑113 (MQ=255)
                                                                                                                                                 gATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCttttt  <  1:582390/151‑1 (MQ=255)
                                                                                                                                                 gATCGCCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCtttt   <  1:1117713/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTT  >  NC_000913/307898‑308192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: