Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,509,356 A→T 13.6% Y152* (TAT→TAA)  yhfW ← phosphopentomutase‑related metalloenzyme superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,509,3560AT13.6% 144.9 / 2.6 66Y152* (TAT→TAA) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (13/43);  new base T (8/2);  total (21/45)
Fisher's exact test for biased strand distribution p-value = 1.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGACGCCATAG‑CCCATATGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCTTCATGTTGCAGCTCTGCC  >  NC_000913/3509120‑3509590
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
tCGACGCCATAG‑CCCATATGCACGACCTGGAAACCTTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGCCCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTtata                                                                                                                                                                                                                                              <  2:923893/239‑1 (MQ=255)
   aCGCCATAG‑CCCATATGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGCCAAAAGAGAT‑ACCAGAGA‑GA‑TTGGC‑GGCAAA‑GTtttaaa                                                                                                                                                                                                                                            >  1:1059723/1‑236 (MQ=255)
    cggcATCG‑CCCATAGGCCCGACATGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCAGTTATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCAACAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GAcc                                                                                                                                                                                                                                         <  2:785305/236‑1 (MQ=255)
    cGCCATAGCCCCATATGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GAc                                                                                                                                                                                                                                          <  1:735186/239‑1 (MQ=255)
    cGCCATAG‑CCCATATGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGGGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GAc                                                                                                                                                                                                                                          <  1:1112932/238‑1 (MQ=255)
    cGCCATAG‑CCCATACGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GAcc                                                                                                                                                                                                                                         <  1:292195/239‑1 (MQ=255)
     gCCATAG‑CCCATATGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATTCCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GAcc                                                                                                                                                                                                                                         <  1:1326586/238‑1 (MQ=255)
      ccATAG‑CCCATATGCACGGCCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCt                                                                                                                                                                                                                                        <  1:488714/238‑1 (MQ=255)
               cATATGCCCGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGTATTGATACCAATAAAGCGCTCTTTTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAAt                                                                                                                                                                                                                               <  1:624526/239‑1 (MQ=255)
                  atGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTCAT‑GTTTTA‑AACCTGGCCCAAATc‑c                                                                                                                                                                                                                            >  1:95738/1‑238 (MQ=255)
                  atGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGCCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAA‑GTTATA‑AACCTGGCCTtaaaac                                                                                                                                                                                                                             >  1:747958/1‑235 (MQ=255)
                     gcaggccTGGAAACCGTTGTCATTAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGc                                                                                                                                                                                                                          <  1:997564/232‑1 (MQ=255)
                       cGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGCCCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTGGTCAAAAGAGAT‑AACAAGAA‑GA‑TTGGC‑GGGAATAGTTAAG‑AACCTGG‑CTACAAC‑CGCCTC‑g                                                                                                                                                                                                                     >  1:1340792/1‑239 (MQ=255)
                       cGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAACGGGAT‑CACAGAGA‑GA‑TTGGCGGGTAAT‑TTTTTA‑GCCTTGGctcaatcccccctc                                                                                                                                                                                                                        >  1:719173/1‑226 (MQ=255)
                       cGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGCCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGCC‑GGGAAA‑GTTTTA‑AACCTGGCCCaaaaaccccccc                                                                                                                                                                                                                       >  2:26989/1‑230 (MQ=255)
                          ccTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATATAACCTGCCCTAAACC‑CCCCCC‑gag                                                                                                                                                                                                                   >  1:27621/1‑239 (MQ=255)
                             ggAAACCGTTGTCATAACCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGCCGTCAAAAGAGAT‑CACAGAGG‑GA‑TTGCC‑GGCAAT‑GTTTTT‑GACCTGGCCCAAATC‑CACCTC‑CGGAttg                                                                                                                                                                                                               >  1:938431/1‑238 (MQ=255)
                                    gTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCCCCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CCCCGAGA‑GA‑TTGGC‑GGGAAT‑TTTTTT‑AACCCGGCCTaaaaaccccccggagttatctttc                                                                                                                                                                                                           >  2:997564/1‑217 (MQ=255)
                                     ttGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTTTT‑GACCCGGCCTAAATC‑CCCCTC‑CAGATTATAGccca                                                                                                                                                                                                        >  1:962774/1‑236 (MQ=255)
                                        tCATAAGCTCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCCCGATCCGCCCAATTTTGATTGCGTCGCCAAAACAGAC‑AACAGA‑A‑GACTTGGCGGGAAAT‑GTTATA‑GACCAGCCTTAAATA‑CCCCCG‑AAGTTAATCGCCaacac                                                                                                                                                                                                     >  2:714125/1‑236 (MQ=255)
                                            aaGCCCCAAAACGCGGCGCTTTGATACCAAAAACGCCCCCTTTTTTGTTTTCGGGCGATTCGAGAATGCGTTGACTGTCTGTTAACAGGCCCCCAAA‑GGTAATTGCCCGCCCGACCTGTCCCTGCTCACGCACGATACGACCATTTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAAcc                                                                                                                                                                                                <  2:1059723/239‑1 (MQ=255)
                                             aGCGCCAGAACGCGGCGCATTGATCCCAATAAAGCCCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCg                                                                                                                                                                                               <  1:1371090/239‑1 (MQ=255)
                                             aGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GC‑TTGGC‑GGTAAT‑GGTATA‑AACCCGGCCTAAATC‑CGCCTC‑GGGGTTTTTGCCAATACCAAcc                                                                                                                                                                                                >  1:62935/1‑238 (MQ=255)
                                                              attttttccaaaaaaaGCGTCTTTCTTTGTTTTCGGGGCCATCGAGAATGGGTTGCCTGTGGGTTAACAGCCCACAAAA‑TGTAATGCCCCGCCCGCCCTGTCCCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CCAAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GCCCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGa                                                                                                                                                                              <  2:1120683/226‑1 (MQ=255)
                                                               aTTGATCCCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGa                                                                                                                                                                              <  1:1066813/238‑1 (MQ=255)
                                                                    tACCAATAAAGCGCCCTTCTTTGCTTTCGGCGGCAACGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGCCCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGc                                                                                                                                                                        <  1:288869/239‑1 (MQ=255)
                                                                         aTAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGCGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑TTTATT‑AGCCTGGTCTAAAAC‑ACCCTC‑CGGATTATTGCCCCAACCAACACCCTGATTGACCCACCGAAATTGcccagc                                                                                                                                                                   >  1:961383/1‑234 (MQ=255)
                                                                         aTAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAAACGCCCGATTTACCCCCAGaaaattcacata                                                                                                                                                                   >  1:519867/1‑230 (MQ=255)
                                                                         aTAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAA‑GTTATA‑AACCCGGCCTAAATC‑CCCCTC‑GAGATTATCGCCAATACGAACCGCCTGATTTGCCCCCAGCAAATACAGatc                                                                                                                                                                   >  1:399270/1‑239 (MQ=255)
                                                                                gcCCTTTTTTGCTTTCTGGGGCATCGGGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCa                                                                                                                                                            <  1:380101/239‑1 (MQ=255)
                                                                                 cccTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTTTA‑GACCTGGCCTAAAAT‑CGCCCC‑GAGATTATCGGCCATCGCAACCCCCCGATTTACCCCCAGAAAATGCAGATCATccccc                                                                                                                                                            >  1:478276/1‑237 (MQ=255)
                                                                                 cccTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑AAACTGGCCTAAATC‑CGCCTC‑GAGATTATCTCCAATTCCAACCGCCCGGTTGAACCACAACAAATTCAGAACATCAcccc                                                                                                                                                           >  1:643109/1‑237 (MQ=255)
                                                                                         tGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCGGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTcc                                                                                                                                                  <  1:608520/239‑1 (MQ=255)
                                                                                                        cGAGAATGCGTTGCCTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGCCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑TACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCg                                                                                                                                   <  1:276811/239‑1 (MQ=255)
                                                                                                         gggaattttttgctggTCGTTTAACAGGCCACAAAA‑TGTAAGGCCCCGACAGACCTGTCCCTGCCACCGCACGAATCGACCAATTTTGATTGGGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAA‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCgg                                                                                                                                  <  2:399270/224‑1 (MQ=255)
                                                                                                         gggaACGCGTTGACTGTCTGTTAACCGGCCACCAAA‑GTAAATGCCCCAACCGCCCTTGACCGGCCCACCCACGATAGTGCCAATTTTGATTGCTGCGT‑AAAAAAGATCCACAGAAA‑GA‑TTGGC‑GGTAAG‑GTTAAA‑GCCTGGGCCAAAACC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCgg                                                                                                                                  <  2:185914/237‑1 (MQ=255)
                                                                                                         gaaaaTGGGTTGACTGGCGGTTAACAGGCCCACAAA‑TGTAAGGCCCCGACCGCCCTGTCCCTGCTCACGCACGATACGCCCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCgg                                                                                                                                  <  1:26989/236‑1 (MQ=255)
                                                                                                                 ttttgttggtggTTAAACGGCCCCCAAA‑TTAAATGACCCGACCGACCTGTCCCTGCTCCCGCACGATACGACCAATTTTGATTGGGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAATCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAAt                                                                                                                          <  2:519867/229‑1 (MQ=255)
                                                                                                                 tttGCCTTTCGTTTAACCGGCCACAAAA‑TGTAATGCCCCGACCGACCGGTACCTGCCCCCGGACGGTACCACCAATTTTATCTGGTTCGTAAAAAAGAAT‑CACAGGGA‑GA‑TTGGC‑GGTAAT‑GTTAAA‑GCCCTGGCCTAAACC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAAt                                                                                                                          <  2:18979/238‑1 (MQ=255)
                                                                                                                      cTGCCGGTTACAAGGCCCACAAA‑TTTAAGGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTg                                                                                                                     <  1:1254648/239‑1 (MQ=255)
                                                                                                                             ttacagtgCCACAAAT‑TGTAAGGCCCCGACCGCCCTTTCCCTCCCCACCCCAGGTAAGACCATTTTTGTTGGGCTCCGC‑AAAGAGAA‑CACAGGGAGGA‑TGGGC‑GGTAAT‑GTTAAA‑GACCTGGCCTAAATC‑CGCCTC‑GGGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAaca                                                                                                              <  2:95738/232‑1 (MQ=255)
                                                                                                                              tAACAGGCCAACAAA‑TGTAATGCCCCCACCGGCCTGTACCTGCTCACGCACGATACGACCAATTTGGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAaca                                                                                                              <  1:802237/238‑1 (MQ=255)
                                                                                                                               taaaagCCCCCAAATTGTATTGCCCCGACCGACCGGTACCTGCTCACGCAGGATACGCCCAATTTGGTTTGCGTCGTAAAAAGAGAC‑CCCAGAGG‑GA‑TTGGC‑GGTAAT‑GTTATA‑GCCCGGGCCTAAACC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAaca                                                                                                              <  2:340445/233‑1 (MQ=255)
                                                                                                                               acccggCCCCAAAA‑TTAAATGCCCCAGCCGACCTTTACCTGCTCACGCAGGAAACGACCAATTTGGTTGGCGTCGTCAAAAGAGAC‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GCCCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACg                                                                                                            <  2:914065/235‑1 (MQ=255)
                                                                                                                                                tAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGat                                                                                                                                                                    <  1:15791/168‑1 (MQ=255)
                                                                                                                                                       ccgaccgaCCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTCAT‑GTTTTA‑GAACTGGCCTAAAAC‑CGCCTC‑GCGATTATCGCCAAACGCACACCCCCGCTTGACCCCCAGCAAATTCAGACAATCCCCACCGCGGCCCCCCCTCCCACCCCCCGACAAC‑AAAGCCTGCTCAAAACGCTCCACCACACCCCCaaaaacaa                                                                                     >  2:102648/1‑235 (MQ=255)
                                                                                                                                                          accgaccTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAAACCCCCGATTGACCCCCAGAAATTGCAGATCATCGCCACGGC‑GCCCCCCCTGCCACCCACCGGGAACTAATTCCTGCTCAACACGGGCAATCCCATCCCGCAAAAGCaac                                                                                   >  1:379459/1‑236 (MQ=255)
                                                                                                                                                               ccTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGa                                                                                                                                                                                          <  1:306837/131‑1 (MQ=255)
                                                                                                                                                               ccTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑AACCTGGCCTAAATC‑CTCCCC‑GAGATTTACGCCAATTCCAACCCCCTCATTTGCCCCCAGCAATTGCAGATCCACCCCACGGC‑GCCCCCCCCTCCCGCCCACCGGAACTAATGCCTCCTCAAAACCGCCAATAACACCGCGAAAAgaaaaaccccac                                                                             >  2:57931/1‑228 (MQ=255)
                                                                                                                                                                cTGTACCTGCTCACAGACGATACGACCAATTTTGATTGCGTCGTAAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTAGA‑GCCCTGGCCTAAATC‑CGCCGC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAgcg                                                                            <  1:539191/239‑1 (MQ=255)
                                                                                                                                                                cTGTACCTGCCCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAgcg                                                                            <  1:812935/239‑1 (MQ=255)
                                                                                                                                                                      ctccccccccacgATTCGCCAATTTTTGATTGCGTCGACAAAAGAGAT‑CACGGAAA‑GA‑TTGGC‑GGTAAT‑GTTATA‑CACCTGGCTTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTa                                                                      <  2:643109/230‑1 (MQ=255)
                                                                                                                                                                             cgcacgATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCCCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAAAGCCTTCTCAAAACGGTCAATCACATCGCGAAAAGGCATCCCCAGCGGCGGTAACCGgc                                                                >  1:32534/1‑238 (MQ=255)
                                                                                                                                                                                        aCCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCCGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCACCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCCCGAAAAAGCACCCGCACCGGCGGGAACGGgggccggccccaaa                                                    >  1:472965/1‑229 (MQ=255)
                                                                                                                                                                                          aaatttttttttgtgcgtcAAAAGAGGC‑CAAAGAAA‑GA‑TTTGG‑GGTAAC‑TTTAAA‑CACCGGGCCGAAATCGCGCC‑C‑GAATTTAACCCCCACCGCAGCCGCCTGATTGCCCCACAGAAATTGCAGACCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCGCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAAtt                                                  <  2:779559/224‑1 (MQ=255)
                                                                                                                                                                                                  tttGGGTCGACAAAAGAGAC‑CACAGAGA‑GG‑TTGGC‑GGAAAT‑GTTTGA‑CCGCCGTCCATAACC‑CGCCTCGGAGA‑TATCGCCAATCGCAACCGCCGGATTGACCCACAGAAATTGCAGAACATCGCACCGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCGGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATGTCCTGATg                                          <  2:780649/238‑1 (MQ=255)
                                                                                                                                                                                                  atgtgctcctcAAAAGGAAC‑CACAGAGA‑GA‑TTGGG‑GGAAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTACCGCAAATCGCAACCGCCTGATTGACCCACAGAAAGTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATg                                          <  1:1145106/230‑1 (MQ=255)
                                                                                                                                                                                                  aTTGCGTCGTCAAAAGAGAT‑CCAAGAGG‑GT‑TTGGC‑GGTAAT‑TTTAAA‑GACCTGGCCTAAATC‑CGCCTC‑GGGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATg                                          <  1:617902/239‑1 (MQ=255)
                                                                                                                                                                                                         ataaagggacagt‑CACAGAGG‑GT‑TTGGC‑GGTAAG‑GTTATA‑ACCCGGGCCTAAACC‑CGCCTC‑GAGATTAACCCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATa                                   <  1:764487/227‑1 (MQ=255)
                                                                                                                                                                                                          gtcAAAAGAGAT‑CACAGAGA‑GT‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCaca                                   <  1:923583/238‑3 (MQ=255)
                                                                                                                                                                                                                 aaagggatcgagaggggtttGGC‑GTTAAT‑GTTATA‑GCCCTGGCCTAAATC‑CGCCTG‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGt                              <  1:348622/221‑1 (MQ=255)
                                                                                                                                                                                                                        tcacgggggatttGGG‑GGTAAT‑TTTATA‑CGCCTGGCAAAAACC‑CCCCTC‑GGAGTTACCGCAAATCCCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGc                         <  2:379459/226‑1 (MQ=255)
                                                                                                                                                                                                                           aaggagaaggggC‑GGTAAT‑GGAAAG‑AACCTGGCATAAATG‑CCCCGG‑TAGATTATCACCAAACGCAAACCACCGAGTGACCCACAGAAATTGCAGATCACCGCCACGGC‑GCTCCACCGGCCAGCCAGCGGAAACTAATGCCTGATCAACACGGTCAATCACAGCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCt                    <  2:957405/229‑1 (MQ=255)
                                                                                                                                                                                                                                              t‑GTTTTA‑GCCCTGGCCTAAGTC‑CGCCTC‑GAGATTATCGCCAATCGCACCCGCCTGATTGACCCACAGAAATTGCAGATCAGCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCTTCATGTTGCAGCTCTGcc  <  1:714125/239‑1 (MQ=255)
                                                                                                                                                                                                                                              t‑GTTATA‑GCCCTGGCCTAAACC‑CGCCCC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCCCCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCTTCATGTTGCAGCTCTGcc  <  1:1156964/239‑1 (MQ=255)
                                                                                                                                                                                                                                               gtttAAA‑GCCCGGGCCTAAATC‑CGCCTC‑GAGATTACCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCTTCATGTTGCAGCTCTGcc  <  1:659802/237‑1 (MQ=255)
                                                                                                                                                                                                                                                gTTAAA‑GCCCTGGCCTAACTC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCTTCATGTTGCAGCTCTGcc  <  1:998132/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttATA‑GCCCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATg                                          <  1:869029/197‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
TCGACGCCATAG‑CCCATATGCACGACCTGGAAACCGTTGTCATAAGCGCCAGAACGCGGCGCATTGATACCAATAAAGCGCCCTTCTTTGCTTTCTGCGGCATCGAGAATGCGTTGACTGTCGGTTAACAGGCCACCAAA‑TGTAATGACCCGACCGACCTGTACCTGCTCACGCACGATACGACCAATTTTGATTGCGTCGTCAAAAGAGAT‑CACAGAGA‑GA‑TTGGC‑GGTAAT‑GTTATA‑GACCTGGCCTAAATC‑CGCCTC‑GAGATTATCGCCAATCGCAACCGCCTGATTGACCCACAGAAATTGCAGATCATCGCCACGGC‑GCTCCACCTGCCAGCCAGCGGAAACTAATGCCTGCTCAACACGGTCAATCACATCGCGAAAAGGCATCCGCAGCGGCGGTAACGGGCGCGTGCCTAAAATTTCCTGATGCCCCATAAAGGTATCGCCACCTTCATGTTGCAGCTCTGCC  >  NC_000913/3509120‑3509590

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑