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.

TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCA  >  NC_000913/3459189‑3459382
                                                                                                     |                                                                                               
tGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggt                                                                                                  <  1:3694116/100‑1 (MQ=255)
 ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGACATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                                                >  2:741452/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggg                                                                                             >  1:2158430/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggg                                                                                             >  2:3521524/1‑100 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                            <  2:2233950/100‑1 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                            <  2:1147618/100‑1 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                            <  2:3658002/100‑1 (MQ=255)
     ggCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCTGATGTCGTCATGGTggggg                                                                                           <  1:1162273/100‑1 (MQ=255)
       cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggga                                                                                         >  2:860744/1‑98 (MQ=255)
       cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggga                                                                                         >  2:630093/1‑98 (MQ=255)
       cATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggga                                                                                         >  2:2726096/1‑98 (MQ=255)
        aTTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggggaa                                                                                        >  2:2311706/1‑97 (MQ=255)
        aTTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCGGGCGGCCTGCGCGCCATGCTTCGCCAGGAGCCGGATGTCGTCAGGGTGGGGGgaa                                                                                        >  1:1131070/1‑100 (MQ=255)
           gggCAAACGCAGGTGCATACACGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCTTGGTGGGGGGAAAtt                                                                                      >  1:772911/1‑100 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                   <  1:3363200/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCTTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGGCAGACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                  <  1:357948/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGTGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                  <  2:3774791/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                  <  1:725863/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                  <  1:2358105/100‑1 (MQ=255)
                agcGCAGGTGCACACCCGTGTGGTTATGTCGTTCGCTCGCGGCCGGCGCGCCATACGTCGCCCGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                  <  1:2094537/97‑1 (MQ=255)
                aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGa                                                                                 <  2:2168483/100‑1 (MQ=255)
                 aCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGAt                                                                                <  2:1138618/100‑1 (MQ=255)
                     aGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAg                                                                            <  2:3746396/100‑1 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                           >  1:485655/1‑100 (MQ=255)
                      ggTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGa                                                                           >  1:2379808/1‑100 (MQ=255)
                        tGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGaaa                                                                         >  2:2251084/1‑100 (MQ=255)
                          aaTACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAcc                                                                       >  2:3400675/1‑100 (MQ=255)
                             aCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgcg                                                                    <  1:1215580/100‑1 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                   <  2:199488/100‑1 (MQ=255)
                                cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAg                                                                 >  2:1656614/1‑100 (MQ=255)
                                    tgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGCGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTg                                                             >  1:973749/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                            >  2:126119/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                            >  1:3487910/1‑100 (MQ=255)
                                      gATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATTCAGAAACCGCGCAGATTGCg                                                           <  1:2311440/100‑1 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                          <  2:756639/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                         <  1:246620/100‑1 (MQ=255)
                                         atGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGCAATTCGACATACAGAAACCGCGCAGATTGCGGt                                                         >  1:2182587/1‑99 (MQ=255)
                                            tCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAg                                                     <  1:2744430/100‑1 (MQ=255)
                                             cTTTCGCTCGTGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGTGGGAGATAAGTGATACAGAAACCGCGAGCATCGCGGTTCAgg                                                    >  1:2829092/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAGATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                    >  1:2556630/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                    <  1:142910/100‑1 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                    >  1:1359434/1‑100 (MQ=255)
                                              tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                   <  1:3475167/100‑1 (MQ=255)
                                              tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                   <  2:2529197/100‑1 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                  <  1:545334/100‑1 (MQ=255)
                                               gtcgctcgcGGCCTGCGCGCCATACGTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                  <  1:1968737/99‑1 (MQ=255)
                                                 cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                                >  1:1507298/1‑100 (MQ=255)
                                                  gctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCg                                               >  2:3087557/1‑100 (MQ=255)
                                                      gcgGCCTGCGCGCCATACTTCGCCAGGACCCGAATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                   <  1:1027948/92‑1 (MQ=255)
                                                       cgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAc                                          >  1:2140844/1‑100 (MQ=255)
                                                       cgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAc                                          <  2:2684421/100‑1 (MQ=255)
                                                        ggCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGAAGATTGCGGTTCAGGCCTCGCTCAcc                                         >  1:2734197/1‑100 (MQ=255)
                                                         gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                        <  1:39769/100‑1 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                       <  1:2853719/100‑1 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGGCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                       <  1:825308/100‑1 (MQ=255)
                                                          ccTGCGCGCCATACTTCGCCAGGACCCGGATGCCGTCAAGGTGGGGGGAAAATCGAGGTAAAGAACCGGACAAGACTGCGGTTCACGCCGCGCt                                             >  1:370224/1‑94 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc                                      >  2:2941333/1‑100 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                     >  2:1830445/1‑100 (MQ=255)
                                                                cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGACTCAGGCCTCGCTCACCGGCCATCt                                 >  1:1958483/1‑100 (MQ=255)
                                                                 gcCATACTGCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATc                                  <  1:1280929/98‑1 (MQ=255)
                                                                  ccATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTg                                >  2:2065935/1‑99 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                             <  1:1208377/100‑1 (MQ=255)
                                                                     tACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAc                            <  2:1450362/100‑1 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          >  1:939349/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          <  2:2791740/100‑1 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          >  2:777461/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                          >  1:1662238/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                       <  2:2844652/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                       <  2:2528290/100‑1 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      >  1:1131927/1‑100 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      >  2:65504/1‑100 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      <  2:2389405/100‑1 (MQ=255)
                                                                           gCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCaa                      >  1:729015/1‑100 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCacc                     >  1:646247/1‑98 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                     <  2:3039074/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                     <  1:1507851/100‑1 (MQ=255)
                                                                              aGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAacac                   >  2:3184094/1‑100 (MQ=255)
                                                                              aGGACCCGGATGTCGTCAGGGTGGGGGGAAATTCGTGATACAGACACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAaca                    <  1:2478262/99‑1 (MQ=255)
                                                                                 aCCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTcc                <  1:2036117/100‑1 (MQ=255)
                                                                                 aCCCGGATCTCGTCATGGTGGGGGGAAATTCGTGATTCAGAACCCGCACAGATCTCGCGTCAGGCCTCGCTCACCGACCATCTGGTACTCTCAACACTcc                <  1:2819616/100‑1 (MQ=255)
                                                                                  cccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca               >  2:2249997/1‑100 (MQ=255)
                                                                                  cccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca               >  2:1943692/1‑100 (MQ=255)
                                                                                   ccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac              >  1:2401382/1‑100 (MQ=255)
                                                                                     ggATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            >  2:1415909/1‑100 (MQ=255)
                                                                                     ggATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            <  1:2856172/100‑1 (MQ=255)
                                                                                     ggATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            >  2:1598747/1‑100 (MQ=255)
                                                                                     ggATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac            >  2:3209727/1‑100 (MQ=255)
                                                                                      ggtGTCGTCATGGTGGGGGGAAATTCGTGGTACAGACACCGCGCAGATTGCGGTCCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac              <  1:3136184/95‑1 (MQ=255)
                                                                                      gATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACt           <  1:3184073/100‑1 (MQ=255)
                                                                                         gtcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc        >  1:2393789/1‑100 (MQ=255)
                                                                                         gtcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc        <  1:1065222/100‑1 (MQ=255)
                                                                                         gtcgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc        <  1:385916/100‑1 (MQ=255)
                                                                                           cgtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg      >  2:1470828/1‑100 (MQ=255)
                                                                                            gtcATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGt     >  1:3796085/1‑100 (MQ=255)
                                                                                              cATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTg    <  1:1170668/99‑1 (MQ=255)
                                                                                              cATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc   <  1:1737874/100‑1 (MQ=255)
                                                                                               aTGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTg    <  1:2301396/98‑1 (MQ=255)
                                                                                               aTGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc   <  1:3674575/99‑1 (MQ=255)
                                                                                               aTGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCa  <  2:772916/100‑1 (MQ=255)
                                                                                               aTGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCa  <  1:630089/100‑1 (MQ=255)
                                                                                                     |                                                                                               
TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCA  >  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: