Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,459,293:1 (G)5→6 100% coding (955/1482 nt) gspE → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,459,2881.G96.4% 473.8 / 9.1 139coding (950/1482 nt)gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (57/77);  total (61/78)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 155 total aligned reads displayed.

TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑GG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCA  >  NC_000913/3459189‑3459382
                                                                                                     |                                                                                               
tGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg‑t                                                                                                  <  1:3694116/100‑1 (MQ=255)
 ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGACATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg‑t‑g                                                                                                >  2:741452/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑g                                                                                             >  1:2158430/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑g                                                                                             >  2:3521524/1‑100 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑gg                                                                                            <  2:2233950/100‑1 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑gg                                                                                            <  2:1147618/100‑1 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑gg                                                                                            <  2:3658002/100‑1 (MQ=255)
     ggCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCTGATGTCGTCATGG‑T‑gg‑ggg                                                                                           <  1:1162273/100‑1 (MQ=255)
       cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑gggga                                                                                         >  2:860744/1‑98 (MQ=255)
       cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑gggga                                                                                         >  2:630093/1‑98 (MQ=255)
       cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑gggga                                                                                         >  2:2726096/1‑98 (MQ=255)
        aTTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gg‑ggggaa                                                                                        >  2:2311706/1‑97 (MQ=255)
        aTTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCGGGCGGCCTGCGCGCCATGCTTCGCCAGGAGCCGGATGTCGTCAGGG‑T‑GG‑GGGgaa                                                                                        >  1:1131070/1‑100 (MQ=255)
           gggCAAACGCAGGTGCATACACGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCTTGG‑TGGG‑GGGAAAtt                                                                                      >  1:772911/1‑100 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGt                                                                                   <  1:3363200/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCTTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGGCAGACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTg                                                                                  <  1:357948/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGTGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTg                                                                                  <  2:3774791/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTg                                                                                  <  1:725863/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTg                                                                                  <  1:2358105/100‑1 (MQ=255)
                agcGCAGGTGCACACCCGTGTGGTTATGTCGTTCGCTCGCGGCCGGCGCGCCATACGTCGCCCGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTg                                                                                  <  1:2094537/97‑1 (MQ=255)
                aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGa                                                                                 <  2:2168483/100‑1 (MQ=255)
                 aCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGAt                                                                                <  2:1138618/100‑1 (MQ=255)
                     aGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAg                                                                            <  2:3746396/100‑1 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGa                                                                           >  1:485655/1‑100 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGa                                                                           >  1:2379808/1‑100 (MQ=255)
                        tGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGaaa                                                                         >  2:2251084/1‑100 (MQ=255)
                          aaTACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAAcc                                                                       >  2:3400675/1‑100 (MQ=255)
                             aCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACcgcg                                                                    <  1:1215580/100‑1 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCgcgc                                                                   <  2:199488/100‑1 (MQ=255)
                                cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAg                                                                 >  2:1656614/1‑100 (MQ=255)
                                    tgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGCG‑GG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTg                                                             >  1:973749/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                            >  2:126119/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                            >  1:3487910/1‑100 (MQ=255)
                                      gATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATTCAGAAACCGCGCAGATTGCg                                                           <  1:2311440/100‑1 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                          <  2:756639/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                         <  1:246620/100‑1 (MQ=255)
                                         atGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGCAATTCGACATACAGAAACCGCGCAGATTGCGGt                                                         >  1:2182587/1‑99 (MQ=255)
                                            tCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAg                                                     <  1:2744430/100‑1 (MQ=255)
                                             cTTTCGCTCGTGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑GGTGGGAGATAAGTGATACAGAAACCGCGAGCATCGCGGTTCAgg                                                    >  1:2829092/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAGATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                    >  1:2556630/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                    <  1:142910/100‑1 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                    >  1:1359434/1‑100 (MQ=255)
                                              tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                   <  1:3475167/100‑1 (MQ=255)
                                              tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                   <  2:2529197/100‑1 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                  <  1:545334/100‑1 (MQ=255)
                                               gtcgctcgcGGCCTGCGCGCCATACGTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                  <  1:1968737/99‑1 (MQ=255)
                                                 cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                                >  1:1507298/1‑100 (MQ=255)
                                                  gctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCg                                               >  2:3087557/1‑100 (MQ=255)
                                                      gcgGCCTGCGCGCCATACTTCGCCAGGACCCGAATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                   <  1:1027948/92‑1 (MQ=255)
                                                       cgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAc                                          >  1:2140844/1‑100 (MQ=255)
                                                       cgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAc                                          <  2:2684421/100‑1 (MQ=255)
                                                        ggCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGAAGATTGCGGTTCAGGCCTCGCTCAcc                                         >  1:2734197/1‑100 (MQ=255)
                                                         gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                        <  1:39769/100‑1 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                       <  1:2853719/100‑1 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGGCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                       <  1:825308/100‑1 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGCCGTCAAGG‑TGGG‑GGGAAAATCGAGGTAAAGAACCGGACAAGACTGCGGTTCACGCCGCGCt                                             >  1:370224/1‑94 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc                                      >  2:2941333/1‑100 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                     >  2:1830445/1‑100 (MQ=255)
                                                                cgcCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGACTCAGGCCTCGCTCACCGGCCATCt                                 >  1:1958483/1‑100 (MQ=255)
                                                                 gcCATACTGCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATc                                  <  1:1280929/98‑1 (MQ=255)
                                                                  ccATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTg                                >  2:2065935/1‑99 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                             <  1:1208377/100‑1 (MQ=255)
                                                                     tACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAc                            <  2:1450362/100‑1 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          >  1:939349/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          <  2:2791740/100‑1 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          >  2:777461/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          >  1:1662238/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                       <  2:2844652/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                       <  2:2528290/100‑1 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      >  1:1131927/1‑100 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      >  2:65504/1‑100 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      <  2:2389405/100‑1 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      >  1:729015/1‑100 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCacc                     >  1:646247/1‑98 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                     <  2:3039074/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                     <  1:1507851/100‑1 (MQ=255)
                                                                              aGGACCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAacac                   >  2:3184094/1‑100 (MQ=255)
                                                                              aGGACCCGGATGTCGTCAGGG‑TGGG‑GGGAAATTCGTGATACAGACACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAaca                    <  1:2478262/99‑1 (MQ=255)
                                                                                 aCCCGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTcc                <  1:2036117/100‑1 (MQ=255)
                                                                                 aCCCGGATCTCGTCATGG‑TGGG‑GGGAAATTCGTGATTCAGAACCCGCACAGATCTCGCGTCAGGCCTCGCTCACCGACCATCTGGTACTCTCAACACTcc                <  1:2819616/100‑1 (MQ=255)
                                                                                  cccGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca               >  2:2249997/1‑100 (MQ=255)
                                                                                  cccGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca               >  2:1943692/1‑100 (MQ=255)
                                                                                   ccGGATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac              >  1:2401382/1‑100 (MQ=255)
                                                                                     ggATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            >  2:1415909/1‑100 (MQ=255)
                                                                                     ggATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            <  1:2856172/100‑1 (MQ=255)
                                                                                     ggATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            >  2:1598747/1‑100 (MQ=255)
                                                                                     ggATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            >  2:3209727/1‑100 (MQ=255)
                                                                                      ggtGTCGTCATGG‑TGGG‑GGGAAATTCGTGGTACAGACACCGCGCAGATTGCGGTCCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac              <  1:3136184/95‑1 (MQ=255)
                                                                                      gATGTCGTCATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACt           <  1:3184073/100‑1 (MQ=255)
                                                                                         gtcgtcATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc        >  1:2393789/1‑100 (MQ=255)
                                                                                         gtcgtcATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc        <  1:1065222/100‑1 (MQ=255)
                                                                                         gtcgtcATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc        <  1:385916/100‑1 (MQ=255)
                                                                                           cgtcATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg      >  2:1470828/1‑100 (MQ=255)
                                                                                            gtcATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGt     >  1:3796085/1‑100 (MQ=255)
                                                                                              cATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTg    <  1:1170668/99‑1 (MQ=255)
                                                                                              cATGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc   <  1:1737874/100‑1 (MQ=255)
                                                                                               aTGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTg    <  1:2301396/98‑1 (MQ=255)
                                                                                               aTGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc   <  1:3674575/99‑1 (MQ=255)
                                                                                               aTGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCa  <  2:772916/100‑1 (MQ=255)
                                                                                               aTGG‑TGGG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCa  <  1:630089/100‑1 (MQ=255)
                                                                                                     |                                                                                               
TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑GG‑GGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCA  >  NC_000913/3459189‑3459382

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