Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,321,962 T→C 5.1% T329A (ACA→GCA)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,321,9620TC 5.1% 241.8 / 2.3 79T329A (ACA→GCA) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAG  >  NC_000913/1321818‑1322107
                                                                                                                                                |                                                                                                                                                 
aaaTGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCg                                                                                                                                               >  8:158000/1‑149 (MQ=255)
 aaTGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCgg                                                                                                                                              <  7:285399/149‑1 (MQ=255)
  aTGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCgg                                                                                                                                              >  3:211648/1‑148 (MQ=255)
  aTGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGc                                                                                                                                             <  6:127413/149‑1 (MQ=255)
       tGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCTATCg                                                                                                                                        <  1:319975/149‑1 (MQ=255)
         ccAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCggg                                                                                                                                      <  6:14056/149‑1 (MQ=255)
           aGATCATTACGGGCGAGGTCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTa                                                                                                                                    <  6:111886/149‑1 (MQ=255)
               cATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGc                                                                                                                                             <  8:151013/136‑1 (MQ=255)
               cATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGc                                                                                                                                             >  7:151013/1‑136 (MQ=255)
                 ttACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGAtct                                                                                                                               >  3:29531/1‑148 (MQ=255)
                    cGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAAt                                                                                                                           <  4:209356/149‑1 (MQ=255)
                       gCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTg                                                                                                                        >  5:193062/1‑149 (MQ=255)
                        cGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTgg                                                                                                                       >  7:43436/1‑149 (MQ=255)
                         gagaTCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTggc                                                                                                                      >  6:182/1‑149 (MQ=255)
                         gagaTCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGGTCCGGCAATCGGGTAGATCTCAATCTggc                                                                                                                      <  4:211648/149‑1 (MQ=255)
                            aTCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTggcggc                                                                                                                   >  6:229902/1‑149 (MQ=255)
                             tCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTggcggc                                                                                                                   <  4:318396/148‑1 (MQ=255)
                               aaCCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtgg                                                                                                                <  6:193062/149‑1 (MQ=255)
                                aCCAGCATCCGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtggt                                                                                                               <  6:199320/149‑1 (MQ=255)
                                aCCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtggt                                                                                                               <  2:1955/149‑1 (MQ=255)
                                 ccAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGCACCATCGGCGCGACGACCGCGTGGGCGTGCTCCGGCAATCGGGGAGATCTCAATCTGGCGGCtggtg                                                                                                              >  2:49928/1‑149 (MQ=255)
                                 ccAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtggtg                                                                                                              <  5:182/149‑1 (MQ=255)
                                 ccAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGCTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtggtg                                                                                                              >  2:44209/1‑149 (MQ=255)
                                       taataTGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTAAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGGGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATGTCAATCTGGCGGCTGGTGGcatca                                                                                                        >  8:356520/5‑149 (MQ=255)
                                            tGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACtt                                                                                                   <  7:182148/149‑1 (MQ=255)
                                             gTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACtt                                                                                                   >  8:347299/1‑148 (MQ=255)
                                             gTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCATGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACtt                                                                                                   >  1:258240/1‑148 (MQ=255)
                                              ttCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTga                                                                                                 <  4:60301/149‑1 (MQ=255)
                                                 aGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCATGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTgagc                                                                                               <  2:258240/148‑1 (MQ=255)
                                                 aGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGAGCTCTTTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTgagc                                                                                               <  2:268524/148‑1 (MQ=255)
                                                   aCAGCTCTTTATGATCGGTACGCATTTCCAGTGCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTgagcgag                                                                                            <  2:357321/149‑1 (MQ=255)
                                                    cAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGCTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTgagcgagc                                                                                           <  1:44209/149‑1 (MQ=255)
                                                    cAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGGGCTCCGGCACTCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTgagcgagc                                                                                           <  1:49928/149‑1 (MQ=255)
                                                     aGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCtgg                                                                                                                >  8:358185/1‑127 (MQ=255)
                                                     aGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCAAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCt                                                                                          >  4:368783/1‑149 (MQ=255)
                                                       ctctTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTACATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCttt                                                                                        >  3:34771/1‑149 (MQ=255)
                                                          tttATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcatc                                                                                                         <  1:399430/129‑1 (MQ=255)
                                                          tttATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcatc                                                                                                         >  2:399430/1‑129 (MQ=255)
                                                          tttATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCg                                                                                     >  5:233804/1‑149 (MQ=255)
                                                            tATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCgg                                                                                    <  5:229902/148‑1 (MQ=255)
                                                             aTGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGc                                                                                   >  3:109500/1‑148 (MQ=255)
                                                               gATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGCCGACCGCGTGGGCGTGTTCCGGAAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCTGCGAc                                                                                >  4:261740/1‑149 (MQ=255)
                                                                    gTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACgcgc                                                                            >  2:191196/1‑148 (MQ=255)
                                                                        gCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGATCCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCg                                                                                     <  8:188808/135‑1 (MQ=255)
                                                                        gCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCg                                                                                     >  7:188808/1‑135 (MQ=255)
                                                                          aTTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAg                                                                     >  3:61838/1‑149 (MQ=255)
                                                                           tttCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAgg                                                                    <  4:109500/149‑1 (MQ=255)
                                                                              ccAGTTCACTACGGCTGTCGAGATCTCTGGCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGt                                                                  <  6:215831/148‑1 (MQ=255)
                                                                              ccAGTTCAATACGGCGGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGt                                                                  <  2:207031/148‑1 (MQ=255)
                                                                                aGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGa                                                                >  3:226823/1‑148 (MQ=255)
                                                                                     aaTACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGAAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCAt                                                          >  7:193008/1‑149 (MQ=255)
                                                                                         cGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATc                                                      <  8:43436/149‑1 (MQ=255)
                                                                                            cTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTg                                                   <  7:158000/149‑1 (MQ=255)
                                                                                              gTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGAcgc                                                                              >  6:134778/1‑120 (MQ=255)
                                                                                              gTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGAcgc                                                                              <  5:134778/120‑1 (MQ=255)
                                                                                               tCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCa                                                 <  7:12284/148‑1 (MQ=255)
                                                                                                      ctctctCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATaaaaaa                                          >  7:144190/1‑148 (MQ=255)
                                                                                                      ctctGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGAc                                                                                <  4:28133/110‑1 (MQ=255)
                                                                                                         tGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCAt                                                          >  6:140440/1‑129 (MQ=255)
                                                                                                         tGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCAt                                                          <  5:140440/129‑1 (MQ=255)
                                                                                                         tGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAAc                                         <  6:233804/146‑1 (MQ=255)
                                                                                                           tCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGt                                     <  3:78023/148‑1 (MQ=255)
                                                                                                           tCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTa                                    >  5:77617/1‑149 (MQ=255)
                                                                                                            ccAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGCCGCGCCAAATAGGGTGCAATCATTATCCTGCATAAAAAACATGTAc                                   >  4:365241/1‑149 (MQ=255)
                                                                                                              aGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACgg                                 >  7:103563/1‑149 (MQ=255)
                                                                                                               gTGAACCATCGGCGCTACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGATGGTGGCATCATACTTGAGCGATCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACggg                                <  5:190766/149‑1 (MQ=255)
                                                                                                                  aaCCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCt                              <  4:61838/148‑1 (MQ=255)
                                                                                                                     cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGa                          >  4:239602/1‑149 (MQ=255)
                                                                                                                     cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGAAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGa                          <  8:193008/149‑1 (MQ=255)
                                                                                                                              gacgacCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCttt                 <  7:347299/149‑1 (MQ=255)
                                                                                                                              gacgacCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCttt                 >  8:167776/1‑149 (MQ=255)
                                                                                                                                cgacCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCGGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCttttt               <  4:29531/149‑1 (MQ=255)
                                                                                                                                cgacCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCGGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCttttt               <  2:263242/149‑1 (MQ=255)
                                                                                                                                 gacCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCaaa                                                                        >  6:30590/1‑91 (MQ=255)
                                                                                                                                  acCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCa             <  8:144190/149‑1 (MQ=255)
                                                                                                                                          ggCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAAc                                         >  3:137713/1‑113 (MQ=255)
                                                                                                                                          ggCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAAc                                         <  4:137713/113‑1 (MQ=255)
                                                                                                                                           gCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCAAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTaa     >  5:111390/1‑148 (MQ=255)
                                                                                                                                           gCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTATCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACTGGCTGGGATTa                       >  6:296484/1‑130 (MQ=255)
                                                                                                                                             gtgtTCCGGCAATCGGGGAGAGCTCAATCGGGCGGCTGGGGGCATCAGCCTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAg  <  2:155634/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATACTTGAGCGAGCTTTCCGGCGACGCGCCAAATAGGGTGAAATCATTATCCTGCATAAAAAACATGTACGGGCTGGGATTACTCTTTTTCAGCACGTAATAG  >  NC_000913/1321818‑1322107

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