Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,459,293 (G)5→6 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.G97.8% 314.6 / 0.3 89coding (950/1482 nt)gspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (41/46);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑GGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGC  >  NC_000913/3459189‑3459381
                                                                                                     |                                                                                             
tGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg‑t                                                                                                <  1:2679346/100‑1 (MQ=255)
 ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg‑t‑g                                                                                              >  1:738515/1‑100 (MQ=255)
 ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtgg‑t‑g                                                                                              >  1:628344/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑ggg                                                                                            >  2:2764590/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑ggg                                                                                            >  2:145761/1‑100 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCTTCATGG‑T‑ggg                                                                                            >  2:891189/1‑99 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gggg                                                                                           <  1:1296527/100‑1 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGACCTTCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gggg                                                                                           <  1:995373/100‑1 (MQ=255)
      gCATTGGGCAACCGCAGGTGAATACCCGTGTGGATATGTCTTTCTCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑ggggg                                                                                          <  2:1479473/99‑1 (MQ=255)
      gCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑gggggg                                                                                         <  2:898504/100‑2 (MQ=255)
           gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGCG‑GGGGGAAAtt                                                                                     >  2:288854/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTc                                                                                    >  2:698954/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTc                                                                                    >  2:706698/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTc                                                                                    >  2:708835/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑ggggggcaattc                                                                                   >  2:349187/1‑93 (MQ=255)
             gCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCg                                                                                   >  2:1361219/1‑100 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGt                                                                                  <  1:1753609/100‑1 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGt                                                                                  <  2:151383/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTg                                                                                 <  1:386316/100‑1 (MQ=255)
                aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTg                                                                                 <  2:163309/99‑1 (MQ=255)
                     aGGTGAATCCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAg                                                                           <  1:1029628/100‑1 (MQ=255)
                       gTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGaa                                                                         <  1:2040030/100‑1 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAAc                                                                       <  1:2219556/100‑1 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAAc                                                                       <  1:1991364/100‑1 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAAc                                                                       >  2:1816018/1‑100 (MQ=255)
                            tACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACcgc                                                                    >  1:3154697/1‑100 (MQ=255)
                             aCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACcgcg                                                                   <  1:2451884/100‑1 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  >  2:2051925/1‑100 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  >  1:426002/1‑100 (MQ=255)
                              cccGTGTGGATATGGCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  <  2:2585197/100‑1 (MQ=255)
                              cccGTGTGGATATATCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  <  2:2594122/100‑1 (MQ=255)
                                cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAg                                                                <  1:997260/100‑1 (MQ=255)
                                 gtgtGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGa                                                               <  1:821271/100‑1 (MQ=255)
                                  tgtgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGAt                                                              >  1:2653666/1‑100 (MQ=255)
                                  tgtgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGAt                                                              >  1:625175/1‑100 (MQ=255)
                                    tgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTg                                                            >  2:205051/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGGAATTCGTGATACAGAAACCGCGCAGAt                                                              >  2:81740/1‑97 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           >  2:1364589/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           <  1:2275076/100‑1 (MQ=255)
                                      gATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCg                                                          <  1:1152490/100‑1 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         >  2:2219328/1‑100 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         <  1:1418340/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                        <  1:1125933/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                        >  1:2956264/1‑100 (MQ=255)
                                           gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                                                     >  2:2538290/1‑100 (MQ=255)
                                           gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                                                     >  1:1241492/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  1:2438152/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  2:2567039/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  1:2577833/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTGCGCCAGGACCTGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgn                                                   <  2:506586/100‑2 (MQ=255)
                                              tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                  <  1:2764625/100‑1 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                 >  1:928561/1‑100 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                 <  2:1598533/100‑1 (MQ=255)
                                                 cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                               >  1:2529648/1‑100 (MQ=255)
                                                 cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                               >  1:2966999/1‑100 (MQ=255)
                                                  gctcgcGGCCTGCGCGCCATACATCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCGGATTGCGGTTCAGGCCTCg                                              >  2:1898220/1‑100 (MQ=255)
                                                     cgcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                                           <  1:1476678/100‑1 (MQ=255)
                                                        ggCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAcc                                        <  2:380005/100‑1 (MQ=255)
                                                         gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                       >  1:2988056/1‑100 (MQ=255)
                                                         gCCTGCGCGCCATACTACGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                       <  2:1912763/100‑1 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                      <  2:2172187/99‑1 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc                                     <  1:753545/100‑1 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    >  1:2999539/1‑100 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    >  2:1181270/1‑100 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    >  2:1181077/1‑100 (MQ=255)
                                                             gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa                                   <  1:1361233/100‑1 (MQ=255)
                                                             gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa                                   >  1:1501376/1‑100 (MQ=255)
                                                                cgcCATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCt                                <  1:2684111/100‑1 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTgg                              <  1:144749/99‑1 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                             >  2:570279/1‑100 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                             <  2:3127667/99‑1 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                            <  2:2248217/100‑1 (MQ=255)
                                                                      aCTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAct                          >  2:1528063/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  2:48951/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  2:358883/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGACTGCGGCTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  2:2200484/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      >  1:85442/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:921940/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:114537/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:737569/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:75505/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      >  1:2806262/1‑100 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  1:2786850/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  1:347576/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    >  1:1411008/1‑100 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGG‑T‑GGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAaca                   >  1:3167098/1‑100 (MQ=255)
                                                                             cAGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  1:2786760/99‑1 (MQ=255)
                                                                              aGGACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  2:161029/98‑1 (MQ=255)
                                                                               ggACCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACt                 >  1:1775483/1‑100 (MQ=255)
                                                                                 aCCCGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTcc               <  2:2778572/100‑1 (MQ=255)
                                                                                  cccGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca              >  1:1348098/1‑100 (MQ=255)
                                                                                   ccGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac             <  2:628335/100‑1 (MQ=255)
                                                                                   ccGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac             <  2:2313180/100‑1 (MQ=255)
                                                                                    cGGATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcaca            <  1:2882212/100‑1 (MQ=255)
                                                                                      gATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac           >  1:2454446/1‑99 (MQ=255)
                                                                                      gATGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACt          <  1:2967505/100‑1 (MQ=255)
                                                                                        tGTCGTCATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTaa        <  2:1241478/100‑1 (MQ=255)
                                                                                            gtcATGG‑TGGGGGGAAATTCGTGATACAGCAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg     <  2:738507/99‑1 (MQ=255)
                                                                                              cATGG‑TGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc  <  1:2538321/100‑1 (MQ=255)
                                                                                              cATGG‑TGGGGGGAAAGTCGTGATACAGAAACCGCGCAGAGTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc  >  2:3102816/1‑100 (MQ=255)
                                                                                                     |                                                                                             
TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGG‑T‑GGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGC  >  NC_000913/3459189‑3459381

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