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

TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGC  >  NC_000913/3459189‑3459381
                                                                                                     |                                                                                             
tGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggt                                                                                                <  1:2679346/100‑1 (MQ=255)
 ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                                              >  1:738515/1‑100 (MQ=255)
 ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCAtggtg                                                                                              >  1:628344/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggg                                                                                            >  2:2764590/1‑100 (MQ=255)
   aaGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggg                                                                                            >  2:145761/1‑100 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCTTCATGGTggg                                                                                            >  2:891189/1‑99 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                           <  1:1296527/100‑1 (MQ=255)
    aGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGACCTTCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggg                                                                                           <  1:995373/100‑1 (MQ=255)
      gCATTGGGCAACCGCAGGTGAATACCCGTGTGGATATGTCTTTCTCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggg                                                                                          <  2:1479473/99‑1 (MQ=255)
      gCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTgggggg                                                                                         <  2:898504/100‑2 (MQ=255)
           gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGCGGGGGGAAAtt                                                                                     >  2:288854/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTc                                                                                    >  2:698954/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTc                                                                                    >  2:706698/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTc                                                                                    >  2:708835/1‑100 (MQ=255)
            ggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTggggggcaattc                                                                                   >  2:349187/1‑93 (MQ=255)
             gCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCg                                                                                   >  2:1361219/1‑100 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                  <  1:1753609/100‑1 (MQ=255)
              cAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGt                                                                                  <  2:151383/100‑1 (MQ=255)
               aaaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                 <  1:386316/100‑1 (MQ=255)
                aaCGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTg                                                                                 <  2:163309/99‑1 (MQ=255)
                     aGGTGAATCCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAg                                                                           <  1:1029628/100‑1 (MQ=255)
                       gTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGaa                                                                         <  1:2040030/100‑1 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                       <  1:2219556/100‑1 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                       <  1:1991364/100‑1 (MQ=255)
                         gAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAAc                                                                       >  2:1816018/1‑100 (MQ=255)
                            tACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgc                                                                    >  1:3154697/1‑100 (MQ=255)
                             aCCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACcgcg                                                                   <  1:2451884/100‑1 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  >  2:2051925/1‑100 (MQ=255)
                              cccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  >  1:426002/1‑100 (MQ=255)
                              cccGTGTGGATATGGCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  <  2:2585197/100‑1 (MQ=255)
                              cccGTGTGGATATATCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCgcgc                                                                  <  2:2594122/100‑1 (MQ=255)
                                cGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAg                                                                <  1:997260/100‑1 (MQ=255)
                                 gtgtGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGa                                                               <  1:821271/100‑1 (MQ=255)
                                  tgtgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGAt                                                              >  1:2653666/1‑100 (MQ=255)
                                  tgtgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGAt                                                              >  1:625175/1‑100 (MQ=255)
                                    tgGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTg                                                            >  2:205051/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGGAATTCGTGATACAGAAACCGCGCAGAt                                                              >  2:81740/1‑97 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           >  2:1364589/1‑100 (MQ=255)
                                     ggATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGc                                                           <  1:2275076/100‑1 (MQ=255)
                                      gATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCg                                                          <  1:1152490/100‑1 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         >  2:2219328/1‑100 (MQ=255)
                                       atatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCgg                                                         <  1:1418340/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                        <  1:1125933/100‑1 (MQ=255)
                                        tatGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGt                                                        >  1:2956264/1‑100 (MQ=255)
                                           gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                                                     >  2:2538290/1‑100 (MQ=255)
                                           gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCa                                                     >  1:1241492/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  1:2438152/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  2:2567039/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgg                                                   >  1:2577833/1‑100 (MQ=255)
                                             cTTTCGCTCGCGGCCTGCGCGCCATACTGCGCCAGGACCTGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAgn                                                   <  2:506586/100‑2 (MQ=255)
                                              tttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGc                                                  <  1:2764625/100‑1 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                 >  1:928561/1‑100 (MQ=255)
                                               ttCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGcc                                                 <  2:1598533/100‑1 (MQ=255)
                                                 cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                               >  1:2529648/1‑100 (MQ=255)
                                                 cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTc                                               >  1:2966999/1‑100 (MQ=255)
                                                  gctcgcGGCCTGCGCGCCATACATCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCGGATTGCGGTTCAGGCCTCg                                              >  2:1898220/1‑100 (MQ=255)
                                                     cgcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTc                                           <  1:1476678/100‑1 (MQ=255)
                                                        ggCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCAcc                                        <  2:380005/100‑1 (MQ=255)
                                                         gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                       >  1:2988056/1‑100 (MQ=255)
                                                         gCCTGCGCGCCATACTACGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCg                                       <  2:1912763/100‑1 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCgg                                      <  2:2172187/99‑1 (MQ=255)
                                                           cTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGc                                     <  1:753545/100‑1 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    >  1:2999539/1‑100 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    >  2:1181270/1‑100 (MQ=255)
                                                            tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                    >  2:1181077/1‑100 (MQ=255)
                                                             gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa                                   <  1:1361233/100‑1 (MQ=255)
                                                             gcgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCa                                   >  1:1501376/1‑100 (MQ=255)
                                                                cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCt                                <  1:2684111/100‑1 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTgg                              <  1:144749/99‑1 (MQ=255)
                                                                   cATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                             >  2:570279/1‑100 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGt                             <  2:3127667/99‑1 (MQ=255)
                                                                    aTACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTa                            <  2:2248217/100‑1 (MQ=255)
                                                                      aCTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTAct                          >  2:1528063/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  2:48951/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  2:358883/1‑100 (MQ=255)
                                                                       cTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGACTGCGGCTCAGGCCTCGCTCACCGGCCATCTGGTActc                         >  2:2200484/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      >  1:85442/1‑100 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:921940/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:114537/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:737569/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      <  2:75505/100‑1 (MQ=255)
                                                                          cGCCAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                      >  1:2806262/1‑100 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  1:2786850/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  1:347576/100‑1 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    >  1:1411008/1‑100 (MQ=255)
                                                                            ccAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAaca                   >  1:3167098/1‑100 (MQ=255)
                                                                             cAGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  1:2786760/99‑1 (MQ=255)
                                                                              aGGACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                    <  2:161029/98‑1 (MQ=255)
                                                                               ggACCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACt                 >  1:1775483/1‑100 (MQ=255)
                                                                                 aCCCGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTcc               <  2:2778572/100‑1 (MQ=255)
                                                                                  cccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCca              >  1:1348098/1‑100 (MQ=255)
                                                                                   ccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac             <  2:628335/100‑1 (MQ=255)
                                                                                   ccGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcac             <  2:2313180/100‑1 (MQ=255)
                                                                                    cGGATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCcaca            <  1:2882212/100‑1 (MQ=255)
                                                                                      gATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac           >  1:2454446/1‑99 (MQ=255)
                                                                                      gATGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACt          <  1:2967505/100‑1 (MQ=255)
                                                                                        tGTCGTCATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTaa        <  2:1241478/100‑1 (MQ=255)
                                                                                            gtcATGGTGGGGGGAAATTCGTGATACAGCAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAg     <  2:738507/99‑1 (MQ=255)
                                                                                              cATGGTGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc  <  1:2538321/100‑1 (MQ=255)
                                                                                              cATGGTGGGGGGAAAGTCGTGATACAGAAACCGCGCAGAGTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc  >  2:3102816/1‑100 (MQ=255)
                                                                                                     |                                                                                             
TGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGC  >  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: