Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,507,063 C→A 5.9% P92P (CCC→CCA cycA → D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,507,0630CA 5.9% 279.7 / 2.7 91P92P (CCC→CCAcycAD‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base A (2/4);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGGCTGCTG  >  NC_002947/5506916‑5507210
                                                                                                                                                   |                                                                                                                                                    
cATCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTg                                                                                                                                                    <  2:44525/150‑1 (MQ=255)
cATCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCCTCACGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTg                                                                                                                                                    <  1:1131692/150‑1 (MQ=255)
 aTCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTgg                                                                                                                                                   >  1:1874438/1‑150 (MQ=255)
 aTCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTgg                                                                                                                                                   >  1:375421/1‑150 (MQ=255)
   cAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGc                                                                                                                                                 >  1:2252561/1‑150 (MQ=255)
    ccccTGGGCGCGCCATCGATCATCTTCATCGACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGcc                                                                                                                                                <  1:1705667/148‑1 (MQ=255)
       cTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCTCCGACCTGCTCGGCCCAGGGGCCgag                                                                                                                                             <  1:187308/150‑3 (MQ=255)
        tGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgagc                                                                                                                                            <  2:2151401/150‑4 (MQ=255)
        tGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTCTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGcc                                                                                                                                                >  2:2348648/1‑146 (MQ=255)
        tGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTCTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGcc                                                                                                                                                <  1:2348648/146‑1 (MQ=255)
         ggCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgagca                                                                                                                                           <  1:242328/150‑5 (MQ=255)
               gCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTAc                                                                                                                                     >  2:1756379/1‑150 (MQ=255)
               gCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTAc                                                                                                                                     <  2:2049501/150‑1 (MQ=255)
                ccATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTAct                                                                                                                                    <  1:1366660/150‑1 (MQ=255)
                   tCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAAAACCACCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggc                                                                                                                                 <  1:1956963/150‑1 (MQ=255)
                     gATCATCTTCGTCTCCATGAGTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggctg                                                                                                                               >  2:101795/1‑150 (MQ=255)
                     gATCATCTTCGTCTACATGATTATCGGCTTCATTCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTActggctg                                                                                                                               >  2:1022282/1‑150 (MQ=255)
                            ttCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTGTCCGCCGACCTGCTCGTCCCCTGGGCCGGCTACTTTACTGGCTGGACCTAc                                                                                                                        >  2:885227/1‑150 (MQ=255)
                                   aCATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTct                                                                                                                >  1:386969/1‑150 (MQ=255)
                                       gATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTctgctg                                                                                                            <  2:1721677/150‑1 (MQ=255)
                                          tATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggt                                                                                                         <  2:1073894/150‑1 (MQ=255)
                                            tCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggtgg                                                                                                       >  1:1407420/1‑150 (MQ=255)
                                               gCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgagc                                                                                                                                            <  1:976135/111‑4 (MQ=255)
                                               gCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCAGGGGCCgagc                                                                                                                                            >  2:976135/1‑108 (MQ=255)
                                                    aTGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGc                                                                                               >  1:453325/1‑150 (MQ=255)
                                                    aTGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGAATTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGc                                                                                               >  1:266694/1‑150 (MQ=255)
                                                      gCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAt                                                                                             <  2:2234690/150‑1 (MQ=255)
                                                       cTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCAtt                                                                                            <  2:608361/150‑1 (MQ=255)
                                                        tGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTg                                                                                           >  2:43224/1‑150 (MQ=255)
                                                          ttcttcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGcc                                                                                         >  1:1529096/1‑150 (MQ=255)
                                                             ttcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGAc                                                                                      <  1:318593/150‑1 (MQ=255)
                                                              tcGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACg                                                                                     <  2:1309276/150‑1 (MQ=255)
                                                                 tCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAc                                                                                                                            <  1:1693602/109‑1 (MQ=255)
                                                                 tCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGAc                                                                                                                            >  2:1693602/1‑109 (MQ=255)
                                                                   aTGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTc                                                                                <  1:778814/150‑1 (MQ=255)
                                                                    tGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTcg                                                                               <  2:1660225/150‑1 (MQ=255)
                                                                    tGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTcg                                                                               >  1:1660225/1‑150 (MQ=255)
                                                                     gcgcgcCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTggg                                                                                                                                                  >  2:1449894/1‑83 (MQ=255)
                                                                     gcgcgcCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTggg                                                                                                                                                  <  1:1449894/83‑1 (MQ=255)
                                                                       gcgcCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGa                                                                            >  1:873354/1‑150 (MQ=255)
                                                                                     ctgctgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACAcc                                                              <  1:1475796/150‑1 (MQ=255)
                                                                                         tgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggt                                                                                                         >  2:2358640/1‑103 (MQ=255)
                                                                                         tgctgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGggt                                                                                                         <  1:2358640/103‑1 (MQ=255)
                                                                                            tgTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCt                                                       <  1:1311772/150‑1 (MQ=255)
                                                                                             gTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTg                                                      <  2:1529096/150‑1 (MQ=255)
                                                                                              tCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTgg                                                     <  2:375421/150‑1 (MQ=255)
                                                                                                gAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGtt                                                   >  1:1445896/1‑150 (MQ=255)
                                                                                                 aaCCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTc                                                  <  1:1754089/150‑1 (MQ=255)
                                                                                                 aaCCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTc                                                  >  1:1106136/1‑150 (MQ=255)
                                                                                                  aCCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcc                                                 >  1:763646/1‑150 (MQ=255)
                                                                                                    cTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGACGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTcccc                                               <  2:1874438/150‑1 (MQ=255)
                                                                                                     tGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCa                                                                                                    >  1:2186465/1‑96 (MQ=255)
                                                                                                     tGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCa                                                                                                    <  2:2186465/96‑1 (MQ=255)
                                                                                                       aaCTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGAc                                            <  1:1124158/150‑1 (MQ=255)
                                                                                                           aCAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGc                                        >  1:27348/1‑150 (MQ=255)
                                                                                                            cAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgc                                                                        <  1:2294445/117‑1 (MQ=255)
                                                                                                            cAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgc                                                                        >  2:2294445/1‑117 (MQ=255)
                                                                                                            cAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGCCCTGcc                                       >  2:399437/1‑150 (MQ=255)
                                                                                                              aGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgccgc                                                                     <  2:1901133/118‑1 (MQ=255)
                                                                                                              aGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgccgc                                                                     >  1:1901133/1‑118 (MQ=255)
                                                                                                                 cGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGt                                  <  1:1756379/150‑1 (MQ=255)
                                                                                                                 cGTTCATCGAATTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGt                                  <  2:266694/150‑1 (MQ=255)
                                                                                                                  gTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTCTTCCGATCTTGACGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTg                                 <  1:1022282/150‑1 (MQ=255)
                                                                                                                  gTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGAGGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTg                                 >  1:2302537/1‑150 (MQ=255)
                                                                                                                       tCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTg                                 >  2:1582940/1‑145 (MQ=255)
                                                                                                                       tCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTg                                 <  1:1582940/145‑1 (MQ=255)
                                                                                                                       tCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATAc                            >  2:536195/1‑150 (MQ=255)
                                                                                                                       tCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATAc                            >  1:222763/1‑150 (MQ=255)
                                                                                                                        cGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATAcc                           >  1:296100/1‑150 (MQ=255)
                                                                                                                         gACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCg                          >  1:985838/1‑150 (MQ=255)
                                                                                                                         gACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCg                          <  2:296100/150‑1 (MQ=255)
                                                                                                                          aCTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCgg                         >  1:397163/1‑150 (MQ=255)
                                                                                                                            ttttCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgcg                       >  2:751668/1‑150 (MQ=255)
                                                                                                                             tttCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgcgc                      >  1:123514/1‑150 (MQ=255)
                                                                                                                              ttCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCt                     <  2:1445896/150‑1 (MQ=255)
                                                                                                                                ccgccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGa                   <  2:453325/150‑1 (MQ=255)
                                                                                                                                  gccgACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGATGGATTTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGAcc                 >  1:561619/1‑150 (MQ=255)
                                                                                                                                     gACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgcc                                                                       <  2:141484/93‑1 (MQ=255)
                                                                                                                                     gACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgcc                                                                       >  1:141484/1‑93 (MQ=255)
                                                                                                                                     gACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGc              <  1:751668/150‑1 (MQ=255)
                                                                                                                                     gACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGc              <  1:101795/150‑1 (MQ=255)
                                                                                                                                     ctccTGCTCGGCCCCTGGGCCGGCTACTTTACTGGGTTGTTTTATTTGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGc              <  2:2371556/148‑1 (MQ=255)
                                                                                                                                     atcttGCTCGGCCACTGGGCGGGCTCATCTACCGGCGCACCCGACTGGGGTTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACCCCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGc              <  2:265161/146‑1 (MQ=255)
                                                                                                                                       ccTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGTATACCGGCGCTGACCTGCGt            <  2:27348/150‑1 (MQ=255)
                                                                                                                                        cTGCTCGGCCCCTGGGCCGGCTACTTTACTGGAAAATCCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCgag           >  1:2371556/1‑148 (MQ=255)
                                                                                                                                          gCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgcgc                      <  1:514075/137‑1 (MQ=255)
                                                                                                                                          gCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGgcgc                      >  2:514075/1‑137 (MQ=255)
                                                                                                                                          gCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCgagat         >  1:265161/1‑146 (MQ=255)
                                                                                                                                           gccGGCCCAGGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAggg        >  1:1806126/3‑150 (MQ=255)
                                                                                                                                           cTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAggg        <  1:536195/150‑1 (MQ=255)
                                                                                                                                           cTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAggg        <  1:2228791/150‑1 (MQ=255)
                                                                                                                                                 cccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgccgc                                                                     >  1:2098594/1‑83 (MQ=255)
                                                                                                                                                 cccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATcgccgc                                                                     <  2:2098594/83‑1 (MQ=255)
                                                                                                                                                 cccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGc                                        >  2:2368399/1‑112 (MQ=255)
                                                                                                                                                 cccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGc                                        <  1:2368399/112‑1 (MQ=255)
                                                                                                                                                 cccTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCATTTGATACCGGCGCTGACCTTCGTAGGgctgctg  >  2:997866/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATCAGCCTGGCCGGGCCATCGATCATCTTCGTCTACATGATTATCGGCTTCATGCTGTTCTTCGTCATGCGCGCCATGGGCGAACTGCTGCTGTCGAACCTGAACTACAAGTCGTTCATCGACTTTTCCGCCGACCTGCTCGGCCCCTGGGCCGGCTACTTTACTGGCTGGACCTACTGGTTCTGCTGGGTGGTCACCGGCATTGCCGACGTGGTCGCGATCGCCGCCTACACCCAGTTCTGGTTCCCCGACCTGCCGCAGTGGATACCGGCGCTGACCTGCGTAGGGCTGCTG  >  NC_002947/5506916‑5507210

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