Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,512,506 C→A 10.3% G9V (GGT→GTT)  lolA ← outer membrane lipoprotein carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,512,5060CA10.3% 274.5 / 14.0 93G9V (GGT→GTT) lolAouter membrane lipoprotein carrier protein
Reads supporting (aligned to +/- strand):  ref base C (40/43);  new base A (6/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGACTTTTCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGC  >  NC_002947/4512359‑4512654
                                                                                                                                                   |                                                                                                                                                    
tGCGACTTTTCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCg                                                                                                                                                    >  1:508249/1‑150 (MQ=255)
     cTTTTCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGgc                                                                                                                                               <  1:553919/150‑1 (MQ=255)
          ccAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcag                                                                                                                                          <  1:358875/150‑1 (MQ=255)
          ccAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag                                                                                                                                          <  2:1837003/150‑1 (MQ=255)
           cAGCAACTGGGTCAGGAGTTGTACGTACAGCAAACCGTTATAATTTGTTGCCGATCGTAGGTTCAGTGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcagg                                                                                                                                         >  2:626735/1‑150 (MQ=255)
            aGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcaggc                                                                                                                                        >  1:1851814/1‑150 (MQ=255)
             gCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcaggca                                                                                                                                       <  1:1361899/150‑1 (MQ=255)
               aaCTGGGTCCGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACCGCATGCGAATCGCGCGCATGGGAATCGTCATGGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcaggcagg                                                                                                                                     <  2:1530975/150‑1 (MQ=255)
                aCTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCAACCGGTGGCGCcaggccttt                                                                                                                                    <  2:1851814/150‑5 (MQ=255)
                     gTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtc                                                                                                                               >  1:222609/1‑150 (MQ=255)
                     gTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtc                                                                                                                               <  1:1420822/150‑1 (MQ=255)
                     gTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtc                                                                                                                               <  1:1396500/150‑1 (MQ=255)
                      tCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAACAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtca                                                                                                                              >  2:1145049/1‑150 (MQ=255)
                      tCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtc                                                                                                                               >  1:1439634/1‑149 (MQ=255)
                      tCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtc                                                                                                                               <  2:1439634/149‑1 (MQ=255)
                       cAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtcat                                                                                                                             <  1:1352571/150‑1 (MQ=255)
                       cAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtcat                                                                                                                             <  2:198827/150‑1 (MQ=255)
                        aGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGTGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag                                                                                                                                          >  2:1829556/1‑136 (MQ=255)
                        aGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGTGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag                                                                                                                                          <  1:1829556/136‑1 (MQ=255)
                         ggCGTTGTGCGTCTGGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCTCCGTTAGCTCTAGGCAGGGTGTtcatca                                                                                                                           >  1:121579/1‑150 (MQ=255)
                                  cGTCCTGCTCTCCGACATAACCAGTAACCGAGCAAAGGGCGATGGCAGAGTCGAGCAGCTGGCGAATCTCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCgg                                                                                                                  <  2:981755/150‑1 (MQ=255)
                                   gTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCggg                                                                                                                 <  2:122342/150‑1 (MQ=255)
                                     cTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcag                                                                                                                                          <  2:1396266/123‑1 (MQ=255)
                                     cTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcag                                                                                                                                          >  1:1396266/1‑123 (MQ=255)
                                         cTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGccc                                                                                                           >  2:1888295/1‑150 (MQ=255)
                                             ccGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGNGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCgngg                                                                                                       <  1:1505445/150‑1 (MQ=255)
                                             ccGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCgggg                                                                                                       <  2:365981/150‑1 (MQ=255)
                                                gCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTa                                                                                                    <  1:2306342/150‑1 (MQ=255)
                                                  aTAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCTCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtcatca                                                                                                                           >  2:320323/1‑125 (MQ=255)
                                                  aTAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCTCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtcatca                                                                                                                           <  1:320323/125‑1 (MQ=255)
                                                  aTAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATc                                                                                                  <  2:508249/150‑1 (MQ=255)
                                                  aTAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag                                                                                                                                          >  2:2262717/1‑110 (MQ=255)
                                                  aTAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag                                                                                                                                          <  1:2262717/110‑1 (MQ=255)
                                                      gCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACtt                                                                                              >  2:1283866/1‑150 (MQ=255)
                                                      gCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACtt                                                                                              >  1:1584852/1‑150 (MQ=255)
                                                       cAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCACGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTc                                                                                             <  1:309475/150‑1 (MQ=255)
                                                            aCCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCgg                                                                                                                  <  2:195924/124‑1 (MQ=255)
                                                            aCCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCgg                                                                                                                  >  1:195924/1‑124 (MQ=255)
                                                            aCCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCa                                                                                                 <  1:2403374/141‑1 (MQ=255)
                                                            aCCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCa                                                                                                 >  2:2403374/1‑141 (MQ=255)
                                                             ccGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGATGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCGGCCCGGGGCTATCACGTCCCGCGa                                                                                       >  2:1309415/1‑150 (MQ=255)
                                                             ccGAGACCAGGGTCAGGGCAGAAACCAACAGCGTGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCCTTGGCGCCAGGTAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGa                                                                                       <  2:1029857/150‑1 (MQ=255)
                                                              cGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGTGAAGGCCTGGCGCAAAGGCCTGTGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCGCGTCCCGCGAg                                                                                      >  1:768126/1‑150 (MQ=255)
                                                                    cAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCggg                                                                                                        <  1:762594/126‑1 (MQ=255)
                                                                    cAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCggg                                                                                                        >  2:762594/1‑126 (MQ=255)
                                                                    cAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTg                                                                                >  1:2313508/1‑150 (MQ=255)
                                                                    cAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCCGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTCCCCGCGAGCCGTTg                                                                                >  2:347020/1‑150 (MQ=255)
                                                                           aGGGCAGAAACCAACATCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCa                                                                         >  2:1529911/1‑150 (MQ=255)
                                                                             ggCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATc                                                                       >  2:432688/1‑150 (MQ=255)
                                                                                aGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGTCTATCACTTCCCGCGAGCCGTTGCTGTTCATCggg                                                                    <  1:1145049/150‑1 (MQ=255)
                                                                                 gAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCTTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCgggg                                                                   <  1:1321062/150‑1 (MQ=255)
                                                                                  aaaCCAACAGCATGTGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGt                                                                  >  2:1751936/1‑150 (MQ=255)
                                                                                    aCCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGa                                                                <  2:2032223/150‑1 (MQ=255)
                                                                                      cAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACa                                                              >  1:632334/1‑150 (MQ=255)
                                                                                            cATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCg                                                        >  1:2352343/1‑150 (MQ=255)
                                                                                                   aTCGCGCGCATGGGAATCCTCATTGAGCGTTGATCAAGTGGTGTGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCAtct                                                 <  2:2018602/150‑1 (MQ=255)
                                                                                                           cATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGAc                                         <  1:1529911/150‑1 (MQ=255)
                                                                                                            aTGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACt                                        <  1:1888295/150‑1 (MQ=255)
                                                                                                              gggAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtcatc                                                                                                                            <  2:1065762/64‑1 (MQ=255)
                                                                                                              gggAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTtcatc                                                                                                                            >  1:1065762/1‑64 (MQ=255)
                                                                                                                gAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATc                                            >  2:712097/1‑142 (MQ=255)
                                                                                                                gAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATc                                            <  1:712097/142‑1 (MQ=255)
                                                                                                                 aaTCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATc                                   >  1:433243/1‑150 (MQ=255)
                                                                                                                   tCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCgggg                                                                   >  2:1683854/1‑116 (MQ=255)
                                                                                                                   tCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCgggg                                                                   <  1:1683854/116‑1 (MQ=255)
                                                                                                                   tCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCAt                                 <  1:58894/150‑1 (MQ=255)
                                                                                                                     cTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATgc                               >  1:1925877/1‑150 (MQ=255)
                                                                                                                      tCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATgcg                              >  1:2289589/1‑150 (MQ=255)
                                                                                                                      tCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATgcg                              >  2:1681778/1‑150 (MQ=255)
                                                                                                                        aTTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATc                                            >  1:1136797/1‑134 (MQ=255)
                                                                                                                        aTTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATc                                            <  2:1136797/134‑1 (MQ=255)
                                                                                                                            aGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATgc                               <  1:1479162/143‑1 (MQ=255)
                                                                                                                            aGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATgc                               >  2:1479162/1‑143 (MQ=255)
                                                                                                                             gCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGcc                       >  2:1731829/1‑149 (MQ=255)
                                                                                                                              cGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCAc                      >  2:2140691/1‑150 (MQ=255)
                                                                                                                                   aaggtgCCTGTCGCGTCCGGTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTg                 <  2:1715943/145‑1 (MQ=255)
                                                                                                                                    aaaGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGt                >  1:1557954/1‑150 (MQ=255)
                                                                                                                                    aaaGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGt                <  1:1283866/150‑1 (MQ=255)
                                                                                                                                      aGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCTTAGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAg              <  1:626735/150‑1 (MQ=255)
                                                                                                                                       ggCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGc             <  1:1681778/150‑1 (MQ=255)
                                                                                                                                        gCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGcc            <  2:1571506/150‑1 (MQ=255)
                                                                                                                                        gCCTGGCGCCACAGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGcc            >  2:404274/1‑150 (MQ=255)
                                                                                                                                         ccTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCg           <  2:768126/150‑1 (MQ=255)
                                                                                                                                          cTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGa                                                                                       <  1:684743/73‑1 (MQ=255)
                                                                                                                                          cTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGa                                                                                       >  2:684743/1‑73 (MQ=255)
                                                                                                                                          cTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGa          <  2:632334/150‑1 (MQ=255)
                                                                                                                                          cTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGa          >  2:1643232/1‑150 (MQ=255)
                                                                                                                                          cTGGCGCCAACGGTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGtt                                                                                 <  2:1772956/79‑1 (MQ=39)
                                                                                                                                          cTGGCGCCAACGGTGGCGCCAGGCAGGGTGTTAATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGtt                                                                                 >  1:1772956/1‑79 (MQ=38)
                                                                                                                                           tGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGAt         >  2:1360642/1‑150 (MQ=255)
                                                                                                                                            ggCGCCAACGGTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGCAGCCGATc        >  2:1310056/1‑150 (MQ=255)
                                                                                                                                                 cACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAg   >  1:1340032/1‑150 (MQ=255)
                                                                                                                                                  aCCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGTCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGc  >  1:1232990/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCGACTTTTCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGC  >  NC_002947/4512359‑4512654

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