Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,331,109 A→T 5.7% V167D (GTC→GAC)  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,1090AT 5.7% 324.4 / 8.2 106V167D (GTC→GAC) PP_4692aspartate/tyrosine/aromatic aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (49/51);  new base T (3/3);  total (52/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 107 total aligned reads displayed.

GCGCGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCCAGCAGCCAG  >  NC_002947/5330960‑5331255
                                                                                                                                                     |                                                                                                                                                  
gcgcGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCGACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCTACTAGGGCACCGa                                                                                                                                                    >  2:1180470/1‑150 (MQ=255)
  gcgTCCATGCCATAGGTGAGGCCGTGTTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACg                                                                                                                                                  <  1:555/150‑1 (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)
       cATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTGCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACAAGGACACCGTCGgtgcc                                                                                                                                             >  1:728379/1‑148 (MQ=255)
        aTGCCATGGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                            <  2:1099972/150‑1 (MQ=255)
        aTGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                            >  1:567315/1‑150 (MQ=255)
        aTGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTcc                                                                                                                                            <  1:1267103/150‑1 (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                                                                                                                                        <  2:153270/150‑1 (MQ=255)
            cATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGt                                                                                                                                        <  1:754990/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                                                                                                                     <  1:132527/146‑1 (MQ=255)
                                   tcGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAg                                                                                                                     >  2:132527/1‑146 (MQ=255)
                                     gTCCACCACCAGAGGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGGCGGGTTGGCCGACGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGc                                                                                                               <  2:632635/150‑1 (MQ=255)
                                           caccaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAg                                                                                                         >  2:2124364/1‑150 (MQ=255)
                                            accaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTc                                                                                                                  >  2:427171/1‑140 (MQ=255)
                                            accaGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTc                                                                                                                  <  1:427171/140‑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)
                                                       aTGTCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTcac                                                                                             <  1:2434075/150‑1 (MQ=255)
                                                         gCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACAt                                                                                           >  2:2172381/1‑150 (MQ=255)
                                                         gCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACAt                                                                                           >  1:1392048/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                                                                                          <  1:801289/136‑1 (MQ=255)
                                                                        ttCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATc                                                                                          >  2:801289/1‑136 (MQ=255)
                                                                             gCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCAcc                                                                                >  2:1686515/1‑141 (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                                                            >  1:209361/1‑150 (MQ=255)
                                                                                        cAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGcc                                                            >  1:1361096/1‑150 (MQ=255)
                                                                                        cAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGcc                                                            >  2:1038820/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                                                                                                                    <  1:2104992/82‑1 (MQ=255)
                                                                                                    ccccAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAgg                                                                                                                    >  2:2104992/1‑82 (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:1435117/150‑1 (MQ=255)
                                                                                                                        ggTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACg                            <  2:1392048/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)
                                                                                                                                                cACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCCAGCAGcc    >  1:629211/1‑150 (MQ=255)
                                                                                                                                                 aCCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCCAGCAGCCa   <  2:1620385/150‑1 (MQ=255)
                                                                                                                                                  ccGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCCAGCAGCCAg  <  1:1038820/150‑1 (MQ=255)
                                                                                                                                                    gACGTTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCa                                                    <  1:2377800/98‑1 (MQ=255)
                                                                                                                                                    gACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCa                                                    >  2:2377800/1‑98 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGCGTCCATGCCATAGGTGAGGCCGTGGTAGATCTCGTCCACCACCAGATGGCCATGCCGCTCGTGGGTGGTTCTGGCCAGGCAGGCCAACGCCTCGCGCCCCAGCACGGTACCGGTCGGGTTGGCCGGCGAAGCGACCAGGGCACCGACGGTGTCCTTGTCCCAGTAGCGCTCAACCAGGTCGGCAGTCAGCTGATAGTTCACATCCGGGCCCACCGGTACCAATTGGGCCCCGCCCTCGACCAGGCGCAGGAAGTGGCGATTGCACGGGTAACCCGGATCGGCCAGCAGCCAG  >  NC_002947/5330960‑5331255

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: