Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,012 A→C 7.9% intergenic (+49/‑24) rpoC → / → PP_mr08 DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947546,0120AC 7.9% 223.0 / 8.6 76intergenic (+49/‑24)rpoC/PP_mr08DNA‑directed RNA polymerase subunit beta'/rpsL_psuedo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base C (3/3);  total (37/39)
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.22e-01

GCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGA  >  NC_002947/545869‑546155
                                                                                                                                               |                                                                                                                                               
gcAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTg                                                                                                                                           <  1:675511/150‑1 (MQ=255)
 cAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAATTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGc                                                                                                                                          <  1:554104/150‑1 (MQ=255)
  aaGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGTCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGcc                                                                                                                                         >  1:1256780/1‑150 (MQ=255)
    gCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCga                                                                                                                                       <  2:652489/150‑1 (MQ=255)
    gCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCga                                                                                                                                       >  1:652489/1‑150 (MQ=255)
     cGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgat                                                                                                                                      <  1:993357/150‑1 (MQ=255)
      gTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgatga                                                                                                                                    >  2:152471/1‑150 (MQ=255)
       tCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCgatga                                                                                                                                    >  1:845272/1‑150 (MQ=255)
         gCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCAACATTTCATGTCGCCGATGATc                                                                                                                                  <  1:770618/150‑1 (MQ=255)
          ccGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCGAGGCCCCGGCCCTCCCCGCCCGATCATCGGCAACATTTCATGTCGCCGATGATcg                                                                                                                                 >  2:238851/1‑150 (MQ=18)
               aTGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATGCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggcgg                                                                                                                            <  2:42569/150‑1 (MQ=255)
                tGCCGACAAACGGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGCAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggcggg                                                                                                                           >  1:652655/1‑150 (MQ=255)
                tGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggcggg                                                                                                                           <  1:203944/150‑1 (MQ=255)
                 gCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCgggg                                                                                                                          <  1:1187842/150‑1 (MQ=255)
                  ccGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGTCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGa                                                                                                                         <  2:886104/150‑1 (MQ=255)
                   cGACAAACCGCTGCGTGTCAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAg                                                                                                                        <  2:329497/150‑1 (MQ=255)
                            gCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                               >  1:640669/1‑150 (MQ=255)
                               gCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTc                                                                                                            <  2:1256780/150‑1 (MQ=255)
                                 gtgtAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGt                                                                                                          <  2:1350716/150‑1 (MQ=255)
                                    tAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCtt                                                                                                       <  2:1343954/150‑1 (MQ=255)
                                       gcgcCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATCCCAATGGCGGTGGAGTACATGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGAc                                                                                                    <  2:1106863/150‑1 (MQ=255)
                                          ccAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                               <  2:105500/136‑1 (MQ=255)
                                          ccAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGc                                                                                                               >  1:105500/1‑136 (MQ=255)
                                          ccAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                 <  1:1368069/150‑1 (MQ=255)
                                          ccAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                 <  2:538619/150‑1 (MQ=255)
                                                   tGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                 <  1:1259546/141‑1 (MQ=255)
                                                   tGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                 >  2:1259546/1‑141 (MQ=255)
                                                     gAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGaaa                                                                                                                                                >  1:1264081/1‑92 (MQ=255)
                                                     gAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGaaa                                                                                                                                                <  2:1264081/92‑1 (MQ=255)
                                                     gAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggcgg                                                                                                                            <  1:1252768/112‑1 (MQ=255)
                                                     gAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATcgggcgg                                                                                                                            >  2:1252768/1‑112 (MQ=255)
                                                        gCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGAt                                                                                        <  1:1144772/145‑1 (MQ=255)
                                                        gCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGAt                                                                                        >  2:1144772/1‑145 (MQ=255)
                                                        gCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCctct                                                                                   <  2:845272/150‑1 (MQ=255)
                                                           gCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTa                                                                                >  2:883644/1‑150 (MQ=255)
                                                            cACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAg                                                                               >  1:283159/1‑150 (MQ=255)
                                                                   cGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                 >  1:1203058/1‑125 (MQ=255)
                                                                   cGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTgg                                                                                                 <  2:1203058/125‑1 (MQ=255)
                                                                   cGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTg                                                                        >  2:959196/1‑150 (MQ=255)
                                                                   cGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTg                                                                        >  2:211389/1‑150 (MQ=255)
                                                                                 ccAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtttt                                                                         <  2:1244612/135‑1 (MQ=255)
                                                                                 ccAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACtttt                                                                         >  1:1244612/1‑135 (MQ=255)
                                                                                  cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCaa                                                                                           >  2:1150871/1‑116 (MQ=255)
                                                                                  cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCaa                                                                                           <  1:1150871/116‑1 (MQ=255)
                                                                                  cAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTgg                                                         <  1:1281893/150‑1 (MQ=255)
                                                                                  cAGCGGTAATTAAGTACAGGGCAAGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGt                                                        <  1:152471/150‑1 (MQ=255)
                                                                                           ttAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTAc                                                                     >  2:759288/1‑129 (MQ=255)
                                                                                           ttAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTAc                                                                     <  1:759288/129‑1 (MQ=255)
                                                                                            tAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCg                                               >  2:785199/1‑150 (MQ=255)
                                                                                             aaGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGt                                              >  1:283298/1‑150 (MQ=255)
                                                                                             aaGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGt                                              <  2:652655/150‑1 (MQ=255)
                                                                                                   aGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGcc                                        >  2:622462/1‑150 (MQ=255)
                                                                                                    gggCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTTTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCa                                       >  2:183233/1‑150 (MQ=255)
                                                                                                     ggCAAGGCCCCGGCAAGCCTCGGGCGGTCATCGGCGACATGAAATGTTGCCGATGATcgggcgg                                                                                                                            >  1:216042/1‑64 (MQ=255)
                                                                                                     ggCAAGGCCCCGGCAAGCCTCGGGCGGTCATCGGCGACATGAAATGTTGCCGATGATcgggcgg                                                                                                                            <  2:216042/64‑1 (MQ=255)
                                                                                                     ggCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCaa                                      >  2:1287357/1‑150 (MQ=255)
                                                                                                     ggCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCaa                                      >  1:150667/1‑150 (MQ=255)
                                                                                                            ccccGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTgg                               >  1:283172/1‑150 (MQ=255)
                                                                                                              ccGGCAAGACTACGGCGATCAGCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCt                             >  2:1059179/1‑150 (MQ=255)
                                                                                                                       cTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCaa                    >  2:223817/1‑150 (MQ=255)
                                                                                                                        ccgcccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAg                   >  1:33289/6‑150 (MQ=38)
                                                                                                                          gggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTAc                                                                     >  1:551972/1‑98 (MQ=255)
                                                                                                                          gggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTAc                                                                     <  2:551972/98‑1 (MQ=255)
                                                                                                                          gggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGAc                 >  1:172064/1‑150 (MQ=255)
                                                                                                                           ggCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACa                <  2:283298/150‑1 (MQ=255)
                                                                                                                            ccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGtctc                                           <  1:730161/121‑1 (MQ=38)
                                                                                                                            ccGATCATCGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGtctc                                           >  2:730161/2‑122 (MQ=38)
                                                                                                                              gATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATa             <  2:283172/150‑1 (MQ=255)
                                                                                                                               atcatcGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATag            <  2:150667/150‑1 (MQ=255)
                                                                                                                                tcatcGGCAACATTTCATGTCGCCGATGATCGGGCGGGGAGGGCCGGGGCCt                                                                                                             <  1:833026/52‑1 (MQ=11)
                                                                                                                                         aCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTTTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGa  <  1:183233/150‑1 (MQ=255)
                                                                                                                                         aCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCtt                        <  2:103987/128‑1 (MQ=255)
                                                                                                                                         aCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCtt                        >  1:103987/1‑128 (MQ=255)
                                                                                                                                         aCATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGa  >  1:1090070/1‑150 (MQ=255)
                                                                                                                                           aTGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTAAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTg                       >  2:1240707/1‑127 (MQ=255)
                                                                                                                                           aTGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTAAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTg                       >  2:529534/1‑127 (MQ=255)
                                                                                                                                           aTGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTAAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTg                       <  1:1240707/127‑1 (MQ=255)
                                                                                                                                           aTGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTAAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTg                       <  1:529534/127‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GCAAGCGTCGCCGTGATGCCGACAAACCGCTGCGTGTAAGCGCCAGTGAAGTGGAAGCTGCACTGACCGAAGCGCTGAATTCCAGCGGTAATTAAGTACAGGGCAAGGCCCCGGCAAGCCTCGGGCGATCATCGGCGACATGAAATGTTGCCGATGATCGGGCGGGGAGGGCCGGGGCCTCGTCTTGACTGGGTGCAAGATCCTCTTTAGACTTTTGTACCCTTAAATTTGGTGAGGCTCCGTCTCGCCAATTTTTGGCTTTCTTGCAAGACAATAGAGTCGCAAGA  >  NC_002947/545869‑546155

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: