Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 428,785 G→A 7.2% S75S (AGC→AGT PP_0351 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947428,7850GA 7.2% 197.5 / 2.1 68S75S (AGC→AGTPP_0351transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (26/37);  new base A (3/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAGGTGCGTAATTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCG  >  NC_002947/428639‑428930
                                                                                                                                                  |                                                                                                                                                   
aaGGTGCGTAATTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTg                                                                                                                                                  <  2:278437/150‑1 (MQ=255)
  ggTGCGTAATTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGgc                                                                                                                                                <  2:328702/150‑1 (MQ=255)
      cGTAATTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGgcgcgc                                                                                                                                            <  1:591587/150‑1 (MQ=255)
         aaTTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGAt                                                                                                                                         <  1:31301/150‑1 (MQ=255)
          aTTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGAtt                                                                                                                                        <  2:823135/150‑1 (MQ=255)
          aTTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGCTGCCAGAGCAGGGGCGCCGCCCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGAtt                                                                                                                                        <  2:519392/150‑1 (MQ=255)
             cACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTgca                                                                                                                                    <  2:88577/150‑1 (MQ=255)
              aCCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGACAACCTTGCAAATGCGTCGCCGCAGACAGGCTCGCGATCgcag                                                                                                                                   >  2:1374206/1‑150 (MQ=255)
                ccGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTgcagca                                                                                                                                 >  2:211932/1‑150 (MQ=255)
                ccGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTgcagca                                                                                                                                 <  2:1376578/150‑1 (MQ=255)
                ccGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTgcagca                                                                                                                                 <  2:117891/150‑1 (MQ=255)
                       ggTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTg                                                                                                                          <  2:459698/150‑1 (MQ=255)
                            gcTGGTCGGCCAGCCAGTAGCTGGGTACCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTg                                                                                                                    <  1:1037172/150‑1 (MQ=255)
                                gTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCGGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCgg                                                                                                                <  1:677649/150‑1 (MQ=255)
                                     ccagccagTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAg                                                                                                           >  1:1068720/1‑150 (MQ=255)
                                      cagccagTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGc                                                                                                          >  2:204249/1‑150 (MQ=255)
                                      cagccagTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGc                                                                                                          >  1:538049/1‑150 (MQ=255)
                                            gTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCcgc                                                                                                    <  2:992781/150‑1 (MQ=255)
                                                cTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTa                                                                                                >  1:786452/1‑150 (MQ=255)
                                                        cAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTgcagc                                                                                                                                  >  2:1320766/1‑109 (MQ=255)
                                                        cAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTgcagc                                                                                                                                  <  1:1320766/109‑1 (MQ=255)
                                                        cAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTg                                                                                        <  2:1160650/150‑1 (MQ=255)
                                                        cAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTg                                                                                        <  2:538049/150‑1 (MQ=255)
                                                         aGCGACGGGGCGCCGACCCATGCCAGCATCGCGCCGGGCAGAAACGGGCGGGCCAGCCCGCCGCCAACCTTGGGAATGCGTGGCCGCAGACTGGCGCGCGATTGTAGACGCTCGTTTCGGGCCAAGCCAGCTGCCGCGCTAACGCGTTgg                                                                                       >  2:691295/1‑150 (MQ=255)
                                                         aGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTCGCGGCGGCGGATACCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTgg                                                                                       >  2:37000/1‑150 (MQ=255)
                                                         aGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGACTGGCGCGCGATTGCAGCAGCGCGCGTCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTgg                                                                                       >  2:691238/1‑150 (MQ=255)
                                                              cGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCAt                                                                                  <  2:1068720/150‑1 (MQ=255)
                                                                    gcCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGcc                                                                            <  2:786452/150‑1 (MQ=255)
                                                                      cGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGa                                                                          <  1:1374206/150‑1 (MQ=255)
                                                                          ccATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGGCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTg                                                                      >  2:126360/1‑150 (MQ=255)
                                                                          ccATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTg                                                                      >  2:412997/1‑150 (MQ=255)
                                                                          ccATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTg                                                                      >  1:464351/1‑150 (MQ=255)
                                                                            aTGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGTCGCCTACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGgc                                                                    >  1:330966/1‑150 (MQ=255)
                                                                               ccAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGcc                                                                 >  2:52819/1‑150 (MQ=255)
                                                                                cagcagCGCGCCGAGCAGAAACGTACCATATAGGCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCa                                                                >  1:398239/1‑150 (MQ=255)
                                                                                              gCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTg                                                                                                                                       >  2:807838/1‑67 (MQ=255)
                                                                                              gCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTg                                                                                                                                       <  1:807838/67‑1 (MQ=255)
                                                                                                  aaaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTc                                              >  1:294777/1‑150 (MQ=255)
                                                                                                  aaaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTc                                              <  1:204249/150‑1 (MQ=255)
                                                                                                   aaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCg                                             <  1:665253/150‑1 (MQ=255)
                                                                                                       tGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTg                                                                                                                          <  1:1389523/70‑1 (MQ=255)
                                                                                                       tGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTg                                                                                                                          >  2:1389523/1‑70 (MQ=255)
                                                                                                         ccGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTgg                                         <  2:1270500/148‑1 (MQ=255)
                                                                                                         ccGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTgg                                         >  1:1270500/1‑148 (MQ=255)
                                                                                                            gCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCa                                    >  1:2772/1‑150 (MQ=255)
                                                                                                             ccAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCAt                                   <  1:437230/150‑1 (MQ=255)
                                                                                                             ccAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCAt                                   <  2:464351/150‑1 (MQ=255)
                                                                                                               tggccggcgCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGc                                 <  2:922511/147‑1 (MQ=255)
                                                                                                                  ccgccgCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAg                              >  1:974417/1‑150 (MQ=255)
                                                                                                                   cgccgcCAACCTTGCCAATGCGTCGCCGCTGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAgc                             <  2:330966/150‑1 (MQ=255)
                                                                                                                   cgccgcCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGt                                                                                          <  2:870075/89‑1 (MQ=255)
                                                                                                                   cgccgcCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGt                                                                                          >  1:870075/1‑89 (MQ=255)
                                                                                                                                gCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGTTCGCCATGCGCTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATACGtt                <  2:760064/150‑1 (MQ=255)
                                                                                                                                 ccAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTg               <  1:561054/150‑1 (MQ=255)
                                                                                                                                 ccAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTAGCGGGCGGCATGCCAGCGCATGAATTCGTTg               <  1:37000/150‑1 (MQ=255)
                                                                                                                                 ccAATGCGTCGCCGCAGACTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGAGTGCGCTGCGCTGGCGGCATGCCAGCGCATGAATTCGTTg               <  1:691238/150‑1 (MQ=255)
                                                                                                                                 ccAATGCGTCGCCGCAGACTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGGACGCGGAACTCGCTGGCGGCATGCCAGCGCAGGAATTCGTTg               <  1:691295/150‑1 (MQ=255)
                                                                                                                                  cAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGc              >  1:1008889/1‑150 (MQ=255)
                                                                                                                                   aaTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCg             >  1:57482/1‑150 (MQ=255)
                                                                                                                                      gCGTCGCCGCAGGCTGGCGCGCGATTGCAGGGGGCGTGCGGTTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCCTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTtc          <  2:825331/150‑1 (MQ=255)
                                                                                                                                         tCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGc              >  2:389140/1‑143 (MQ=255)
                                                                                                                                         tCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGc              <  1:389140/143‑1 (MQ=255)
                                                                                                                                         tCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTtcatc       >  1:1231367/1‑150 (MQ=255)
                                                                                                                                           gccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCt                                           >  1:917715/1‑112 (MQ=255)
                                                                                                                                           gccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCt                                           <  2:917715/112‑1 (MQ=255)
                                                                                                                                           gccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAg     <  1:412997/150‑1 (MQ=255)
                                                                                                                                            ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAgg    <  2:294777/150‑1 (MQ=255)
                                                                                                                                              gcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCg  >  2:224084/1‑150 (MQ=255)
                                                                                                                                                gggCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCa                                                            >  1:92826/2‑90 (MQ=255)
                                                                                                                                                gggCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCa                                                            <  2:92826/89‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGGTGCGTAATTCACCCGTGGCGGTGCGCTGGTCGGCCAGCCAGTAGCTGGGTGCCAGCGACGGGGCGCCGACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCTGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCG  >  NC_002947/428639‑428930

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