Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,641,826 T→A 7.9% Q460L (CAG→CTG)  PP_3209 ← isoxanthopterin deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,641,8260TA 7.9% 220.3 / 8.1 76Q460L (CAG→CTG) PP_3209isoxanthopterin deaminase
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base A (2/4);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCGGGG  >  NC_002947/3641679‑3641971
                                                                                                                                                    |                                                                                                                                                 
gAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCt                                                                                                                                                   <  5:426680/148‑1 (MQ=255)
gAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTg                                                                                                                                                  <  2:95533/149‑1 (MQ=255)
 aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTgg                                                                                                                                                 <  2:260402/149‑1 (MQ=255)
 aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTgg                                                                                                                                                 >  2:319443/1‑149 (MQ=255)
    ggCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTg                                                                                                                                              <  5:114679/149‑1 (MQ=255)
     gCAGGGGAGCAGGGCCAGGCGTGGCCAGTTCCACACCCAGCACGTCGTCTGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGCACGCCAGGCTgt                                                                                                                                             <  8:381977/149‑2 (MQ=255)
     gCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGc                                                                                                                                             >  5:275432/1‑149 (MQ=255)
     gCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGCACGCCAGGCTgt                                                                                                                                             <  8:376588/149‑2 (MQ=255)
      cAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGcc                                                                                                                                            <  7:298888/149‑1 (MQ=255)
      cAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGcc                                                                                                                                            <  6:79569/149‑1 (MQ=255)
       aGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGccc                                                                                                                                           <  4:387682/149‑1 (MQ=255)
        ggggAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGACTGCCCa                                                                                                                                          <  1:293300/149‑1 (MQ=255)
           gAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGGCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGtt                                                                                                                                       <  5:346220/149‑1 (MQ=255)
             gCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCg                                                                                                                                     <  5:140397/149‑1 (MQ=255)
             gCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCg                                                                                                                                     <  3:286320/149‑1 (MQ=255)
              cAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCgg                                                                                                                                    <  8:13239/149‑1 (MQ=255)
                 gCCCCGGCGGGGCCAGTTCCACACCCAGCTCGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGt                                                                                                                                        <  2:389941/142‑1 (MQ=255)
                 gCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGt                                                                                                                                        >  1:389941/1‑142 (MQ=255)
                  ccccGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAg                                                                                                                                <  1:319443/149‑1 (MQ=255)
                             ccAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                                     <  3:125505/149‑1 (MQ=255)
                                   ccACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATc                                                                                                                <  5:333299/148‑1 (MQ=255)
                                     acacCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCa                                                                                                                      >  6:215844/1‑140 (MQ=255)
                                     acacCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCa                                                                                                                      <  5:215844/140‑1 (MQ=255)
                                        cccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCg                                                                                                          >  3:390555/1‑149 (MQ=255)
                                         ccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGt                                                                                                         >  7:70048/1‑149 (MQ=255)
                                                gtcgtcGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCAc                                                                                                  >  1:92766/1‑149 (MQ=255)
                                                 tcgtcgGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACg                                                                                                 >  8:27173/1‑149 (MQ=255)
                                                    tcgGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCt                                                                                               >  5:369816/1‑148 (MQ=255)
                                                        cTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggc                                                                                          >  3:240360/1‑149 (MQ=255)
                                                         ttCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggcg                                                                                         >  1:374373/1‑149 (MQ=255)
                                                           cAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggcgg                                                                                        <  3:292020/148‑1 (MQ=255)
                                                             aGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGgccg                                                                                     >  5:71699/1‑149 (MQ=255)
                                                              ggCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGgccgc                                                                                    >  3:3274/1‑149 (MQ=255)
                                                                cgcAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCt                                                                                  <  3:254496/149‑1 (MQ=255)
                                                                 gcAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTa                                                                                 >  1:322844/1‑149 (MQ=255)
                                                                  cAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTAc                                                                                <  4:240360/149‑1 (MQ=255)
                                                                   aGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTAcc                                                                               <  2:92766/149‑1 (MQ=255)
                                                                    gggCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCa                                                                              <  6:275432/149‑1 (MQ=255)
                                                                             ccTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGccc                                                                     <  2:322844/149‑1 (MQ=255)
                                                                              cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCg                                                                    >  6:128882/1‑149 (MQ=255)
                                                                              cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCg                                                                    >  6:112605/1‑149 (MQ=255)
                                                                                 cAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGc                                                                   <  4:3274/147‑1 (MQ=255)
                                                                                  agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                                     >  3:19978/1‑96 (MQ=255)
                                                                                  agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                                     <  4:19978/96‑1 (MQ=255)
                                                                                  agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGc                                                                <  6:369816/149‑1 (MQ=255)
                                                                                   gagGCGTCCAGCTCATGCCCGTACCCCCGCGCCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTCCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATACTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCg                                                               >  6:128258/1‑149 (MQ=255)
                                                                                          ccAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAgg                                                                                                                               <  5:6577/78‑1 (MQ=255)
                                                                                          ccAGCTCATGCGCGTACCACCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAgg                                                                                                                               >  6:6577/1‑78 (MQ=255)
                                                                                           cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                                       >  6:332058/1‑149 (MQ=255)
                                                                                           cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                                       >  8:20678/1‑149 (MQ=255)
                                                                                            aGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGc                                                      >  5:183155/1‑149 (MQ=255)
                                                                                               tCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCAc                                                   >  7:247914/1‑149 (MQ=255)
                                                                                               tCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCAc                                                   >  6:22677/1‑149 (MQ=255)
                                                                                                cATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCAGTCCCTGGCTGCCCAGTTCGGCGAGGTCCAGGCCAGGAATCGCATCGTTTTCAACGCTCAGGCGGCCGCCTACCAGCAATCCCCCCAGCGTCGGACGGCCAc                                                   >  1:89395/1‑149 (MQ=255)
                                                                                                 aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAg                                                                                                                     >  6:315343/1‑81 (MQ=255)
                                                                                                 aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGtttt                                                                                                      >  5:90216/1‑96 (MQ=255)
                                                                                                 aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGtttt                                                                                                      <  6:90216/96‑1 (MQ=255)
                                                                                                      cgTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAAc                                            <  6:183155/149‑1 (MQ=255)
                                                                                                           ccccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTg                                       >  8:19408/1‑149 (MQ=255)
                                                                                                           ccccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGACGGACGGCCACCGCCAACCACTg                                       >  4:231024/1‑149 (MQ=255)
                                                                                                           ccccGCGGCCTGATGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCAGCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTg                                       <  5:128258/149‑1 (MQ=255)
                                                                                                             ccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTggg                                     <  6:71699/149‑1 (MQ=255)
                                                                                                                     tGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGcc                             >  3:404604/1‑149 (MQ=255)
                                                                                                                       ttGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGccc                            <  4:390555/148‑1 (MQ=255)
                                                                                                                       ttGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGcc                             <  7:20678/147‑1 (MQ=255)
                                                                                                                            ggCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCa                      <  2:89395/149‑1 (MQ=255)
                                                                                                                              cGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATg                    <  8:247914/149‑1 (MQ=255)
                                                                                                                                 tgggcAGCCAGGCACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAg                 >  5:425426/5‑149 (MQ=255)
                                                                                                                                 tgggcAGCCAGGCACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAg                 >  7:355771/5‑149 (MQ=255)
                                                                                                                                  gTACAGCCTGGCGTGCCTTGCTGCCCAGTTCGGCGATGGCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGa                <  8:70048/149‑1 (MQ=255)
                                                                                                                                  gTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGa                <  7:27173/149‑1 (MQ=255)
                                                                                                                                       gCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGaa           <  5:22677/149‑1 (MQ=255)
                                                                                                                                               acgCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCgg    <  6:425426/147‑1 (MQ=255)
                                                                                                                                               acgCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCgg    <  8:355771/147‑1 (MQ=255)
                                                                                                                                                  ccTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCgggg  <  2:374373/148‑1 (MQ=255)
                                                                                                                                                   cTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCggg   <  4:263411/146‑1 (MQ=255)
                                                                                                                                                   cTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCggg   >  3:263411/1‑146 (MQ=255)
                                                                                                                                                    tGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCTTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACTGCCACCGCCAACCACTGGGGc                                   <  2:243494/113‑1 (MQ=255)
                                                                                                                                                    tGGCTGCCCAGTTCGGCGACGTCCAAGCCAAGAATCGCATCGATTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGc                                   >  1:243494/1‑113 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAGTAGCGGGG  >  NC_002947/3641679‑3641971

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: