Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,319,210 C→A 6.4% S193S (TCG→TCT trpC ← indole‑3‑glycerolphosphate synthetase and N‑(5‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,319,2100CA 6.4% 256.4 / 4.0 86S193S (TCG→TCTtrpCindole‑3‑glycerolphosphate synthetase and N‑(5‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base A (3/3);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCCATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGGTCGTC  >  NC_000913/1319061‑1319359
                                                                                                                                                     |                                                                                                                                                     
tGGGCCATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATc                                                                                                                                                       <  1:1533639/150‑1 (MQ=255)
  ggCCATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGa                                                                                                                                                     <  2:726129/150‑1 (MQ=255)
   gCCATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGAc                                                                                                                                                    >  1:771463/1‑150 (MQ=255)
   gCCATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGAc                                                                                                                                                    >  1:961929/1‑150 (MQ=255)
     cATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACaa                                                                                                                                                  >  1:746281/1‑150 (MQ=255)
      aTCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACaaa                                                                                                                                                 >  2:1375609/1‑150 (MQ=255)
          aCGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCAc                                                                                                                                             >  2:1488078/1‑150 (MQ=255)
            gCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGc                                                                                                                                           >  1:342846/1‑150 (MQ=255)
             ccGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCa                                                                                                                                          >  1:1112462/1‑150 (MQ=255)
             ccGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCa                                                                                                                                          >  1:853828/1‑150 (MQ=255)
             ccGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCa                                                                                                                                          >  1:1548187/1‑150 (MQ=255)
                  ccAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATcg                                                                                                                                     >  1:428693/1‑150 (MQ=255)
                   cAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATcgc                                                                                                                                    <  1:102516/150‑1 (MQ=255)
                   cAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATcgc                                                                                                                                    >  1:934036/1‑150 (MQ=255)
                     aTCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCgg                                                                                                                                  <  2:342846/150‑1 (MQ=255)
                     aTCAGAAAACCGTTAGAGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGAATCCGGATTAGATGATTACCGTCACGTTGTGACCCAGTTTCTGCGCAAGCTCGCGGGTAAGGTTGAGATCAATCGACAAATCACGCATATCGCgg                                                                                                                                  >  1:667092/1‑150 (MQ=255)
                       cAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTTTTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGgtt                                                                                                                                <  1:349041/150‑1 (MQ=255)
                         gAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGgttgt                                                                                                                              >  1:530455/1‑150 (MQ=255)
                           aaaCCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTCATTACCGTCACGTTCTGCCCCAGTTTCGGCGCAAGATAGCGGCTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGatgttg                                                                                                                            >  1:1494433/1‑150 (MQ=255)
                              ccGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATg                                                                                                                         <  2:1569679/150‑1 (MQ=255)
                               cGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGc                                                                                                                        >  2:564039/1‑150 (MQ=255)
                                   aGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAAc                                                                                                                    >  2:1524060/1‑150 (MQ=255)
                                    gCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACg                                                                                                                   <  1:1524060/150‑1 (MQ=255)
                                          tGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCttt                                                                                                             <  1:1475237/150‑1 (MQ=255)
                                          gggCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCttt                                                                                                             <  1:564161/149‑1 (MQ=255)
                                           ggCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTg                                                                                                            >  1:15702/1‑150 (MQ=255)
                                           ggCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTg                                                                                                            >  1:740342/1‑150 (MQ=255)
                                                 aCTCGCGCACCTGTGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCTACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCca                                                                                                      >  2:532816/1‑150 (MQ=255)
                                                 aCTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCca                                                                                                      >  1:720622/1‑150 (MQ=255)
                                                          ccTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGCCCTTTGCTCCCAATGCAATgg                                                                                             >  1:521793/1‑150 (MQ=255)
                                                            tGAGCGTAAGTTTGTATGCCGGAGTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGgcg                                                                                           <  1:24524/150‑1 (MQ=255)
                                                               gCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGgcgcgc                                                                                        >  1:1316298/1‑150 (MQ=255)
                                                                cGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCt                                                                                       >  1:1515325/1‑150 (MQ=255)
                                                                     gTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGt                                                                                  <  2:521793/150‑1 (MQ=255)
                                                                        ttGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTcc                                                                               <  2:808888/150‑1 (MQ=255)
                                                                          gATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCctc                                                                             >  2:1495359/1‑150 (MQ=255)
                                                                             gCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTc                                                                          >  2:734040/1‑150 (MQ=255)
                                                                                    tCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATTCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTg                                                                   <  1:1170716/150‑1 (MQ=255)
                                                                                     cGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGa                                                                  >  2:376938/1‑150 (MQ=255)
                                                                                     cGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATAGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGa                                                                  >  1:1573036/1‑150 (MQ=255)
                                                                                     cGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATAGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGa                                                                  >  1:1579442/1‑150 (MQ=255)
                                                                                      gCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGAc                                                                 >  2:557176/1‑150 (MQ=255)
                                                                                      gCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGAc                                                                 <  2:513588/150‑1 (MQ=255)
                                                                                       cTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCATGTTCCTCTTCATTACTGACt                                                                >  2:604835/1‑150 (MQ=255)
                                                                                        tGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACtt                                                               >  2:1550720/1‑150 (MQ=255)
                                                                                        tGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACtt                                                               <  2:428693/150‑1 (MQ=255)
                                                                                         gATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTc                                                              <  1:1619537/150‑1 (MQ=255)
                                                                                         gATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATAGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGAAATGGCGCGCTCCTGTTAGTCTTCTTTACTGACtta                                                              >  2:1687319/1‑149 (MQ=255)
                                                                                           ttACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCgg                                                            <  2:934036/150‑1 (MQ=255)
                                                                                           ttACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCgg                                                            <  2:530455/150‑1 (MQ=255)
                                                                                             aCCGTCACGTTGTGCCCCAGTTTCGGCGGATGCTCGCGGGTTCGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTc                                                          <  1:1449893/150‑1 (MQ=255)
                                                                                              ccGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCa                                                         <  2:1112462/150‑1 (MQ=255)
                                                                                                gTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCACTCGACAAAACACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAAAGCAATGGCGCGCTCCTGTTCCTCTTCAATACTGACTACGGTCAGc                                                       >  2:1066961/1‑150 (MQ=255)
                                                                                                 tCACGATGTGCACCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCa                                                      <  2:1063499/150‑1 (MQ=255)
                                                                                                     gTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCAcccc                                                  <  2:1404792/150‑1 (MQ=255)
                                                                                                      ttGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCAccccc                                                 >  2:805973/1‑150 (MQ=255)
                                                                                                       tgtgCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCa                                                <  2:720622/150‑1 (MQ=255)
                                                                                                          gCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCAtct                                             >  1:1434345/1‑150 (MQ=255)
                                                                                                           ccccAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCAtctc                                            >  1:1414541/1‑150 (MQ=255)
                                                                                                            cccAGTTTCGGCGCAAGCTCGAGGGTACGTGTGAGATCAATCGACAAATCACGCAGATCTCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGTTCAGCACCCCCATCTcc                                           <  1:894581/150‑1 (MQ=255)
                                                                                                                     ggCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGa                                  >  1:342146/1‑150 (MQ=255)
                                                                                                                        gcAATCTCGAGGGTACGGTTGAGATCAATAGACAAATCAAGCAGATCGCGGTTGTTGATGCCAACTACCTTTTCTACCAATGCAATGGCGCGATCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCg                               <  2:667092/150‑1 (MQ=255)
                                                                                                                          aaGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGAc                             <  2:1436835/150‑1 (MQ=255)
                                                                                                                              tCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGAcggcg                         >  1:127577/1‑150 (MQ=255)
                                                                                                                              tCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGAcggcg                         <  1:734040/150‑1 (MQ=255)
                                                                                                                               cgcgGGTACGGTTGAGATCAATAGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGAcggcgg                        <  2:1579442/150‑1 (MQ=255)
                                                                                                                               cgcgGGTACGGTTGAGATCAATAGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGAcggcgg                        <  2:1573036/150‑1 (MQ=255)
                                                                                                                                gcgGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACggcggc                       >  1:688834/1‑150 (MQ=255)
                                                                                                                                 cgGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCa                      <  1:564039/150‑1 (MQ=255)
                                                                                                                                   ggTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAg                    >  2:885683/1‑150 (MQ=255)
                                                                                                                                    gTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGc                   <  2:710834/150‑1 (MQ=255)
                                                                                                                                      aCGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTg                 <  1:557176/150‑1 (MQ=255)
                                                                                                                                      aCGGTGGAGATCAATCGACAAATCACGCAGATGGCGGTTGTTGATGCCAACGACCTTTGCTCCGAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCAACTCCAGACTGTGAGCGACGGCGGCAAGCTg                 <  1:1066961/150‑1 (MQ=255)
                                                                                                                                       gggTTGAGATCAATCGACAAATCAGGCATATCGGGGTTGTGGATGCCAACGACCTTTGCTCCCAATGCAATGGGGGGCGCCTGTTCGTCTTCATTCTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGc               <  2:746281/149‑1 (MQ=255)
                                                                                                                                        ggTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGTTCAGCACCCCCAGCTCCAGACTGTGCGCGACGGCGGCAAGCTGGc               <  2:260619/150‑1 (MQ=255)
                                                                                                                                        ggTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGc               >  1:1482272/1‑150 (MQ=255)
                                                                                                                                         gTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCg              <  2:1316298/150‑1 (MQ=255)
                                                                                                                                            gagaTCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGata           >  2:1136776/1‑150 (MQ=255)
                                                                                                                                            gagaTCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGata           >  1:703594/1‑150 (MQ=255)
                                                                                                                                              gaTCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATAtt         >  2:538546/1‑150 (MQ=255)
                                                                                                                                              gaTCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATAtt         <  2:1414541/150‑1 (MQ=255)
                                                                                                                                                tCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGAGATATTgg       >  2:1248689/1‑150 (MQ=255)
                                                                                                                                                 cAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGgt      <  2:853828/150‑1 (MQ=255)
                                                                                                                                                  aaTCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGgtc     <  2:611930/150‑1 (MQ=255)
                                                                                                                                                  aaTCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGgtc     >  1:749546/1‑150 (MQ=255)
                                                                                                                                                   aTCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGgtcg    <  1:1488078/150‑1 (MQ=255)
                                                                                                                                                   aTCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGgtcg    >  1:113494/1‑150 (MQ=255)
                                                                                                                                                     cGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTACCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGCCTGTGAGCGACGGCGGCAAGCTGGCGATATTGgtcgtc  >  2:9941/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGCCATCAACGCCGAACCAATCAGAAAACCGTTAGCGAAGTGGCTTAACTCGCGCACCTGAGCGTAAGTATTGATGCCGGATTCGCTGATTACCGTCACGTTGTGCCCCAGTTTCGGCGCAAGCTCGCGGGTACGGTTGAGATCAATCGACAAATCACGCAGATCGCGGTTGTTGATGCCAACGACCTTTGCTCCCAATGCAATGGCGCGCTCCTGTTCCTCTTCATTACTGACTTCGGTCAGCACCCCCATCTCCAGACTGTGAGCGACGGCGGCAAGCTGGCGATATTGGTCGTC  >  NC_000913/1319061‑1319359

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