Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,433 T→G 6.8% 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.8% 227.8 / 5.0 83H110P (CAT→CCT) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base G (4/2);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTG  >  NC_002947/1509292‑1509581
                                                                                                                                               |                                                                                                                                                    
cATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGGCTGCGTGGTGGCTa                                                                                                                                              >  1:472219/1‑150 (MQ=255)
   gTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAccc                                                                                                                                           <  2:1801597/150‑1 (MQ=255)
      cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACTCGCAACAACGACCAACCCGCGATACATAGGGCAAAAAAAACGAAGAGATTTCCGCATTTGAAGGCGTGTTTCCGCAAGAGAAGCGGTCTGCCTGAATGGGGGATACCCTgt                                                                                                                                        >  2:1073486/1‑149 (MQ=255)
        tCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAg                                                                                                                                      <  2:84655/150‑1 (MQ=255)
          tGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAggg                                                                                                                                    <  1:986616/150‑1 (MQ=255)
            gTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGGAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCg                                                                                                                                  >  1:2127307/1‑150 (MQ=255)
                 cTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCAc                                                                                                                             <  1:2159771/150‑1 (MQ=255)
                 cTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTTGTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCGGGGCGACCAc                                                                                                                             >  1:908505/1‑150 (MQ=255)
                 cTTAAATGGGGCGGTGGCCTTTGCTACGCGCAAAAAAGCTTTGATGGAGCGGCACAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTGTTAAGCCATGAAAGCGGTCGGGCAGCAGCGGGGCTACCCTGCAGGGCGACCAc                                                                                                                             >  1:1615604/1‑150 (MQ=255)
                  ttAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACg                                                                                                                            >  2:1824265/1‑150 (MQ=255)
                   tAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGGGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGt                                                                                                                           >  2:1751449/1‑150 (MQ=255)
                         gggCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                            <  2:2369313/127‑1 (MQ=255)
                         gggCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                            >  1:2369313/1‑127 (MQ=255)
                           gCGGTGGCCTCTGTTTCGCGTAACAACGACCAGCCCGCGATACATAGAGCAAAACAAACGAAGTAGTGTTAGGTTTCGAAGGCGTTTTTCAGACCTTGAAACGATATCCCTGTCCGTTGGCTACCCTGCAGGGCGACCACGTCGAAGCgg                                                                                                                   <  1:1751449/150‑1 (MQ=255)
                            cGGTGGCATTTGACACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAATAACGAAGTGGGTTCACCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGACGCGGc                                                                                                                  <  1:1073486/150‑1 (MQ=255)
                              gTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                  <  1:1216592/116‑1 (MQ=255)
                              gTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggtg                                                                                                                                                  >  2:1216592/1‑116 (MQ=255)
                               tGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAgg                                                                                                               <  1:471423/150‑1 (MQ=255)
                                 gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAgggg                                                                                                             <  1:1956498/150‑1 (MQ=255)
                                  ccTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGt                                                                                                            >  1:1414149/1‑150 (MQ=255)
                                    tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGTAAGTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAg                                                                                                          <  1:289167/150‑1 (MQ=255)
                                    tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTc                                                                                                                          <  2:337072/134‑1 (MQ=255)
                                    tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTc                                                                                                                          >  1:337072/1‑134 (MQ=255)
                                    tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAACGAAGTGGGTTCAGCATTCGAGGGCGTTTTTCGGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGt                                                                                                         >  2:271507/1‑150 (MQ=255)
                                    tttGCTACGCGCAACAACGACAAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAg                                                                                                          <  2:2345192/150‑1 (MQ=255)
                                        actggggcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACTACGTCGAAGCGGCAGGGGTAGTTGc                                                                                                      >  2:1193377/7‑150 (MQ=255)
                                         tACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggt                                                                                                                                                   <  2:1339411/104‑1 (MQ=255)
                                         tACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAtggt                                                                                                                                                   >  1:1339411/1‑104 (MQ=255)
                                           cgcgCAACAACGACCAGCCCGCGATACATAGAGCAAAAACAACGACGTGGCCTCAGCCTTACCCATGGTGTATCGGGCCTTAAAGCGGTCGGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGcccc                                                                                                   <  1:168795/150‑1 (MQ=255)
                                            gcgaAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGCTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCa                                                                                                  >  2:1053963/1‑150 (MQ=255)
                                              gcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                <  2:1479840/150‑1 (MQ=255)
                                              gcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                <  1:1193377/150‑1 (MQ=255)
                                               caacaaCGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTg                                                                                               <  1:2223919/150‑1 (MQ=255)
                                               aaacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCTAGTg                                                                                               <  2:1897779/149‑1 (MQ=255)
                                                    acGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCa                                                                                                                                       >  2:2087155/1‑105 (MQ=255)
                                                    acGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCa                                                                                                                                       <  1:2087155/105‑1 (MQ=255)
                                                      gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCa                                                                                                  >  2:892753/1‑140 (MQ=255)
                                                      gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCa                                                                                                  <  1:892753/140‑1 (MQ=255)
                                                        ccAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCtt                                                                                      <  1:2007059/150‑1 (MQ=255)
                                                            cccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGc                                                                                  >  2:1844548/1‑150 (MQ=255)
                                                             ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCa                                                                                 <  1:2416290/150‑1 (MQ=255)
                                                                     cATAGAGCAAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAg                                                                          <  1:569759/150‑1 (MQ=255)
                                                                     cATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                         <  1:2212618/150‑1 (MQ=255)
                                                                      aTAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAggg                                                                        >  2:479339/1‑150 (MQ=255)
                                                                      aTAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAggg                                                                        >  2:122023/1‑150 (MQ=255)
                                                                      aTAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCGCCTTCTGCAGGAACAggg                                                                        >  1:121076/1‑150 (MQ=255)
                                                                           gCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcg                                                                   <  1:2057034/150‑1 (MQ=255)
                                                                                aaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                               <  1:1844548/150‑1 (MQ=255)
                                                                                  aaaaCGACGTGGGTTCAGCAAGCGAATGGGGTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGATTCGCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGc                                                             <  1:1732027/150‑1 (MQ=255)
                                                                                       gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGa                                                                                                                                 <  1:809907/77‑1 (MQ=255)
                                                                                       gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAGGGTGGCTACCCTGCAGGGCGa                                                                                                                                 >  2:809907/1‑77 (MQ=255)
                                                                                          gTGGGTTCAGCATTCGAAGGCGTTTTTTTTCTTTTTTAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCg                                                     >  2:485572/1‑150 (MQ=255)
                                                                                          gTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCg                                                     >  2:1644340/1‑150 (MQ=255)
                                                                                               ttCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCt                                                >  2:93109/1‑150 (MQ=255)
                                                                                               ttCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCt                                                >  1:398219/1‑150 (MQ=255)
                                                                                                   gCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCtt                                            >  1:1691450/1‑150 (MQ=255)
                                                                                                    cATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCTCCAGTgg                                                                                              >  2:213468/1‑99 (MQ=255)
                                                                                                    cATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCTCCAGTgg                                                                                              <  1:213468/99‑1 (MQ=255)
                                                                                                           aaGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCg                                     >  1:737960/1‑150 (MQ=255)
                                                                                                              gCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGGGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCg                                  >  2:1679975/1‑150 (MQ=255)
                                                                                                                 tttttCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATa                               <  2:2127307/150‑1 (MQ=255)
                                                                                                                 tttttCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATa                               <  1:1824265/150‑1 (MQ=255)
                                                                                                                       aGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                            <  2:512715/83‑1 (MQ=255)
                                                                                                                       aGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                            >  1:512715/1‑83 (MQ=255)
                                                                                                                        gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGTACAGGGTGGCGGCAAGCACCACACGTTTCTTCTTGCGCCCGTCGATACTGTCga                        >  2:896064/1‑150 (MQ=255)
                                                                                                                        gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                        <  1:271507/150‑1 (MQ=255)
                                                                                                                        gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                        >  1:1964529/1‑150 (MQ=255)
                                                                                                                         ccATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgag                       >  1:1677397/1‑150 (MQ=255)
                                                                                                                          cATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgaga                      >  1:2363469/1‑150 (MQ=255)
                                                                                                                            ttGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                <  2:1899480/74‑1 (MQ=255)
                                                                                                                            ttGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                >  1:1899480/1‑74 (MQ=255)
                                                                                                                               aaGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGcc                 <  1:2419295/150‑1 (MQ=255)
                                                                                                                                    gtcgccctgcAGGGTAGCCACCATGCAGGCAGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                        <  2:2019672/133‑1 (MQ=38)
                                                                                                                                    gtcgccctgcAGGGTAGCCACCATGCAGGCAGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                        >  1:2019672/6‑138 (MQ=38)
                                                                                                                                     tCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGaa           <  1:1422142/150‑1 (MQ=255)
                                                                                                                                         cctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAg                                                                                <  1:1288986/77‑1 (MQ=255)
                                                                                                                                         cctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAg                                                                                >  2:1288986/1‑77 (MQ=255)
                                                                                                                                          ctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGACGCCcg      >  2:1363615/1‑150 (MQ=255)
                                                                                                                                          ctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCcg      >  1:2352442/1‑150 (MQ=255)
                                                                                                                                          ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcg                                                                   >  2:570088/1‑89 (MQ=38)
                                                                                                                                          ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcg                                                                   <  1:570088/89‑1 (MQ=38)
                                                                                                                                           tgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCcgc     <  2:1677397/150‑1 (MQ=255)
                                                                                                                                             cATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTc                                   >  2:680571/1‑118 (MQ=255)
                                                                                                                                             cATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTc                                   <  1:680571/118‑1 (MQ=255)
                                                                                                                                             cATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGt   <  2:398219/150‑1 (MQ=255)
                                                                                                                                             cATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGt   >  2:238252/1‑150 (MQ=255)
                                                                                                                                              aTGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTg  >  2:1635624/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CATGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTG  >  NC_002947/1509292‑1509581

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