Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,131,860 C→A 6.5% L72M (CTG→ATG) ‡ PP_4517 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,131,8600CA 6.5% 196.1 / 2.9 61L72M (CTG→ATG) ‡PP_4517hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

GACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGG  >  NC_002947/5131719‑5132006
                                                                                                                                               |                                                                                                                                                  
gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGCCACATCGTGCCCGAGCGCGAAAGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCtggtg                                                                                                                                               <  1:119530/149‑1 (MQ=255)
gACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGgc                                                                                                                                             >  6:59218/1‑149 (MQ=255)
     tGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCAGGGTGgcctc                                                                                                                                          <  6:72432/147‑4 (MQ=255)
         gccgcACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAgg                                                                                                                                     <  5:338100/148‑1 (MQ=255)
          ccgcACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGca                                                                                                                                   >  6:183191/1‑149 (MQ=255)
            gcACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcaca                                                                                                                                 >  3:155053/1‑149 (MQ=255)
              acacGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcacaca                                                                                                                               >  6:120341/1‑149 (MQ=255)
              acacGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGcacaca                                                                                                                               <  6:3639/149‑1 (MQ=255)
                    cTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGccc                                                                                                                         <  7:196840/149‑1 (MQ=255)
                       ccTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGc                                                                                                                       >  3:117489/1‑148 (MQ=255)
                        cTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGc                                                                                                                       >  7:293844/1‑147 (MQ=255)
                        cTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCAGGGTGgcctc                                                                                                                                          >  7:97136/1‑125 (MQ=255)
                        cTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCAGGGTGgcctc                                                                                                                                          <  8:97136/128‑4 (MQ=255)
                                tctTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGc                                                                                                               >  3:7096/1‑147 (MQ=255)
                                    tCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGc                                                                                                           <  7:12281/147‑1 (MQ=255)
                                        tCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACgaa                                                                                                     <  6:172492/149‑1 (MQ=255)
                                        tCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACgaa                                                                                                     <  5:183191/149‑1 (MQ=255)
                                          aaCGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACgaag                                                                                                    <  5:120341/148‑1 (MQ=255)
                                           aCGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACgaagaa                                                                                                  >  7:136518/1‑149 (MQ=255)
                                             gCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAAtt                                                                                                <  3:251586/149‑1 (MQ=255)
                                              ccG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTc                                                                                               >  3:288757/1‑149 (MQ=255)
                                                    ggACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGAcc                                                                                           <  6:293494/149‑1 (MQ=255)
                                                      aagTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGAcctg                                                                                         >  2:328645/3‑149 (MQ=255)
                                                        gTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGAcctgcc                                                                                       >  5:204319/1‑149 (MQ=255)
                                                               aGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgac                                                                                >  8:244398/1‑149 (MQ=255)
                                                                    gAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgacctcga                                                                           <  7:244398/149‑1 (MQ=255)
                                                                     aaTGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgacctcgac                                                                          <  5:59218/149‑1 (MQ=255)
                                                                      aTGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgacctcgacc                                                                         >  4:19297/1‑149 (MQ=255)
                                                                             gACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCg                                                                   <  8:136518/148‑1 (MQ=255)
                                                                                 cccTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAg                                                              >  4:58813/1‑149 (MQ=255)
                                                                                     gTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCg                                                           >  5:25760/1‑148 (MQ=255)
                                                                                      ttCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAACGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcga                                                                                 >  4:207525/1‑125 (MQ=255)
                                                                                      ttCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAACGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcga                                                                                 <  3:207525/125‑1 (MQ=255)
                                                                                        cGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGcc                                                                                                                      >  4:143979/1‑86 (MQ=255)
                                                                                         gCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCACCAGCTGGCCAGCGACgaa                                                                                                     <  8:3726/102‑1 (MQ=255)
                                                                                           cGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGAtt                                                    >  7:106371/1‑149 (MQ=255)
                                                                                           cGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGAtt                                                    >  7:111768/1‑149 (MQ=255)
                                                                                               gagCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGa                                                <  6:204319/149‑1 (MQ=255)
                                                                                                agCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGGGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGAc                                               <  8:266051/149‑1 (MQ=255)
                                                                                                 gCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGAcc                                              >  1:188161/1‑149 (MQ=255)
                                                                                                         aCGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGc                                      >  6:112938/1‑149 (MQ=255)
                                                                                                             cAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCt                                     >  1:216614/1‑146 (MQ=255)
                                                                                                             cAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCt                                     <  2:216614/146‑1 (MQ=255)
                                                                                                             cAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGt                                  <  8:293844/149‑1 (MQ=255)
                                                                                                             cAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGt                                  <  2:188161/149‑1 (MQ=255)
                                                                                                               gCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGc                                <  3:19297/149‑1 (MQ=255)
                                                                                                               gCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGc                                <  3:58813/149‑1 (MQ=255)
                                                                                                                cGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGCCCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCa                               >  5:290211/1‑149 (MQ=255)
                                                                                                                    gCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCGCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTTCATCgg                           <  2:3985/149‑1 (MQ=255)
                                                                                                                         gcgcACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCa                      <  5:10247/149‑1 (MQ=255)
                                                                                                                          cgcACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCAt                     <  6:25760/149‑1 (MQ=255)
                                                                                                                              ccGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGaa                 <  4:117489/149‑1 (MQ=255)
                                                                                                                               cGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAg                >  5:110831/1‑149 (MQ=255)
                                                                                                                               cGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAg                <  1:328645/149‑1 (MQ=255)
                                                                                                                                 agagCGAGGCCACCCTGGTGGCGCTGATGGCACACATTGCCCGCCAGCTGGCTAGCGACGAAGAATTCGACCTGCCCTTCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGc              <  4:7096/149‑1 (MQ=255)
                                                                                                                                 agaTCGAGGCCACCCTGGTGGCGCTCCAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGc              >  6:254842/1‑149 (MQ=255)
                                                                                                                                 agaTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGc              <  4:288757/149‑1 (MQ=255)
                                                                                                                                   aTCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCgt            >  8:179728/1‑149 (MQ=255)
                                                                                                                                       gcgCCACCAGGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTc        >  7:122799/3‑149 (MQ=255)
                                                                                                                                           cACCCTGGTGGCGCTCAAGGCACAAATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGcgaactcgac                                                                          <  7:285863/79‑1 (MQ=255)
                                                                                                                                            aCCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGGCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCg   >  5:305731/1‑149 (MQ=255)
                                                                                                                                              ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCgg  >  7:141604/1‑148 (MQ=255)
                                                                                                                                              ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCAGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCgg  >  3:300661/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCG‑‑AGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGGCCACCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGCTGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGG  >  NC_002947/5131719‑5132006

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