Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,666,986 A→G 5.2% T287T (ACT→ACC hcaT ← putative 3‑phenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,666,9860AG 5.2% 312.0 / 4.3 95T287T (ACT→ACChcaTputative 3‑phenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (50/40);  new base G (2/3);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCC  >  NC_000913/2666838‑2667129
                                                                                                                                                    |                                                                                                                                               
cgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTa                                                                                                                                                 <  3:752836/149‑1 (MQ=255)
cgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTa                                                                                                                                                 <  1:679567/149‑1 (MQ=255)
  cccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGtt                                                                                                                                              <  2:869946/149‑1 (MQ=255)
   ccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGtt                                                                                                                                              <  1:1048850/149‑1 (MQ=255)
   ccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGtt                                                                                                                                              >  4:847209/1‑149 (MQ=255)
        gcgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTcc                                                                                                                                         >  1:117393/1‑149 (MQ=255)
          gACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTgatagtg                                                                                                                                       <  1:474220/149‑7 (MQ=255)
              gCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAAt                                                                                                                                   <  3:974323/149‑1 (MQ=255)
                agagTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGc                                                                                                                                 >  3:856121/1‑149 (MQ=255)
                agagTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGc                                                                                                                                 <  4:621320/149‑1 (MQ=255)
                agagTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGc                                                                                                                                 >  2:679494/1‑149 (MQ=255)
                  agTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGccc                                                                                                                               >  2:1067149/1‑149 (MQ=255)
                    tAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCa                                                                                                                             >  1:886275/1‑149 (MQ=255)
                      aaCCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAgc                                                                                                                           <  2:917031/149‑1 (MQ=255)
                          gCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCAc                                                                                                                       <  2:786348/149‑1 (MQ=255)
                           ccTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACt                                                                                                                      >  1:588383/1‑149 (MQ=255)
                              gTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcg                                                                                                                   >  4:383337/1‑149 (MQ=255)
                               tAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcgc                                                                                                                  <  2:208101/149‑1 (MQ=255)
                                aaaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcgcc                                                                                                                 >  1:41219/1‑149 (MQ=255)
                                aaaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcgcc                                                                                                                 <  1:864214/149‑1 (MQ=255)
                                aaaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcgcc                                                                                                                 <  4:688461/149‑1 (MQ=255)
                                 aaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTc                                                                                                                                          <  3:525935/123‑1 (MQ=255)
                                 aaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTc                                                                                                                                          >  4:525935/1‑123 (MQ=255)
                                        gACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCg                                                                                                         >  2:638758/1‑149 (MQ=255)
                                         aCTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGc                                                                                                        >  3:207498/1‑149 (MQ=255)
                                          cTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGcc                                                                                                       <  3:383337/149‑1 (MQ=255)
                                          cTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGcc                                                                                                       <  1:166047/149‑1 (MQ=255)
                                          cTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGcc                                                                                                       >  2:415923/1‑149 (MQ=255)
                                           ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCg                                                                                                      >  3:82892/1‑149 (MQ=255)
                                           ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCg                                                                                                      <  2:397135/149‑1 (MQ=255)
                                           ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCg                                                                                                      >  2:530242/1‑149 (MQ=255)
                                           ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCg                                                                                                      >  4:1058462/1‑149 (MQ=255)
                                                  cccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCaa                                                                                               >  3:214771/1‑149 (MQ=255)
                                                   ccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAAc                                                                                              <  2:773113/149‑1 (MQ=255)
                                                    cTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACa                                                                                             <  2:117393/149‑1 (MQ=255)
                                                    cTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACa                                                                                             <  2:462988/149‑1 (MQ=255)
                                                     tGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAg                                                                                            >  2:138145/1‑149 (MQ=255)
                                                         gAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAAGAGCCGCAATGGAGAATTTGCACCACTATCAACCACGGCACCGCCGTAGTTGCTCCCATAATGCGCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCata                                                                                        >  2:27415/1‑149 (MQ=255)
                                                               caATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGt                                                                                  >  1:668884/1‑149 (MQ=255)
                                                                 atatAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGc                                                                                >  4:943182/1‑149 (MQ=255)
                                                                     aaCGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTga                                                                                                                                            <  1:247798/85‑2 (MQ=38)
                                                                     aaCGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTga                                                                                                                                            >  4:503589/1‑84 (MQ=38)
                                                                     aaCGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTga                                                                                                                                            >  2:247798/1‑84 (MQ=38)
                                                                     aaCGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTga                                                                                                                                            <  3:503589/85‑2 (MQ=38)
                                                                        gCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                                                         <  1:638758/149‑1 (MQ=255)
                                                                        gCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                                                         >  3:618894/1‑149 (MQ=255)
                                                                         cATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAc                                                                        <  1:311760/149‑1 (MQ=255)
                                                                            ggcggcCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                                                     <  4:571337/149‑1 (MQ=255)
                                                                            ggcggcCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                                                     <  1:415923/149‑1 (MQ=255)
                                                                             gcggcCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcg                                                                    <  3:943182/149‑1 (MQ=255)
                                                                                  cAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaaa                                                               <  2:1007491/149‑1 (MQ=255)
                                                                                   aGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAg                                                              >  4:114134/1‑149 (MQ=255)
                                                                                   aGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAg                                                              >  4:996122/1‑149 (MQ=255)
                                                                                    ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGt                                                             <  4:445678/149‑1 (MQ=255)
                                                                                    ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGt                                                             >  3:339359/1‑149 (MQ=255)
                                                                                       ggCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTtta                                                          <  1:803952/149‑1 (MQ=255)
                                                                                          aGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                                                       <  3:996122/149‑1 (MQ=255)
                                                                                          aGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                                                       <  3:114134/149‑1 (MQ=255)
                                                                                          aGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                                                       >  3:456390/1‑149 (MQ=255)
                                                                                          aGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAACATATCGCGTGCACTACAACGGCGGAAAAGTttatta                                                       <  1:641148/149‑1 (MQ=255)
                                                                                                    ggTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaa                                             >  1:535932/1‑149 (MQ=255)
                                                                                                     gTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaaa                                            <  3:1058462/149‑1 (MQ=255)
                                                                                                     gTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaaa                                            <  3:341698/149‑1 (MQ=255)
                                                                                                          gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                       <  4:613250/149‑1 (MQ=255)
                                                                                                          gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                       <  2:846425/149‑1 (MQ=255)
                                                                                                           cAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt                                      <  1:1067149/149‑1 (MQ=255)
                                                                                                             aTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCaa                                                                                               <  1:979061/90‑1 (MQ=255)
                                                                                                             aTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCaa                                                                                               >  2:979061/1‑90 (MQ=255)
                                                                                                              tGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTAATTACTCAGCGCAAAGATAATGAc                                   >  2:719420/1‑149 (MQ=255)
                                                                                                                  gAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTcc                               >  3:310543/1‑149 (MQ=255)
                                                                                                                       tGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGccac                          <  2:535932/149‑1 (MQ=255)
                                                                                                                         caccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGccacca                        >  2:94813/1‑149 (MQ=255)
                                                                                                                          accacTATGAACTACGGCAACGCCGTAGTTGCTCCCCTAATGCCCCAGCGCCCTAGGCCGCACATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGATTATTACTCAGCGCAAAAATAACGACTTCCGCcaccac                       >  4:949105/1‑149 (MQ=255)
                                                                                                                          accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCcaccac                       >  3:58734/1‑149 (MQ=255)
                                                                                                                          accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCcaccac                       >  1:71956/1‑149 (MQ=255)
                                                                                                                          accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCcaccac                       >  2:827118/1‑149 (MQ=255)
                                                                                                                            cacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGc                     >  4:243025/1‑149 (MQ=255)
                                                                                                                               tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCa                  >  2:987581/1‑149 (MQ=255)
                                                                                                                               tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCa                  >  4:734847/1‑149 (MQ=255)
                                                                                                                                aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAg                 >  3:213110/1‑149 (MQ=255)
                                                                                                                                   aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa              <  4:213110/149‑1 (MQ=255)
                                                                                                                                   aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa              >  3:209872/1‑149 (MQ=255)
                                                                                                                                    aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGAc             >  3:819045/1‑149 (MQ=255)
                                                                                                                                    aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGAc             >  4:971207/1‑149 (MQ=255)
                                                                                                                                    aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGAc             <  3:107918/149‑1 (MQ=255)
                                                                                                                                    aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGAc             >  1:516775/1‑149 (MQ=255)
                                                                                                                                    aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGAc             >  3:911140/1‑149 (MQ=255)
                                                                                                                                      cACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACca           >  1:904245/1‑149 (MQ=255)
                                                                                                                                           cAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAAt      <  4:791755/149‑1 (MQ=255)
                                                                                                                                           cAACGCCGTAGTTGCTCCCATATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaaa                                                               >  1:164231/1‑91 (MQ=255)
                                                                                                                                           cAACGCCGTAGTTGCTCCCATATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaaa                                                               <  2:164231/91‑1 (MQ=255)
                                                                                                                                            gttgCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATa     >  4:390410/4‑149 (MQ=255)
                                                                                                                                            aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATa     <  4:819045/149‑1 (MQ=255)
                                                                                                                                            aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATa     >  4:488321/1‑149 (MQ=255)
                                                                                                                                            aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATa     >  4:275423/1‑149 (MQ=255)
                                                                                                                                               gCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAccc  <  4:82892/149‑1 (MQ=255)
                                                                                                                                               gCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAccc  >  2:671413/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCC  >  NC_000913/2666838‑2667129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: