Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 561,696:1 (G)6→7 100% coding (2000/2604 nt) sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6901.G65.1% 343.9 / 208.1 152coding (1994/2604 nt)sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base . (29/24);  new base G (49/50);  total (78/74)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561594‑561784
                                                                                                  |                                                                                                 
 ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggg                                                                                            >  1:4162862/1‑100 (MQ=255)
 ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggg                                                                                            >  1:1532599/1‑100 (MQ=255)
  cacaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCACTTACGGTACGCT‑‑ggggg                                                                                           <  2:2001551/100‑1 (MQ=255)
    caACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑ggggggg                                                                                         >  1:3507555/1‑99 (MQ=255)
     aaCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggggggc                                                                                        <  2:2659826/100‑3 (MQ=255)
      aCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggggggc                                                                                        >  1:2412453/1‑97 (MQ=255)
      aCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCgg                                                                                       <  2:4384914/100‑1 (MQ=255)
       ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggggggcgg                                                                                      <  2:859539/100‑5 (MQ=255)
        cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggggggcggg                                                                                     >  2:2933678/1‑95 (MQ=255)
        cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gggggggcggg                                                                                     <  1:4249760/100‑6 (MQ=255)
         tGAGTTATCACATAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGt                                                                                     <  2:3505330/100‑1 (MQ=255)
            gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATaa                                                                                 >  1:4267308/1‑100 (MQ=255)
             ttATCACGTAAGCCAGGGTGATACCTCGAATCATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑ggggggggcgggtct                                                                                 >  1:3517383/1‑90 (MQ=255)
              tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                                >  2:1414646/1‑100 (MQ=255)
              tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                                <  2:583703/100‑1 (MQ=255)
              tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                                <  1:4016505/100‑1 (MQ=255)
               aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACt                                                                               >  2:3717539/1‑100 (MQ=255)
               aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTa                                                                              <  2:4401531/100‑1 (MQ=255)
                tCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTa                                                                              <  2:373455/100‑1 (MQ=255)
                tCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTa                                                                              <  2:62515/100‑1 (MQ=255)
                 cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTAc                                                                             >  1:3402833/1‑100 (MQ=255)
                 cACGTAAGCCACGGGGATACCTCGAATAATGGGTACACGGGGAGCGCCACGGCCAACGGGTAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACt                                                                               <  1:586718/98‑1 (MQ=255)
                  aCGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                            >  1:4252119/1‑100 (MQ=255)
                  aCGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                            >  2:494468/1‑100 (MQ=255)
                  aCGTAAGCCAGGGTGATACCTCGAATAATGGGAACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                            >  2:4185294/1‑100 (MQ=255)
                    gTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                            <  2:4218552/98‑1 (MQ=255)
                    gTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGAcc                                                                         >  2:3014181/1‑100 (MQ=255)
                     tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAcc                                                                         >  2:4107250/1‑100 (MQ=255)
                      aaGCCATGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcg                                                                        >  1:786622/1‑100 (MQ=255)
                      aaGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACcgc                                                                       <  1:763202/100‑1 (MQ=255)
                       aGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACcgcg                                                                      <  2:2432630/100‑1 (MQ=255)
                        gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGa                                                                     >  2:1538780/1‑100 (MQ=255)
                         ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGa                                                                     >  1:1089314/1‑100 (MQ=255)
                         ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGa                                                                     >  2:2185637/1‑100 (MQ=255)
                               tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               >  1:2388761/1‑100 (MQ=255)
                               tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               >  1:825221/1‑100 (MQ=255)
                               tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  2:3215307/100‑1 (MQ=255)
                                 aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATg                                                             >  2:359747/1‑100 (MQ=255)
                                   aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  1:3613104/1‑100 (MQ=255)
                                    ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  2:3070476/1‑100 (MQ=255)
                                    ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          <  2:287457/100‑1 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  1:140357/1‑100 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         <  2:1356347/100‑1 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                        <  2:3556183/100‑1 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                        >  2:3305922/1‑100 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTAAGCTGGGGGGGCGGGTATAACTACGACCGGGATAAACATGACGt                                                         >  1:3955988/1‑100 (MQ=255)
                                       cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         <  2:2128678/98‑1 (MQ=255)
                                            aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTgg                                                  <  2:2801473/100‑1 (MQ=255)
                                            aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTgg                                                  >  2:4163773/1‑100 (MQ=255)
                                             atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                                <  2:2277702/100‑1 (MQ=255)
                                             atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                                <  1:2204448/100‑1 (MQ=255)
                                              tGGGTACACGGGCAGCGCCACGGCAAACGGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                                <  1:3015991/100‑1 (MQ=255)
                                               gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                               <  2:238486/100‑1 (MQ=255)
                                               gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                               >  1:3235720/1‑100 (MQ=255)
                                               gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                              <  2:2705617/100‑1 (MQ=255)
                                               gggTACACGGGCAGCGCCACGGCAAACTGCCAGGCCGCTTACGGTACGCTT‑‑GGGGGGCGGTTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                               <  1:3304986/100‑1 (MQ=255)
                                                ggTACACGGGCAGCGCCACGGCAAACTGGCGGGCCGCTGACGGTACGCTGGGGGGGCGGGTAGAACGACGACCGCGATCAACATGACGTTAACTGGCa                                                <  1:2870521/98‑1 (MQ=255)
                                                  tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                            <  2:1513147/100‑1 (MQ=255)
                                                  tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                            >  1:1276234/1‑100 (MQ=255)
                                                      cGGGCAGCGCCACGGCAACCTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTAGCAg                                               >  1:491075/1‑92 (MQ=255)
                                                          cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                    >  2:2048937/1‑100 (MQ=255)
                                                          cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                    <  1:2185630/100‑1 (MQ=255)
                                                          cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                    >  1:3136932/1‑100 (MQ=255)
                                                          cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                    >  1:1714437/1‑100 (MQ=255)
                                                          cAGCGCCACAGCAAACGGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                    <  2:1260601/99‑1 (MQ=255)
                                                           aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                                  <  1:570600/100‑1 (MQ=255)
                                                             cgcCACGGCAAACTTGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                                  <  1:3950850/99‑1 (MQ=255)
                                                             cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                 <  2:4313759/100‑1 (MQ=255)
                                                              gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTCTCTGGCGGtgtg                                >  1:1640991/1‑100 (MQ=255)
                                                                cACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGt                              <  1:292487/100‑1 (MQ=255)
                                                                cACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GCGGGGGCGTGTATAACTACGACCGCGATCAACAGGACGTTAACTGGCAGCTGGCTGGCGGggtgtt                              >  1:1730433/1‑98 (MQ=255)
                                                                cACGGCAAACTGGCAGGCCGCTTACGCTACGCTGGGGGGGGGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGggtggt                              >  2:2688098/1‑100 (MQ=255)
                                                                 aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                             >  1:1141798/1‑100 (MQ=255)
                                                                  cGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgg                           >  2:2493370/1‑100 (MQ=255)
                                                                   ggCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                          <  1:3763413/100‑1 (MQ=255)
                                                                    gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                          <  2:3277376/100‑1 (MQ=255)
                                                                    gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                          <  2:1714444/100‑1 (MQ=255)
                                                                     cAAACTGGCAGGCCGCCTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                         <  1:1068447/100‑1 (MQ=255)
                                                                      aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCa                        >  2:2375052/1‑100 (MQ=255)
                                                                      aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCa                        >  1:1511636/1‑100 (MQ=255)
                                                                        aCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATg                      >  2:3646848/1‑100 (MQ=255)
                                                                        aCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGACGGCGAGGTCGGGCATg                      >  1:3352044/1‑100 (MQ=255)
                                                                        aCTGGCAGGCCGATTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATTAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATg                      <  2:1918897/100‑1 (MQ=255)
                                                                            gCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa                  <  2:746549/100‑1 (MQ=255)
                                                                               ggCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc              >  1:1977489/1‑100 (MQ=255)
                                                                                gCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc              <  2:1662609/100‑1 (MQ=255)
                                                                                gCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc              <  2:1845477/100‑1 (MQ=255)
                                                                                gCCGCTTACGGTACGCT‑‑GGTGGGGCGGGTATCACTACGACCGGGATCAACAGGACGTTGACTGGCTGCCGTCTGGGGGTGCGGTCGGGCCTGAAAATTGc              >  1:1387569/1‑100 (MQ=255)
                                                                                 ccGCTTACGGTACGCTGGGGGGGGCGAGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc              >  2:2689950/1‑100 (MQ=255)
                                                                                  cggttACGGTACGCT‑‑GGGGGGGGGGTATAACGGCGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg                <  1:2790648/92‑1 (MQ=255)
                                                                                    cTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc         >  2:2523640/1‑100 (MQ=255)
                                                                                      tACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGc       >  1:1609529/1‑100 (MQ=255)
                                                                                       aCGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGc       <  2:3620066/100‑1 (MQ=255)
                                                                                       aCGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGc       <  1:1149896/100‑1 (MQ=255)
                                                                                        cGGTACGCTGGGGGGGCGGGTATAACTATGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCt      <  1:1473455/100‑1 (MQ=255)
                                                                                         ggTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTg     >  1:110186/1‑100 (MQ=255)
                                                                                           tACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGCGTGGTCGGGCATGAAAATGGCATAACGCTGag   <  2:4191611/100‑1 (MQ=255)
                                                                                           tACGCT‑‑GGGGGGCGGGTATAACTACGACCTCGATCAACATGACTTTAACTGGCAGCTGTCTGGCGGTGTTGTCGCGCATGAAAATGGCATAACGCTGagc  <  1:1432539/100‑1 (MQ=255)
                                                                                            aCGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc  <  2:2788916/100‑1 (MQ=255)
                                                                                            aCGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCt      >  1:3694175/1‑95 (MQ=255)
                                                                                                  |                                                                                                 
GCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561594‑561784

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