Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,565 A→T 5.2% L352Q (CTG→CAG)  ldcC ← lysine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,677,5650AT 5.2% 297.9 / 4.2 93L352Q (CTG→CAG) ldcClysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (46/42);  new base T (3/2);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGGCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTCGCGGCGATCGCC  >  NC_002947/4677416‑4677714
                                                                                                                                                     |                                                                                                                                                     
aGGCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTgca                                                                                                                                                       <  1:1803581/150‑1 (MQ=255)
aGGCCTGGCTGAACGACGCGCGCCGATTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTgca                                                                                                                                                       <  2:1840117/150‑1 (MQ=255)
 ggCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTgcag                                                                                                                                                      <  1:2110125/150‑1 (MQ=255)
 ggCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTgcag                                                                                                                                                      <  1:2049188/150‑1 (MQ=255)
         tGAACGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGa                                                                                                                                              <  2:571244/150‑1 (MQ=255)
           aaCGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGaca                                                                                                                                            <  1:589543/150‑1 (MQ=255)
              gccgcGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTg                                                                                                                                         <  2:1932171/150‑1 (MQ=255)
               ccgcGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTgg                                                                                                                                        >  2:2058273/1‑150 (MQ=255)
               ccgcGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTgg                                                                                                                                        >  1:57053/1‑150 (MQ=255)
                 gcgAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCa                                                                                                                                      <  1:1509140/150‑1 (MQ=255)
                   gAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCAcc                                                                                                                                    >  2:348743/1‑150 (MQ=255)
                      cagcTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAg                                                                                                                                 >  2:1052/1‑150 (MQ=255)
                      cagcTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCCGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAg                                                                                                                                 <  2:1234592/150‑1 (MQ=255)
                       agcTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAgg                                                                                                                                >  2:2115500/1‑150 (MQ=255)
                        gcTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAggg                                                                                                                               >  2:1943219/1‑150 (MQ=255)
                        gcTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGACGGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAggg                                                                                                                               <  1:21849/150‑1 (MQ=255)
                         cTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACACAGGGc                                                                                                                              >  2:1088013/1‑150 (MQ=255)
                          ttGTGGGTGGAATGGGTACTGAACACCGGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAgg                                                                                                                                <  2:859143/147‑1 (MQ=255)
                          ttGTGGGTGGAATGGGTACTGAACACCGGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAgg                                                                                                                                >  1:859143/1‑147 (MQ=255)
                             tgggtggAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGc                                                                                                                          <  2:111168/150‑1 (MQ=255)
                               ggtggAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCtt                                                                                                                        >  2:725593/1‑150 (MQ=255)
                                gtggAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTg                                                                                                                       >  1:28552/1‑150 (MQ=255)
                                         gTACTGAACACCAGCGAGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGGCTGCTTGAGCAGCCCg                                                                                                              >  2:1223081/1‑150 (MQ=255)
                                           aCTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGc                                                                                                            >  1:1034315/1‑150 (MQ=255)
                                           aCTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCCTAGCGCCCGGTGAAAAAGTCGTGGAACGAGGACTCGGCAAACCAGGATTGGTTATTGGACTGCACCTCGTCACTGTCACCCAGGGCCTGCTTGATCAGCCCGGc                                                                                                            <  1:1863192/150‑1 (MQ=255)
                                             tGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGt                                                                                                          >  1:595267/1‑150 (MQ=255)
                                             tGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGt                                                                                                          >  1:1822194/1‑150 (MQ=255)
                                              gAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTg                                                                                                         <  2:2147559/150‑1 (MQ=255)
                                               aaCACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTgg                                                                                                        <  2:592176/150‑1 (MQ=255)
                                                acacCAGCGGGCTGTCTGCCGCGCCGGCGGGGCCCATGGTGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGTAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGt                                                                                                       <  1:2089485/150‑1 (MQ=255)
                                                 cacCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTa                                                                                                      >  1:331076/1‑150 (MQ=255)
                                                  acCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTaa                                                                                                     >  2:712430/1‑150 (MQ=255)
                                                  acCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTaa                                                                                                     >  2:1045534/1‑150 (MQ=255)
                                                    cAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAaca                                                                                                   <  2:1822194/150‑1 (MQ=255)
                                                         ggCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGAACCTCGACACTTGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGcc                                                                                              >  1:1117538/1‑150 (MQ=255)
                                                            aggcTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGc                                                                                                                          <  1:59487/116‑1 (MQ=255)
                                                            aggcTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGc                                                                                                                          >  2:59487/4‑119 (MQ=255)
                                                             gTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTNGCACCCNGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCcgtc                                                                                          <  1:2414425/150‑1 (MQ=255)
                                                             gTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCcgtc                                                                                          >  1:2080884/1‑150 (MQ=255)
                                                              tCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTg                                                                                                                                         <  1:561198/102‑1 (MQ=255)
                                                              tCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTg                                                                                                                                         >  2:561198/1‑102 (MQ=255)
                                                              tCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCTCGGCGTGGTAACACAGCCcgtcg                                                                                         >  1:1555618/1‑150 (MQ=255)
                                                                tGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTa                                                                                       <  2:331076/150‑1 (MQ=255)
                                                                tGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTa                                                                                       <  2:718245/150‑1 (MQ=255)
                                                                 gCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCTTAg                                                                                      >  1:1298383/1‑150 (MQ=255)
                                                                 gCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAg                                                                                      <  1:712430/150‑1 (MQ=255)
                                                                 gCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAg                                                                                      <  1:1943219/150‑1 (MQ=255)
                                                                 gCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAg                                                                                      <  2:1034315/150‑1 (MQ=255)
                                                                  ccGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAgg                                                                                     >  1:1221545/1‑150 (MQ=255)
                                                                    gcgcAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTg                                                                                   <  2:1117538/150‑1 (MQ=255)
                                                                     cgcAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTgg                                                                                  <  2:57053/150‑1 (MQ=255)
                                                                     cgcAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTgg                                                                                  <  1:2058273/150‑1 (MQ=255)
                                                                      gcAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGa                                                                                 <  1:1698929/150‑1 (MQ=255)
                                                                        aGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGt                                                                               <  1:1064627/150‑1 (MQ=255)
                                                                          gCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTg                                                                             <  2:1530854/150‑1 (MQ=255)
                                                                          gCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTg                                                                             <  1:2341023/150‑1 (MQ=255)
                                                                           cGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTgg                                                                            >  2:628364/1‑150 (MQ=255)
                                                                            ggTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGt                                                                           <  1:628364/150‑1 (MQ=255)
                                                                                 ccATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGACATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTa                                                                      >  1:2222289/1‑150 (MQ=255)
                                                                                  cATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATCGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTAc                                                                     >  1:777767/1‑150 (MQ=255)
                                                                                  cATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTAc                                                                     >  2:53980/1‑150 (MQ=255)
                                                                                   aTGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAggg                                                                                                                               >  1:779001/1‑91 (MQ=255)
                                                                                   aTGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAggg                                                                                                                               <  2:779001/91‑1 (MQ=255)
                                                                                   aTGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTAcc                                                                    >  2:1875183/1‑150 (MQ=255)
                                                                                          aGCGCCCGGTGAAAAAGTCGTGGAACGATTCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCt                                                             <  1:1052/150‑1 (MQ=255)
                                                                                                   tGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCt                                                    >  1:1067916/1‑150 (MQ=255)
                                                                                                     aaaaaGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGc                                                  >  1:1909875/1‑150 (MQ=255)
                                                                                                       aaaGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTaa                                                                                                     >  1:315326/1‑97 (MQ=255)
                                                                                                       aaaGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTaa                                                                                                     <  2:315326/97‑1 (MQ=255)
                                                                                                             gTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATc                                          >  2:130328/1‑150 (MQ=255)
                                                                                                                 aaCGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACgg                                      <  1:358067/150‑1 (MQ=255)
                                                                                                                       gCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAgg                                >  1:1037505/1‑150 (MQ=255)
                                                                                                                       gCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAgg                                <  2:2080884/150‑1 (MQ=255)
                                                                                                                         aTAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAg                                                                                <  1:2284195/100‑1 (MQ=255)
                                                                                                                         aTAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAg                                                                                >  2:2284195/1‑100 (MQ=255)
                                                                                                                           aGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGt                            <  1:1875183/150‑1 (MQ=255)
                                                                                                                            ggCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGtt                           >  2:25418/1‑150 (MQ=255)
                                                                                                                             gCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAaca                                                                                                   >  2:1378143/1‑77 (MQ=255)
                                                                                                                             gCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAaca                                                                                                   <  1:1378143/77‑1 (MQ=255)
                                                                                                                             gCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTg                          >  2:1343836/1‑150 (MQ=255)
                                                                                                                               aaaCCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGc                        >  1:384972/1‑150 (MQ=255)
                                                                                                                               aaaCCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGACAGGGGGTTGGc                        >  1:1686573/1‑150 (MQ=255)
                                                                                                                                       cagtgtcgAGGTGCTGCACTTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCt                >  1:993658/5‑150 (MQ=255)
                                                                                                                                        cTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGGTCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCtt               >  2:308117/1‑150 (MQ=255)
                                                                                                                                        cTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGc                                                     >  2:525723/1‑112 (MQ=255)
                                                                                                                                        cTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGc                                                     <  1:525723/112‑1 (MQ=255)
                                                                                                                                        cTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCtt               >  1:417280/1‑150 (MQ=255)
                                                                                                                                        agtgtcgAGGTGCTGCACTTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATa                                                        >  2:2157534/4‑109 (MQ=255)
                                                                                                                                        agtgtcgAGGTGCTGCACTTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCtt               >  2:2250121/4‑150 (MQ=255)
                                                                                                                                        agtgtcgAGGTGCTGCACTTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTGGGTCACTACCGCCAGCTTGATa                                                        <  1:2157534/106‑1 (MQ=255)
                                                                                                                                         tcgtcgAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTc              <  1:1045534/150‑1 (MQ=255)
                                                                                                                                         tcgtcgAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGCCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTc              >  2:2413356/1‑150 (MQ=255)
                                                                                                                                           gtcgAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTctc            <  2:595267/150‑3 (MQ=255)
                                                                                                                                           gtcgAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTcgc            <  2:1221545/150‑1 (MQ=255)
                                                                                                                                             cgAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTCgcgg          <  1:25418/150‑1 (MQ=255)
                                                                                                                                                aGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGCGTTACTCACTACCGTCCGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTCGCGGCGa       >  1:282385/1‑150 (MQ=255)
                                                                                                                                                  tGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTCGCGGCGATc     <  2:1298383/150‑1 (MQ=255)
                                                                                                                                                     agcaCCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCACGGGGGTTGGCCTGAATCTTCGCGGCGATCGcc  >  1:41474/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCCTGGCTGAACGCCGCGAGCAGCTTGTGGGTGGAATGGGTACTGAACACCAGCGGGCTGTCTGCCGCGCAGGCGGTGCCCATGGCGTAGCGCCCGGTGAAAAAGTCGTGGAACGCGGCATAGGCAAACCAGGCCTCGTCGAAGTGCAGCACCTCGACACTGGCACCCAGGGCCTGCTTGATCAGCCCGGCGTGGTAACACAGCCCGTCGTAGGTGGAGTTGGTCACTACCGCCAGCTTGATACGCTGCCCACGATCACGGGCAAGGGGGTTGGCCTGAATCTTCGCGGCGATCGCC  >  NC_002947/4677416‑4677714

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: ‑