Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,437 G→A 6.2% H109Y (CAC→TAC)  PP_1324 ← endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,509,4370GA 6.2% 180.8 / 2.7 64H109Y (CAC→TAC) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base A (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

ATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAC  >  NC_002947/1509288‑1509585
                                                                                                                                                       |                                                                                                                                                    
aTTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                      <  1:1430869/150‑1 (MQ=255)
    cATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTa                                                                                                                                                  >  1:1842231/1‑150 (MQ=255)
     aTGTCCATCTGGTGACCTTAAATGGGTCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGCTACATAGAGCAAACAAAACGACGTGGGCTTAGAAT‑‑GACGGGAATTTTGAGATGCATTGAAGGGGTCTGCCTGCATGGTGGCTAc                                                                                                                                                 >  2:1390996/1‑150 (MQ=255)
        tcaaTCTGGTAACTTTAAAGGGGGCGGGGGCCTTTGCTACGCCAAACAACGACAAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCt                                                                                                                                              <  2:1348337/147‑1 (MQ=255)
           aTCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCa                                                                                                                                           >  1:73757/1‑150 (MQ=255)
           aTCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCa                                                                                                                                           <  2:1704748/150‑1 (MQ=255)
            tCTGGTGACCTTAAATGGGGCGGTAGCCTTTGCTACGCGCAACAACGCCCAGCCCGCGATCTAGAGAAGGTAGAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAg                                                                                                                                          <  1:526714/150‑1 (MQ=255)
                  gACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                                    >  1:817182/1‑150 (MQ=255)
                  gACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                                    <  1:1673331/150‑1 (MQ=255)
                    ccTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTGGTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCa                                                                                                                                  >  1:1166615/1‑150 (MQ=255)
                     cTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCAc                                                                                                                                 >  1:838037/1‑150 (MQ=255)
                      ttAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACg                                                                                                                                >  2:682621/1‑150 (MQ=255)
                                     gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGt                                                                                                                >  1:1374234/1‑150 (MQ=255)
                                       cTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTa                                                                                                               >  1:1503795/1‑150 (MQ=255)
                                           gCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTg                                                                                                           <  2:1533365/150‑1 (MQ=255)
                                                gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                    >  1:1710486/1‑136 (MQ=255)
                                                gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                    <  2:1710486/136‑1 (MQ=255)
                                                 cgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAg                                                                                                     <  2:1503795/150‑1 (MQ=255)
                                                  acaacaaCGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                    <  1:677028/149‑1 (MQ=255)
                                                    aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                  >  1:418229/1‑150 (MQ=255)
                                                     acaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTggg                                                                                                 <  1:315448/150‑1 (MQ=255)
                                                      caacCACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                <  1:1380936/150‑1 (MQ=255)
                                                       aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcc                                                                                               >  1:1378239/1‑150 (MQ=255)
                                                       aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcc                                                                                               >  1:751833/1‑150 (MQ=255)
                                                          gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcctcc                                                                                            <  2:1842231/150‑1 (MQ=255)
                                                             cAGCCGGCGATACAAAGGGCAAAAAGGCCGGAGTGGTTTCAGCGTTGGAAGGGGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTc                                                                                         <  1:1884195/150‑1 (MQ=255)
                                                             cAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTc                                                                                         <  1:192497/150‑1 (MQ=255)
                                                                              agCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggc                                                                        <  1:48423/150‑1 (MQ=255)
                                                                                cAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcgg                                                                      >  2:773571/1‑150 (MQ=255)
                                                                                  aaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                    <  2:1049666/102‑1 (MQ=255)
                                                                                  aaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                    >  1:1049666/1‑102 (MQ=255)
                                                                                      aaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCa                                                                >  2:1428991/1‑150 (MQ=255)
                                                                                      aaaCGAAGTGGGTTCAGCATTCGAAGACGTTTTTCAG‑‑CCATTGAAGCGCTCTGCCTGCATGGTGGCTACCCTGCGGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGTGCTCCTACTGCAAGAACAGGTTGGCGGCaacct                                                                >  2:251822/1‑147 (MQ=255)
                                                                                             gTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCg                                                         >  1:1520182/1‑150 (MQ=255)
                                                                                              tGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGt                                                        >  1:377646/1‑150 (MQ=255)
                                                                                                ggTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCCAGTACCAACCGttt                                                      <  2:73757/150‑1 (MQ=255)
                                                                                                  ttCAGCATTCGAAGGCGTTTTTCTT‑‑CTTTTTAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTTGCGGCAAGCACCAACCGTTTCt                                                    >  2:97056/1‑150 (MQ=255)
                                                                                                      gCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCCGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCtt                                                <  1:413315/150‑1 (MQ=255)
                                                                                                       cATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGc                                                                                                          >  2:584275/1‑91 (MQ=255)
                                                                                                       cATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGc                                                                                                          <  1:584275/91‑1 (MQ=255)
                                                                                                        aTTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCACTGACTTGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCATTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTgc                                              <  1:1477900/150‑1 (MQ=255)
                                                                                                         ttCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTgcg                                             >  1:293254/1‑150 (MQ=255)
                                                                                                             aaGGCGTTTTTCAG‑‑CCCTTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCg                                         >  2:163155/1‑150 (MQ=255)
                                                                                                                gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                                   >  1:276997/1‑121 (MQ=255)
                                                                                                                gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                                   <  2:276997/121‑1 (MQ=255)
                                                                                                                gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAgg                                                                                                                                         >  2:1463641/1‑51 (MQ=18)
                                                                                                                gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAgg                                                                                                                                         <  1:1463641/51‑1 (MQ=18)
                                                                                                                  gTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGAt                                    >  1:996169/1‑150 (MQ=255)
                                                                                                                  gTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGAt                                    >  1:69368/1‑150 (MQ=255)
                                                                                                                   tttttCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATa                                   <  2:1378239/150‑1 (MQ=255)
                                                                                                                        cAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTc                              <  1:103002/150‑1 (MQ=255)
                                                                                                                         aG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAAcc                                                          <  1:1714539/121‑1 (MQ=255)
                                                                                                                         aG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAAcc                                                          >  2:1714539/1‑121 (MQ=255)
                                                                                                                          g‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                            <  2:1374234/150‑1 (MQ=255)
                                                                                                                               aTTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGagag                         <  2:293254/150‑1 (MQ=255)
                                                                                                                                ttGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgag                           <  1:803289/147‑1 (MQ=255)
                                                                                                                                ttGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgag                           >  2:803289/1‑147 (MQ=255)
                                                                                                                                ttGAAGCGGTCTGCCTGCATGGAGGCTCAGTTGCAGGACGAGCACGTCGAGGTGGCAGGGCTAGTTGCCCGAGTGGGCCTCCTTCTGTAGCAACAGGGTGGGGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCAATACTGGCgcgcgc                        >  2:1730502/1‑145 (MQ=255)
                                                                                                                                     gCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTc                   >  1:635281/1‑150 (MQ=255)
                                                                                                                                        gtcgccctgcAGGGTAGCCACCATGCAGGCAGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGACAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGa                >  2:1700465/6‑150 (MQ=255)
                                                                                                                                        gTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGa                >  1:816553/1‑150 (MQ=255)
                                                                                                                                         tCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCTGTCGATACTGTCGAGAGCGCCTCCGaa               <  1:773571/150‑1 (MQ=255)
                                                                                                                                             cctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGCCGGCAAGCACCAACGGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGccc           >  1:33467/1‑150 (MQ=255)
                                                                                                                                              ctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGATAGCGCCTCCGAAGCCcg          >  1:409559/1‑150 (MQ=255)
                                                                                                                                                 cATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGt       <  2:838037/150‑1 (MQ=255)
                                                                                                                                                      tggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAc  <  1:682621/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAC  >  NC_002947/1509288‑1509585

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