Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,331,103 G→A 5.9% A169V (GCC→GTC)  PP_4692 ← aspartate/tyrosine/aromatic aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,331,1030GA 5.9% 302.0 / 6.1 102A169V (GCC→GTC) PP_4692aspartate/tyrosine/aromatic aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (48/48);  new base A (3/3);  total (51/51)
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.06e-01

Only 100 of 103 total aligned reads displayed.

CGCTCGGCGCGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCC  >  NC_002947/5330954‑5331246
                                                                                                                                                     |                                                                                                                                               
cGCTCGGCGCGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAggg                                                                                                                                                 <  1:1130652/150‑1 (MQ=255)
 gCTCGGCGCGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCTCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGc                                                                                                                                                >  1:1052806/1‑150 (MQ=255)
      gcgcGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCGACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCTACTAGGGCACCGa                                                                                                                                           >  2:1180470/1‑150 (MQ=255)
        gcgTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACg                                                                                                                                         <  1:1189161/150‑1 (MQ=255)
            ccATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGgtgt                                                                                                                                     <  2:231869/150‑1 (MQ=255)
             cATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTc                                                                                                                                    >  2:597632/1‑150 (MQ=255)
             cATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTc                                                                                                                                    >  1:1678339/1‑150 (MQ=255)
             cATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGGCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGGGTc                                                                                                                                    >  1:853165/1‑150 (MQ=255)
             cATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTGCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACAAGGACACCGTCGgtgcc                                                                                                                                    >  1:728379/1‑148 (MQ=255)
              aTGCCATGGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                   <  2:1099972/150‑1 (MQ=255)
              aTGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                   <  1:1267103/150‑1 (MQ=255)
              aTGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                   >  1:567315/1‑150 (MQ=255)
              aTGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                   >  1:1873618/1‑150 (MQ=255)
               tGCCATAGGTTAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGACACCGTCGGTGCCCt                                                                                                                                  >  1:546776/1‑150 (MQ=255)
                gCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCtt                                                                                                                                 <  2:2360776/150‑1 (MQ=255)
                  cATAGGTTAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGACACCGTCGGTGCCCTGGt                                                                                                                               <  2:546776/150‑1 (MQ=255)
                  cATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGt                                                                                                                               <  1:754990/150‑1 (MQ=255)
                  cATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGt                                                                                                                               <  2:153270/150‑1 (MQ=255)
                  cATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACAAGGACACCGTCGGTGCCCTGGt                                                                                                                               <  2:2265071/150‑1 (MQ=255)
                   aTAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTc                                                                                                                              >  1:1090717/1‑150 (MQ=255)
                      ggTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCa                                                                                                                           >  1:1511809/1‑150 (MQ=255)
                       gTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAg                                                                                                                          >  2:480576/1‑150 (MQ=255)
                           ggCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAgc                                                                                                                      <  2:673939/150‑1 (MQ=255)
                             ccGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAgcgc                                                                                                                    >  1:1435117/1‑150 (MQ=255)
                                tGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCa                                                                                                                 >  2:1726679/1‑150 (MQ=255)
                                  gTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAAc                                                                                                               >  2:1049541/1‑150 (MQ=255)
                                  gTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCAACGACGGTGGCCTTGTCCCAGTAgcgccccac                                                                                                               >  2:16983/1‑145 (MQ=255)
                                   tAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAAcc                                                                                                              >  1:586497/1‑150 (MQ=255)
                                   tAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAAcc                                                                                                              <  1:1525339/150‑1 (MQ=255)
                                     gATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGACGCGACCAGTGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAg                                                                                                            >  2:1103051/1‑150 (MQ=255)
                                     gATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAg                                                                                                            <  1:904299/150‑1 (MQ=255)
                                     gATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGACACCGTCGGTGCCCtggt                                                                                                                               >  2:118671/1‑128 (MQ=255)
                                     gATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGACACCGTCGGTGCCCTGGt                                                                                                                               <  1:118671/131‑1 (MQ=255)
                                       tctcGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGt                                                                                                          >  2:1458397/1‑150 (MQ=255)
                                         tcGTCCACCACTAGCTGGGTATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCg                                                                                                        >  2:2434075/1‑150 (MQ=255)
                                         tcGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAg                                                                                                            >  2:132527/1‑146 (MQ=255)
                                         tcGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAg                                                                                                            <  1:132527/146‑1 (MQ=255)
                                           gTCCACCACCAGAGGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGGCGGGTTGGCCGACGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGc                                                                                                      <  2:632635/150‑1 (MQ=255)
                                               accaccaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGAGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTc                                                                                                  >  1:377961/1‑150 (MQ=255)
                                                 caccaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAg                                                                                                >  2:2124364/1‑150 (MQ=255)
                                                  accaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTc                                                                                                         <  1:427171/140‑1 (MQ=255)
                                                  accaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTc                                                                                                         >  2:427171/1‑140 (MQ=255)
                                                  accaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGc                                                                                               <  2:1873618/150‑1 (MQ=255)
                                                     aGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGa                                                                                            >  2:561325/1‑150 (MQ=255)
                                                      gATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGAt                                                                                           <  1:1726679/150‑1 (MQ=255)
                                                      gATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGAt                                                                                           <  1:1043737/150‑1 (MQ=255)
                                                           ccATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGAACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTc                                                                                      <  2:1238916/150‑1 (MQ=255)
                                                           ccATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGAGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTc                                                                                      <  2:377961/150‑1 (MQ=255)
                                                             aTGTCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTcac                                                                                    <  1:2434075/150‑1 (MQ=255)
                                                               gCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACAt                                                                                  >  1:1392048/1‑150 (MQ=255)
                                                               gCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACAt                                                                                  >  2:2172381/1‑150 (MQ=255)
                                                                ccGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATc                                                                                 >  1:592301/1‑150 (MQ=255)
                                                                      gtgggtggTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGcc                                                                           <  2:1293310/150‑1 (MQ=255)
                                                                       tgggtggTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGccc                                                                          >  1:2071117/1‑150 (MQ=255)
                                                                         ggtggtTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCTAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCAc                                                                        >  1:1800302/1‑150 (MQ=255)
                                                                          gtggtTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCAcc                                                                       >  1:814840/1‑150 (MQ=255)
                                                                              ttCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATc                                                                                 >  2:801289/1‑136 (MQ=255)
                                                                              ttCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATc                                                                                 <  1:801289/136‑1 (MQ=255)
                                                                                   gCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCAcc                                                                       <  1:1686515/141‑1 (MQ=255)
                                                                                    ccAGGTATGCCAACGCCTCGGGCCCCAGAACGGGACGGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTCGCGCTCAACCAGGGCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAAtt                                                             <  2:853165/150‑1 (MQ=255)
                                                                                     caggcaggCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTg                                                            <  2:1252681/150‑1 (MQ=255)
                                                                                        gcaggcCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAAc                                                                                                               >  2:2182827/1‑96 (MQ=255)
                                                                                        gcaggcCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAAc                                                                                                               <  1:2182827/96‑1 (MQ=255)
                                                                                              cAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                   >  1:951305/1‑70 (MQ=255)
                                                                                              cAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGcc                                                   >  2:1038820/1‑150 (MQ=255)
                                                                                              cAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGcc                                                   >  1:209361/1‑150 (MQ=255)
                                                                                              cAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGcc                                                   >  1:1361096/1‑150 (MQ=255)
                                                                                               aaCGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGGACCGACGGTGACCTTGTCCCAGTAGCGCTCAAGCAGGTCGGGAGTCAGCTGATAGTTCACATGCGGGCCCACCGGTAATAATTGGGCCCGGcgc                                                  >  2:429929/1‑148 (MQ=255)
                                                                                               aaCGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCGCATCCGGGCCCACCGGTACCAATTGGGCCCCGccc                                                  >  2:1688576/1‑150 (MQ=255)
                                                                                                 cGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTc                                                <  2:2235075/150‑1 (MQ=255)
                                                                                                  gCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCg                                               >  2:759231/1‑150 (MQ=255)
                                                                                                   ccTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGa                                              >  1:980690/1‑150 (MQ=255)
                                                                                                    cTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCCCATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGAc                                             <  2:586497/150‑1 (MQ=255)
                                                                                                    cTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCCGCATTAATTGGGCCCCGCCCTCGAc                                             >  2:1136248/1‑150 (MQ=255)
                                                                                                     tCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGa                                                                                            >  1:2195042/1‑102 (MQ=255)
                                                                                                     tCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGa                                                                                            <  2:2195042/102‑1 (MQ=255)
                                                                                                          ccccAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAgg                                                                                                           >  2:2104992/1‑82 (MQ=255)
                                                                                                          ccccAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAgg                                                                                                           <  1:2104992/82‑1 (MQ=255)
                                                                                                             cAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGAACAGGGGAAAGACGGTGTCCTTGTCCCAGGAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTCCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAg                                    <  2:314125/150‑1 (MQ=255)
                                                                                                                   ggTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTg                              <  2:1755501/150‑1 (MQ=255)
                                                                                                                   ggTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGGCGGCAGGCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTg                              <  2:1090717/150‑1 (MQ=255)
                                                                                                                    gTACCGGTCGGGTTGGCCGGCGAAGAACCGAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTgg                             >  1:1882015/1‑150 (MQ=255)
                                                                                                                          gTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGNTTg                       <  2:1511809/150‑1 (MQ=255)
                                                                                                                              ggTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACg                   <  2:1392048/150‑1 (MQ=255)
                                                                                                                              ggTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACg                   <  2:1435117/150‑1 (MQ=255)
                                                                                                                               gTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACgg                  <  2:940041/150‑1 (MQ=255)
                                                                                                                               gTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACgg                  >  2:1010131/1‑150 (MQ=255)
                                                                                                                                      ggCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGNCCCCGCNCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAAncc           <  1:1688576/150‑1 (MQ=255)
                                                                                                                                        cGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCgg         <  1:597632/150‑1 (MQ=255)
                                                                                                                                         gAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCggg                                                                             <  2:389173/81‑1 (MQ=255)
                                                                                                                                         gAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCggg                                                                             >  1:389173/1‑81 (MQ=255)
                                                                                                                                         gAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGTAGTGGCGATTGCACGGGTAACCCGGa        >  1:2341617/1‑150 (MQ=255)
                                                                                                                                         gAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGa        <  1:561325/150‑1 (MQ=255)
                                                                                                                                          aaGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGAt       <  1:2172381/150‑1 (MQ=255)
                                                                                                                                           aGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATc      <  2:1052806/150‑1 (MQ=255)
                                                                                                                                            gCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCg     <  1:1180470/150‑1 (MQ=255)
                                                                                                                                              gACCAGTGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGc   <  1:1103051/150‑1 (MQ=255)
                                                                                                                                              gACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGc   <  2:367760/150‑1 (MQ=255)
                                                                                                                                              gACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGc   <  2:1361096/150‑1 (MQ=255)
                                                                                                                                               aCCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGcc  >  2:2086860/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGCTCGGCGCGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCC  >  NC_002947/5330954‑5331246

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: