Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,462,943 G→A 5.9% G7D (GGC→GAC)  gspK → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,462,9430GA 5.9% 318.0 / 6.3 102G7D (GGC→GAC) gspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (46/50);  new base A (3/3);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 106 total aligned reads displayed.

CTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAG  >  NC_000913/3462797‑3463088
                                                                                                                                                  |                                                                                                                                                 
cTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                                                                                                                                                 >  2:55787/1‑149 (MQ=255)
cTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                                                                                                                                                 >  2:63055/1‑149 (MQ=255)
 ttCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                                                                                                                                                 >  6:82648/1‑148 (MQ=255)
 ttCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGt                                                                                                                                                >  8:129456/1‑149 (MQ=255)
 ttCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGt                                                                                                                                                >  1:489903/1‑149 (MQ=255)
     tctctcAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCTAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTcgcg                                                                                                                                            >  5:149093/1‑149 (MQ=255)
     tctctcAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTcgcg                                                                                                                                            >  1:478537/1‑149 (MQ=255)
           aGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAGCAGCGCGACGCCGCgttgt                                                                                                                                        <  7:8536/147‑2 (MQ=255)
           aGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCgctgct                                                                                                                                       <  8:1207/148‑1 (MQ=255)
            gTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCa                                                                                                                                     >  8:221647/1‑149 (MQ=255)
             tCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCAt                                                                                                                                    >  3:153524/1‑149 (MQ=255)
             tCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCGCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCAt                                                                                                                                    >  6:148845/1‑149 (MQ=255)
               aCGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCAtt                                                                                                                                   <  6:350315/148‑1 (MQ=255)
               aCGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCAtt                                                                                                                                   <  8:369487/148‑1 (MQ=255)
                cGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGt                                                                                                                                 <  5:82648/149‑1 (MQ=255)
                 gTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCTCGCTGCTCATTGTa                                                                                                                                >  4:322684/1‑149 (MQ=255)
                 gTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTa                                                                                                                                <  2:267471/149‑1 (MQ=255)
                    aCATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACt                                                                                                                              >  3:92277/1‑148 (MQ=255)
                      aTCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCAt                                                                                                                           <  6:59400/149‑1 (MQ=255)
                      aTCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCAt                                                                                                                           <  5:480159/149‑1 (MQ=255)
                      aTCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCAt                                                                                                                           <  5:480113/149‑1 (MQ=255)
                           cAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATg                                                                                                                          <  2:478537/145‑1 (MQ=255)
                              ccTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTgg                                                                                                                   >  5:294959/2‑149 (MQ=255)
                               cTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGgc                                                                                                                  >  5:61122/1‑149 (MQ=255)
                                tCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGgcg                                                                                                                 >  8:410496/1‑149 (MQ=255)
                                tCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGgcg                                                                                                                 >  1:407219/1‑149 (MQ=255)
                                  aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGCCGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCt                                                                                                               >  5:477534/1‑149 (MQ=255)
                                  aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCa                                                                                                                                     >  8:247646/1‑127 (MQ=255)
                                  aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCa                                                                                                                                     <  7:247646/127‑1 (MQ=255)
                                  aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGgcgc                                                                                                                <  8:474368/148‑1 (MQ=255)
                                  aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCt                                                                                                               <  5:148845/149‑1 (MQ=255)
                                    cAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCa                                                                                                             <  2:489903/149‑1 (MQ=255)
                                     aaTATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCAt                                                                                                            <  6:70595/149‑1 (MQ=255)
                                       tatGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGGAGAGTCGCCCGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCGTGTTACTGGCGCTCATgt                                                                                                          >  8:259233/1‑148 (MQ=255)
                                       tatGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGTGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATgg                                                                                                          <  5:57879/149‑1 (MQ=255)
                                       tatGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATgg                                                                                                          >  4:158728/1‑149 (MQ=255)
                                       tatGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATgg                                                                                                          <  4:23800/149‑1 (MQ=255)
                                           gggggTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGgctg                                                                                                       <  2:55312/148‑1 (MQ=255)
                                              ggTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGctct                                                                                                   <  7:449768/149‑1 (MQ=255)
                                                   cAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGgcag                                                                                              <  1:55787/149‑1 (MQ=255)
                                                   cAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGgcag                                                                                              <  1:63055/149‑1 (MQ=255)
                                                     gcgcCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGgcagca                                                                                            <  6:418946/149‑1 (MQ=255)
                                                         cGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACa                                                                                        >  8:335785/1‑149 (MQ=255)
                                                          gTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACAt                                                                                       <  6:61122/149‑1 (MQ=255)
                                                                     ccTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTgg                                                                                                                   >  2:136647/1‑110 (MQ=255)
                                                                     ccTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTgg                                                                                                                   <  1:136647/110‑1 (MQ=255)
                                                                     ccTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtt                                                                            >  3:24149/1‑149 (MQ=255)
                                                                      cTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTgg                                                                                                 <  8:303394/127‑1 (MQ=255)
                                                                      cTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTgg                                                                                                 >  7:303394/1‑127 (MQ=255)
                                                                      cTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGttt                                                                           >  3:300220/1‑149 (MQ=255)
                                                                      cTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGttt                                                                           <  6:149093/149‑1 (MQ=255)
                                                                       tGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCAt                                                                                                            <  2:423583/115‑1 (MQ=255)
                                                                       tGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCAt                                                                                                            >  1:423583/1‑115 (MQ=255)
                                                                       tGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                          >  7:477207/1‑149 (MQ=255)
                                                                         aGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                          <  7:9034/147‑1 (MQ=255)
                                                                         aGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                          <  7:129456/147‑1 (MQ=255)
                                                                           cAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTcaca                                                                      <  7:263892/149‑1 (MQ=255)
                                                                                 gCGAGAGAAGAGTCGCCAGCTCAGACGCACGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCAcagccagc                                                                >  4:425093/1‑149 (MQ=255)
                                                                                   gagagaAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCtt                                                              >  2:439994/1‑149 (MQ=255)
                                                                                     gagaAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAACCAGCtt                                                              >  5:428992/1‑147 (MQ=255)
                                                                                      agaagaGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAg                                                           <  4:300220/149‑1 (MQ=255)
                                                                                       gaagaGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGc                                                          <  3:158728/149‑1 (MQ=255)
                                                                                             tCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAAccc                                                    <  7:221647/149‑1 (MQ=255)
                                                                                                     tCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACACCGCGGCGTCGCGCTGCTCATTGTCCTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGCCATGCCGCTCAGTTTTCACAGCCCGCTTCAGCGAACCCGCCAggg                                             >  6:37083/1‑147 (MQ=255)
                                                                                                     tCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGt                                             >  2:401529/1‑148 (MQ=255)
                                                                                                       aGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTaaa                                          >  1:270845/1‑149 (MQ=255)
                                                                                                       aGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTaaa                                          >  6:410511/1‑149 (MQ=255)
                                                                                                         aCGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGc                                                  <  8:419151/139‑1 (MQ=255)
                                                                                                               gCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCatcatc                                  <  1:401529/149‑1 (MQ=255)
                                                                                                                     aataatAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAAc                            >  1:189744/1‑149 (MQ=255)
                                                                                                                      ataataACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACg                           <  2:189744/149‑1 (MQ=255)
                                                                                                                      ataataACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTAAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACg                           >  2:462388/1‑149 (MQ=255)
                                                                                                                          taaCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCaa                       >  7:313165/1‑149 (MQ=255)
                                                                                                                          taaCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCaa                       <  4:24149/149‑1 (MQ=255)
                                                                                                                            aCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAata                     <  3:425093/149‑1 (MQ=255)
                                                                                                                            aCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAat                      <  6:217644/148‑1 (MQ=255)
                                                                                                                             ccATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAatat                    <  7:335785/149‑1 (MQ=255)
                                                                                                                              cATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAatat                    >  4:179208/1‑148 (MQ=255)
                                                                                                                              cATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATg                   <  2:407219/149‑1 (MQ=255)
                                                                                                                                tGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGat                 <  5:410511/149‑1 (MQ=255)
                                                                                                                                  tacaatGAGCAGCGCGACGCCGCGTTGTTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAg                                               >  5:314026/4‑117 (MQ=38)
                                                                                                                                  tacaatGAGCAGCGCGACGCCGCGTTGTTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAg                                               <  6:314026/114‑1 (MQ=38)
                                                                                                                                  tacaatGAGCAGCGCGACGCCGCGTTGTTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGatat               >  7:108777/4‑149 (MQ=255)
                                                                                                                                  aataatGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAACCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGata                <  6:428992/148‑1 (MQ=255)
                                                                                                                                   ataatGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATAtt              <  4:92277/149‑1 (MQ=255)
                                                                                                                                   ataatGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATAtt              <  6:294959/149‑1 (MQ=255)
                                                                                                                                      atGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGaa           >  1:33080/1‑149 (MQ=255)
                                                                                                                                      atGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGaa           <  8:321257/149‑1 (MQ=255)
                                                                                                                                       tGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGa                  >  2:265516/1‑141 (MQ=255)
                                                                                                                                       tGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGa                  <  1:265516/141‑1 (MQ=255)
                                                                                                                                       tGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAAc          >  3:144046/1‑149 (MQ=255)
                                                                                                                                       tGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAAc          <  6:477534/149‑1 (MQ=255)
                                                                                                                                        gAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACCTGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACt         >  1:409534/1‑149 (MQ=255)
                                                                                                                                        gAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACt         >  6:324765/1‑149 (MQ=255)
                                                                                                                                          acaacGCGGCGTCGCGCTGCTCTTTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTgg       <  3:322684/149‑1 (MQ=255)
                                                                                                                                           cagcgcgACGCCGCGTTGTTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGc      <  1:26076/146‑1 (MQ=255)
                                                                                                                                            aacGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGCCATGCCGCTCAGTGTTCACAGCCAGCTTCCGCGACCCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCg     >  2:366495/1‑149 (MQ=255)
                                                                                                                                            aacGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCg     <  4:153524/149‑1 (MQ=255)
                                                                                                                                             acGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCg     >  2:41424/1‑148 (MQ=255)
                                                                                                                                               gcggcgTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGgag  <  1:439994/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAG  >  NC_000913/3462797‑3463088

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