Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,433 T→G 6.3% H110P (CAT→CCT)  PP_1324 ← endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,509,4330TG 6.3% 177.0 / 4.8 64H110P (CAT→CCT) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base G (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 = 3.65e-01

GATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGT  >  NC_002947/1509287‑1509580
                                                                                                                                                    |                                                                                                                                                   
gATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAtggt                                                                                                                                                  <  1:1216516/150‑1 (MQ=255)
 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:1673331/150‑1 (MQ=255)
                   gACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                               >  1:817182/1‑150 (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)
                                      gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGC‑AAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑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                                                                                                               <  2:1710486/136‑1 (MQ=255)
                                                 gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                               >  1:1710486/1‑136 (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:751833/1‑150 (MQ=255)
                                                        aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcc                                                                                          >  1:1378239/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                                                              <  2:276997/121‑1 (MQ=255)
                                                                                                                 gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                              >  1:276997/1‑121 (MQ=255)
                                                                                                                 gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAgg                                                                                                                                    <  1:1463641/51‑1 (MQ=18)
                                                                                                                 gCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAgg                                                                                                                                    >  2:1463641/1‑51 (MQ=18)
                                                                                                                   gTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGAt                               >  1:69368/1‑150 (MQ=255)
                                                                                                                   gTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGAt                               >  1:996169/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)
                                                                                                                                                    |                                                                                                                                                   
GATTGCATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAG‑‑CCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGT  >  NC_002947/1509287‑1509580

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