Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,243,472 T→C 5.1% M61V (ATG→GTG)  folE ← GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,243,4720TC 5.1% 248.7 / 2.8 79M61V (ATG→GTG) folEGTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

GTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGA  >  NC_000913/2243330‑2243620
                                                                                                                                              |                                                                                                                                                    
gTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACt                                                                                                                                                >  1:370465/1‑149 (MQ=255)
gTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACt                                                                                                                                                <  8:246005/149‑1 (MQ=255)
 tCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTg                                                                                                                                               >  8:54487/1‑149 (MQ=255)
   agagTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCAt                                                                                                                                                      <  6:300380/140‑1 (MQ=255)
   agagTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCAt                                                                                                                                                      >  5:300380/1‑140 (MQ=255)
    gagTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtca                                                                                                                                            >  6:219749/1‑149 (MQ=255)
     agTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcat                                                                                                                                           <  1:100829/149‑1 (MQ=255)
     agTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcat                                                                                                                                           <  3:156701/149‑1 (MQ=255)
           atCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtc                                                                                                                                             <  8:62810/141‑1 (MQ=255)
           atCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtc                                                                                                                                             >  7:62810/1‑141 (MQ=255)
                gcACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTc                                                                                                                                 >  6:110644/1‑148 (MQ=255)
                gcACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCg                                                                                                                                <  7:368973/149‑1 (MQ=255)
                   cGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCgaga                                                                                                                             <  7:54487/149‑1 (MQ=255)
                       gACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGAtt                                                                                                                           >  4:5739/1‑147 (MQ=255)
                          cATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTcag                                                                                                                        >  5:235603/1‑147 (MQ=255)
                          cATTTCATCGAAATTCATTTTGGTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATAGACATACATTTTATCGATGAGATGCGGCGTTTCCAGCAAACTGACATAAGCCAGGTCGACATTcag                                                                                                                        >  6:272479/1‑147 (MQ=255)
                            tttCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCagcagc                                                                                                                    >  2:186952/1‑149 (MQ=255)
                            tttCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCagcagc                                                                                                                    >  8:156307/1‑149 (MQ=255)
                            tttCATCGACCTACATATTGTTTACAATGAGGGTGATTTTAGGGAAATTGGCGTAATCAAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCagcagc                                                                                                                    >  2:220499/1‑149 (MQ=255)
                                aTCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCga                                                                                                                               >  5:275823/1‑134 (MQ=255)
                                aTCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCa                                                                                                                <  5:108980/149‑1 (MQ=255)
                                aTCGACCTTCATTTTGTTTTCAATGAGGGGGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCga                                                                                                                               <  6:275823/134‑1 (MQ=255)
                                 tCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcatc                                                                                                                                          <  7:84749/122‑1 (MQ=255)
                                 tCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcatc                                                                                                                                          >  8:84749/1‑122 (MQ=255)
                                    aCCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGAt                                                                                                            >  1:48572/1‑149 (MQ=255)
                                    aCCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGAt                                                                                                            >  6:346410/1‑149 (MQ=255)
                                      cTTCATTTTGTTTTCAATGCGGGGGATTTTCGGGAAATGGGCGTAATCCAGACCGGACAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGAttt                                                                                                          <  4:105901/149‑1 (MQ=255)
                                      cTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGa                                                                                                             <  6:263985/146‑1 (MQ=255)
                                      cTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGa                                                                                                             >  5:263985/1‑146 (MQ=255)
                                         cATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTc                                                                                                         <  7:156307/147‑1 (MQ=255)
                                          aTTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCgg                                                                                                       <  8:370282/148‑1 (MQ=255)
                                          aTTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGt                                                                                                      <  2:370465/149‑1 (MQ=255)
                                           ttttGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTc                                                                                                     <  5:295193/149‑1 (MQ=255)
                                           ttttGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTc                                                                                                     <  3:5739/149‑1 (MQ=255)
                                                  ttCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGa                                                                                               <  6:103390/148‑1 (MQ=255)
                                                    cAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGAcc                                                                                             <  7:407807/148‑1 (MQ=255)
                                                     aaTGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcatc                                                                                                                                          <  1:277618/102‑1 (MQ=255)
                                                     aaTGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcatc                                                                                                                                          >  2:277618/1‑102 (MQ=255)
                                                     aaTGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAg                                                                                           >  3:66693/1‑149 (MQ=255)
                                                     aaTGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAg                                                                                           <  6:235603/149‑1 (MQ=255)
                                                         aGGGTGATTTTCGGGAAATTGGCTTAATCCAGACCTGATAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAAt                                                                                       <  2:48572/149‑1 (MQ=255)
                                                           ggTGATTTTCGGGAAATTGGCGTAATCCAGAGCGGAGAAAATTTCGGCGGGATACAGTGTAGCGGTGCGAGGCGGCGGTTCCATCATACGGTCAGCAGCCAGGTGGAGATTAGGCAGCGGCATGATTTCGGTCTTGTGCCCGGCAAt                                                                                       >  4:215705/1‑147 (MQ=255)
                                                                   tCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCAATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCtttt                                                                              <  5:90325/148‑1 (MQ=255)
                                                                     ggtaaaTTGGCGTATTCCAGACGGGAGAAAATTTCTTCGACATACTTTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGGCGAGATTCAGCAGCTGCCTGATTTCGGTCATATGACCAGCAATAAGGCGTTTgcg                                                                           <  6:140276/146‑1 (MQ=255)
                                                                     gggAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGAtt                                                                                                                           >  4:342304/1‑101 (MQ=255)
                                                                     gggAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGAtt                                                                                                                           <  3:342304/101‑1 (MQ=255)
                                                                         aaTTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGATGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGtt                                                                       <  5:319939/149‑1 (MQ=255)
                                                                           ttGGCGTAATCCAGACCGGATAAAATTTCATCGACATACATTTTAGCGATGCTATGCGGCGTTTCCATCAAACTGGCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTc                                                                     <  5:105886/149‑1 (MQ=255)
                                                                              gCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGtt                                                                  <  1:186952/149‑1 (MQ=255)
                                                                               cGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTa                                                                 <  5:346410/149‑1 (MQ=255)
                                                                                 tAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATc                                                               <  5:110644/149‑1 (MQ=255)
                                                                                   aTCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCa                                                             >  3:26183/1‑149 (MQ=255)
                                                                                    tCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCAt                                                            <  2:1340/149‑1 (MQ=255)
                                                                                     ccAGACCGGAGAAAATTTCATCGACATACAGTTTAGCGTTGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCGGCAGCTGCATGATTTCGGTCATATGCCGCGCAGTAAGGCTGTTGCGCGTTTCGTTAtgcgtt                                                           >  2:408834/1‑144 (MQ=255)
                                                                                         aCCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCAt                                                       >  1:320644/1‑149 (MQ=255)
                                                                                         aCCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCAt                                                       >  5:109673/1‑149 (MQ=255)
                                                                                            ggAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCa                                                    >  1:6750/1‑149 (MQ=255)
                                                                                               gaAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTTTCTTCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCAcgg                                                 <  3:198768/149‑1 (MQ=255)
                                                                                                aaaaTTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCAcggg                                                >  8:412890/1‑149 (MQ=255)
                                                                                                  aaTTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTcagca                                                                                                                      >  4:228614/1‑77 (MQ=255)
                                                                                                    tttCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCCGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCAcgggcggg                                            >  5:81039/1‑149 (MQ=255)
                                                                                                    tttCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCAcgggcggg                                            >  5:86769/1‑149 (MQ=255)
                                                                                                        aTCGACATACATTTTAGCGATGCGCTGCGGCGTTTCCATCACACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCCTTTCATGCACGTGCGGGCGCa                                        >  6:325187/1‑149 (MQ=255)
                                                                                                        aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGggcgggcg                                          >  7:306655/1‑147 (MQ=255)
                                                                                                        aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCACCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCtttt                                                                              >  6:231564/1‑111 (MQ=255)
                                                                                                        aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCACCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCtttt                                                                              <  5:237278/111‑1 (MQ=255)
                                                                                                        aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCACCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCtttt                                                                              <  5:231564/111‑1 (MQ=255)
                                                                                                        aTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCACCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCtttt                                                                              >  6:237278/1‑111 (MQ=255)
                                                                                                              atacatTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAATGCTTTTGCGCGTTTCGTAATCCATTTCATGCACGGGCGGGCGCAGCGgtg                                  >  5:161447/1‑149 (MQ=255)
                                                                                                              atacatTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGgtg                                  >  2:223067/1‑149 (MQ=255)
                                                                                                               tacatTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGgtgt                                 <  5:219749/149‑1 (MQ=255)
                                                                                                                   ttttAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTc                              >  3:324912/1‑148 (MQ=255)
                                                                                                                     ttAGCGATGCGCTGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAg                           >  8:160758/1‑149 (MQ=255)
                                                                                                                     ttAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCa                            <  4:26183/148‑1 (MQ=255)
                                                                                                                     ttAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCa                            <  1:283615/148‑1 (MQ=255)
                                                                                                                          gatgcgatgcTGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATTAGGGTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTc                      <  1:220499/149‑1 (MQ=255)
                                                                                                                          gatgcgatgcGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTc                      <  1:223067/149‑1 (MQ=255)
                                                                                                                                 tgcgTCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGAATCAGCAGCTGCATGATTTCGGGCATATGAGCATCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCaa                <  5:272479/148‑1 (MQ=255)
                                                                                                                                       gTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATTTGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACTGGCGGGCGCAGCGGGGTATCCAGTCCTCGCGCAACTAACGc         >  6:396425/1‑149 (MQ=255)
                                                                                                                                         ttCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAg                           <  7:53802/129‑1 (MQ=255)
                                                                                                                                         ttCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAg                           >  8:53802/1‑129 (MQ=255)
                                                                                                                                              tCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGa  <  6:109673/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGA  >  NC_000913/2243330‑2243620

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