Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,321,616 C→T 5.6% R444H (CGT→CAT)  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,6160CT 5.6% 228.4 / 2.3 71R444H (CGT→CAT) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACA  >  NC_000913/1321468‑1321762
                                                                                                                                                    |                                                                                                                                                  
cTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGcgac                                                                                                                                                    <  3:388807/149‑1 (MQ=255)
   cGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCgacgac                                                                                                                                                 >  6:23856/1‑149 (MQ=255)
     gAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCt                                                                                                                                               >  8:297672/1‑149 (MQ=255)
        cAGAATCAAGTACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCATCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGAGGTGAAATAACCTACCGCGCCGCCGTATCTGCCGCGGCGACGACCTTcc                                                                                                                                            <  4:245374/149‑1 (MQ=255)
        cAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTcc                                                                                                                                            <  4:338879/149‑1 (MQ=255)
           aaTCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCCGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGcc                                                                                                                                         <  8:327457/149‑1 (MQ=255)
             tCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGTCGACGACCTTCCGCCTc                                                                                                                                       <  2:273396/149‑1 (MQ=255)
               aaGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACAAGCGCCGAGCTGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGATGACCTTCCGCCTCgg                                                                                                                                     >  8:359469/1‑149 (MQ=255)
                aGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGc                                                                                                                                    <  6:185088/149‑1 (MQ=255)
                  gACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACTAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGATGACCTTCCGCCTCGGCaa                                                                                                                                  >  4:150526/1‑149 (MQ=255)
                    cTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAAtt                                                                                                                                <  8:35877/149‑1 (MQ=255)
                     tACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAAtt                                                                                                                                >  8:87350/1‑148 (MQ=255)
                     tACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAAtt                                                                                                                                >  4:98211/1‑148 (MQ=255)
                     tACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGACATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGATGCCGCGGCGACGACCTTCCGCCTCGGCAAtt                                                                                                                                >  3:171560/1‑148 (MQ=255)
                          cAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACACTGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGc                                                                                                                          <  1:189617/149‑1 (MQ=255)
                          cAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGc                                                                                                                          <  7:87350/149‑1 (MQ=255)
                              aCCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGcgacga                                                                                                                                                  <  8:58446/121‑1 (MQ=255)
                                ccGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAgcg                                                                                                                    >  5:374721/1‑149 (MQ=255)
                                      gCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTa                                                                                                                <  5:23856/147‑1 (MQ=255)
                                         cGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTc                                                                                                           <  7:113771/149‑1 (MQ=255)
                                          ggTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCg                                                                                                          <  5:260607/149‑1 (MQ=255)
                                          ggTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCg                                                                                                          <  4:254614/149‑1 (MQ=255)
                                           gTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCgg                                                                                                         <  4:181710/149‑1 (MQ=255)
                                                gATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCAc                                                                                                    <  3:96581/149‑1 (MQ=255)
                                                gATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCAGTAGCAGCCTCGGCGACGACCTTCCGCCTCGGCAATTAACAGCATAGCTCGTACTTTCGGAGCAc                                                                                                    >  1:192795/1‑149 (MQ=255)
                                                    ccGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGc                                                                                                 <  3:98211/148‑1 (MQ=255)
                                                       ttttCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTaa                                                                                             >  8:406496/1‑149 (MQ=255)
                                                       ttttCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTaa                                                                                             <  2:333868/149‑1 (MQ=255)
                                                         ttCCACCAGGGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGACTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTatc                                                                                            <  3:83591/148‑3 (MQ=255)
                                                            caccaGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTaa                                                                                             >  5:349682/1‑144 (MQ=255)
                                                            caccaGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTaa                                                                                             <  6:349682/144‑1 (MQ=255)
                                                            caccaGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTcc                                                                                        >  8:42346/1‑149 (MQ=255)
                                                            caccaGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCACGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTcc                                                                                        >  1:418185/1‑149 (MQ=255)
                                                              ccaGCGCCGAGCGGATCACAATGCAGGTGGCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTcccc                                                                                      <  5:360974/149‑1 (MQ=255)
                                                                   gcCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATAtt                                                                                 <  3:103293/149‑1 (MQ=255)
                                                                   gcCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATAtt                                                                                 >  3:344352/1‑149 (MQ=255)
                                                                       aGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGATGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATa                                                                             <  7:359469/149‑1 (MQ=255)
                                                                       aGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATa                                                                             <  7:297672/149‑1 (MQ=255)
                                                                          ggATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAg                                                                          <  8:118490/149‑1 (MQ=255)
                                                                                 aaTGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGAACTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGa                                                                    >  8:127986/1‑148 (MQ=255)
                                                                                  aTGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGAt                                                                   >  8:318504/1‑148 (MQ=255)
                                                                                   tGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGAt                                                                   >  3:279901/1‑147 (MQ=255)
                                                                                   tGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATaa                                                                 <  7:42346/149‑1 (MQ=255)
                                                                                        gtgtCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTg                                                            >  3:345608/1‑149 (MQ=255)
                                                                                         tgtCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTg                                                            >  8:261227/1‑148 (MQ=255)
                                                                                           tCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAg                                                         <  7:406496/149‑1 (MQ=255)
                                                                                             gagaTCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAggg                                                       >  7:400653/1‑149 (MQ=255)
                                                                                              agaTCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGATGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGc                                                      <  3:150526/149‑1 (MQ=255)
                                                                                                aTCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGctt                                                    >  6:367848/1‑147 (MQ=255)
                                                                                                 tCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTc                                                   <  5:226598/149‑1 (MQ=255)
                                                                                                     cATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGa                                               <  7:318504/149‑1 (MQ=255)
                                                                                                       tGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCACTTAACTGCATAGCGCCTACTTTCGGCGCACCGCGTACCGTCCCCATATTCATACCGGCGCGATAAGCGTGCAGGGCGTc                                                   >  1:254242/1‑143 (MQ=255)
                                                                                                       tGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCa                                                          >  7:352876/1‑136 (MQ=255)
                                                                                                       tGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCa                                                          >  8:418941/1‑136 (MQ=255)
                                                                                                       tGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCa                                                          <  8:352876/136‑1 (MQ=255)
                                                                                                       tGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCa                                                          <  7:418941/136‑1 (MQ=255)
                                                                                                           cgGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGa                                         <  6:374721/149‑1 (MQ=255)
                                                                                                                aaaTAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAg                                    >  6:40466/1‑149 (MQ=255)
                                                                                                                 aaTAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGt                                   <  7:261227/149‑1 (MQ=255)
                                                                                                                   tAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTc                                 <  5:373469/149‑1 (MQ=255)
                                                                                                                    aaCCTAACGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTAGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACAGCTTAACGTCCCAATATTAATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCg                                >  2:13991/1‑149 (MQ=255)
                                                                                                                     aCCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCaa                                                 >  3:88736/1‑131 (MQ=255)
                                                                                                                     aCCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCaa                                                 <  4:88736/131‑1 (MQ=255)
                                                                                                                     aCCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGc                               >  1:172719/1‑149 (MQ=255)
                                                                                                                      ccTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGcc                              >  5:161259/1‑149 (MQ=255)
                                                                                                                        tACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGa                            >  4:384549/1‑149 (MQ=255)
                                                                                                                          ccGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGAc                           <  2:161792/148‑1 (MQ=255)
                                                                                                                            gcgcCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTAc                        >  2:364253/1‑149 (MQ=255)
                                                                                                                              gccgccGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGa                            <  2:388959/143‑1 (MQ=255)
                                                                                                                              gccgccGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGa                            >  1:388959/1‑143 (MQ=255)
                                                                                                                                      aGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGAc               >  7:320610/1‑148 (MQ=255)
                                                                                                                                           ccGCGGCGACGTCCTTCCGCCTCGGCATTTAACTGCATAGCGCGTACTTTCTGCGCACGGCTAAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGt          <  2:192795/148‑1 (MQ=255)
                                                                                                                                             gcggcgACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCa       >  4:298673/1‑149 (MQ=255)
                                                                                                                                              cggcgACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATTAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCAt      >  1:94283/1‑149 (MQ=255)
                                                                                                                                                  gacgacCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATcaca  >  1:361599/1‑149 (MQ=255)
                                                                                                                                                    cgacCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAgg                                                        <  2:136644/93‑1 (MQ=255)
                                                                                                                                                    cgacCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAgg                                                        >  1:136644/1‑93 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACA  >  NC_000913/1321468‑1321762

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