Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,891,842 C→T 5.1% G465D (GGT→GAT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,891,8420CT 5.1% 388.2 / 7.3 118G465D (GGT→GAT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base C (60/52);  new base T (2/4);  total (62/56)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

Only 100 of 125 total aligned reads displayed.

AGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGCCTG  >  NC_000913/1891695‑1891990
                                                                                                                                                   |                                                                                                                                                    
agCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGAcc                                                                                                                                                     <  2:1019356/149‑1 (MQ=255)
agCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGAcc                                                                                                                                                     <  1:1077174/149‑1 (MQ=255)
     aTTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttgt                                                                                                                                                >  1:403212/1‑149 (MQ=255)
      ttACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttgtt                                                                                                                                               <  3:423632/149‑1 (MQ=255)
      ttACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttgtt                                                                                                                                               <  4:382554/149‑1 (MQ=255)
       tACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgttt                                                                                                                                              <  2:645089/149‑2 (MQ=255)
       tACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTa                                                                                                                                              <  2:67558/149‑1 (MQ=255)
          cGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCt                                                                                                                                           >  3:385889/1‑149 (MQ=255)
              gcTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCga                                                                                                                                       <  1:741496/149‑1 (MQ=255)
              gcTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCga                                                                                                                                       <  4:83227/149‑1 (MQ=255)
              gcTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCga                                                                                                                                       <  1:1035847/149‑1 (MQ=255)
                tGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatg                                                                                                                                     >  3:645501/1‑149 (MQ=255)
                  aCAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCTGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatgat                                                                                                                                   <  3:114070/149‑1 (MQ=255)
                        tgctgcAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCg                                                                                                                             <  2:457376/149‑1 (MQ=255)
                            gcAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcagc                                                                                                                         >  1:1069904/1‑149 (MQ=255)
                             cAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcagca                                                                                                                        >  2:200736/1‑149 (MQ=255)
                              aaCAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCAt                                                                                                                       >  3:617794/1‑149 (MQ=255)
                               aCAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATc                                                                                                                                  <  2:954738/137‑1 (MQ=255)
                                  gTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGcc                                                                                                                   <  3:984253/149‑1 (MQ=255)
                                  gTTGCCCTTTGCTGGTTACCCCCTGCATCATTACAGGAAGCGTCATGGTAGCGCGTAACTGCTCGGCCAGATCGCGCATCATGGCGTTGGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGcc                                                                                                                   <  2:882158/149‑1 (MQ=255)
                                   ttGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCg                                                                                                                  >  2:594903/1‑149 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGAc                                                                                                                                                      >  1:891763/1‑109 (MQ=255)
                                        cTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttta                                                                                                                                            >  2:28304/1‑115 (MQ=255)
                                         tttGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACt                                                                                                            >  1:341303/1‑149 (MQ=255)
                                            gCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGc                                                                                                         <  4:385889/149‑1 (MQ=255)
                                              tGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGa                                                                                                       >  4:915805/1‑149 (MQ=255)
                                              tGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGa                                                                                                       <  2:249897/149‑1 (MQ=255)
                                              tGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGa                                                                                                       >  3:1063819/1‑149 (MQ=255)
                                                gTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGc                                                                                                     >  1:358613/1‑149 (MQ=255)
                                                gTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGc                                                                                                     >  3:285731/1‑149 (MQ=255)
                                                gTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGc                                                                                                     >  3:684930/1‑149 (MQ=255)
                                                 tttCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCa                                                                                                    <  3:192847/149‑1 (MQ=255)
                                                 tttCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCa                                                                                                    <  3:75386/149‑1 (MQ=255)
                                                    cccccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGaa                                                                                                 <  2:839317/149‑1 (MQ=255)
                                                      cccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAcc                                                                                               >  3:75802/1‑149 (MQ=255)
                                                       ccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACct                                                                                              >  3:956126/1‑149 (MQ=255)
                                                       ccTGCAACAATACAGGAAGCGTCATAGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACCACCGTTGTTAGCTTCGATGATCAGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACct                                                                                              <  4:195781/149‑1 (MQ=255)
                                                             acaaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggctgg                                                                                        <  1:228617/149‑1 (MQ=255)
                                                              caaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggt                                                                                       <  2:142434/149‑1 (MQ=255)
                                                               aaTCCAGGAAGCGTCATGGTGGCGCGGAACTTGTCGGCCAGACCGGGCAGCATGGCGTGCGCGGTACACACCATAAAGCAACAACCGATATTAGCTGCGCTGATCGGTCGCAGCATTGCCGCTAACAGGCGAGCAGAGCCTGGCgggtt                                                                                      <  2:984485/149‑1 (MQ=255)
                                                               aaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggtt                                                                                      <  4:51771/149‑1 (MQ=255)
                                                                 tACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGa                                                                                                       >  4:730818/1‑130 (MQ=255)
                                                                 tACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACct                                                                                              >  4:91193/1‑139 (MQ=255)
                                                                    aGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGttt                                                                                 >  3:1047825/1‑149 (MQ=255)
                                                                      gAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGc                                                                               >  2:1000129/1‑149 (MQ=255)
                                                                       aaGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCg                                                                              <  1:455995/149‑1 (MQ=255)
                                                                         gcatCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGc                                                                            <  4:1019156/146‑1 (MQ=255)
                                                                           gTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAg                                                                          >  4:55238/1‑149 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGa                                                                         >  4:336148/1‑149 (MQ=255)
                                                                              aTGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGAtt                                                                       >  2:47419/1‑149 (MQ=255)
                                                                                   aGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgca                                                                  <  4:645501/149‑1 (MQ=255)
                                                                                   aGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgca                                                                  <  1:100387/149‑1 (MQ=255)
                                                                                      gcgGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGAGCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgcagca                                                               >  4:857716/1‑149 (MQ=255)
                                                                                       cgGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcac                                                              >  1:459174/1‑149 (MQ=255)
                                                                                         gAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacac                                                            >  1:128152/1‑149 (MQ=255)
                                                                                         gAACTGCTAGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcacac                                                            >  3:997674/1‑149 (MQ=255)
                                                                                                    ccAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGcc                                                 <  2:112140/149‑1 (MQ=255)
                                                                                                       gATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttg                                                                                     >  1:1036509/1‑110 (MQ=255)
                                                                                                       gATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttg                                                                                     <  2:1036509/110‑1 (MQ=255)
                                                                                                        aTCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAAGTGGCGAGCAGAACCTGGCttgttg                                                                                     <  1:938711/109‑1 (MQ=255)
                                                                                                        aTCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttg                                                                                     >  2:938711/1‑109 (MQ=255)
                                                                                                           gcgcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTggttggtt                                                                                  >  3:574582/1‑109 (MQ=255)
                                                                                                            cgcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCt                                         <  1:624520/149‑1 (MQ=255)
                                                                                                             gcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTg                                        <  4:617794/149‑1 (MQ=255)
                                                                                                             gcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTg                                        >  2:828567/1‑149 (MQ=255)
                                                                                                                tcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTaa                                     <  1:639215/149‑1 (MQ=255)
                                                                                                                tcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTaa                                     >  3:867157/1‑149 (MQ=255)
                                                                                                                tcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCAGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTaa                                     >  1:532437/1‑149 (MQ=255)
                                                                                                                  atGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATc                                   >  4:773783/1‑149 (MQ=255)
                                                                                                                     gCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaaa                                <  4:308542/149‑1 (MQ=255)
                                                                                                                       gTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtg                              <  3:857716/149‑1 (MQ=255)
                                                                                                                          cGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggc                           <  2:358613/149‑1 (MQ=255)
                                                                                                                           gAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggct                          >  2:736027/1‑149 (MQ=255)
                                                                                                                               tACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCACCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggctgggt                      >  3:96594/1‑148 (MQ=255)
                                                                                                                                 cAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggctggg                       >  4:265036/1‑146 (MQ=255)
                                                                                                                                 cAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggctggg                       <  3:265036/146‑1 (MQ=255)
                                                                                                                                  aaaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCTGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAAc                   >  2:649454/1‑149 (MQ=255)
                                                                                                                                  aaaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAAc                   <  4:1009140/149‑1 (MQ=255)
                                                                                                                                  aaaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTAGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGAGCCAc                   >  4:778341/1‑149 (MQ=255)
                                                                                                                                   aaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtg                              <  2:1057837/137‑1 (MQ=255)
                                                                                                                                   aaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACa                  >  3:627547/1‑149 (MQ=255)
                                                                                                                                   aaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACa                  <  1:232809/149‑1 (MQ=255)
                                                                                                                                       agctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCt                                         <  4:581798/118‑1 (MQ=255)
                                                                                                                                       agctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCt                                         >  3:581798/5‑122 (MQ=255)
                                                                                                                                        tAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCATATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAAc             >  2:372607/1‑149 (MQ=255)
                                                                                                                                        tAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGAtt                                                                       <  1:107468/91‑1 (MQ=255)
                                                                                                                                        tAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGAtt                                                                       >  2:107468/1‑91 (MQ=255)
                                                                                                                                        tAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAAc             <  4:994308/149‑1 (MQ=255)
                                                                                                                                          aaaCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGa           >  1:725850/1‑149 (MQ=255)
                                                                                                                                           aacaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGAc          >  1:1093208/1‑149 (MQ=255)
                                                                                                                                            acaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACt         >  1:1024833/1‑149 (MQ=255)
                                                                                                                                              aacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCg       >  4:1080148/1‑149 (MQ=255)
                                                                                                                                              aacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCg       <  3:1080148/149‑1 (MQ=255)
                                                                                                                                                cGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGGTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCGGTAATCAAATGGGCTGGGCAACACCAACGACTCGGc     <  4:913698/149‑1 (MQ=255)
                                                                                                                                                 gACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTGGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGcc    >  4:1070398/1‑149 (MQ=255)
                                                                                                                                                  aCCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGCCt   >  3:155666/1‑149 (MQ=255)
                                                                                                                                                  aCCGTTGTTAGCTTCGATGAGCGGTCGCAGCATTGCCGCTCACTGGCGAGAAGAACCTGGATGGTTGGTGTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTGATCAAATGGGCGGGGCAACAGCAACGACTCGGcc  >  2:233949/1‑148 (MQ=255)
                                                                                                                                                   ccGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggctgggc                      >  2:230071/1‑129 (MQ=255)
                                                                                                                                                   ccGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggctgggc                      <  1:230071/129‑1 (MQ=255)
                                                                                                                                                   ccGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGCCTg  >  3:288076/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGACTCGGCCTG  >  NC_000913/1891695‑1891990

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: