Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,533,038 G→C 6.2% R74P (CGT→CCT) ‡ pck → phosphoenolpyruvate carboxykinase [ATP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,533,0380GC 6.2% 582.5 / 21.2 164R74P (CGT→CCT) ‡pckphosphoenolpyruvate carboxykinase [ATP]
Reads supporting (aligned to +/- strand):  ref base G (72/81);  new base C (2/8);  total (75/89)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 164 total aligned reads displayed.

ACCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGTT  >  NC_000913/3532894‑3533177
                                                                                                                                                |                                                                                                                                           
aCCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGAc                                                                                                                                         >  2:992941/1‑149 (MQ=255)
 cccAAGCTATGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACg                                                                                                                                        >  2:187341/1‑149 (MQ=255)
 cccAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACg                                                                                                                                        >  3:1048788/1‑149 (MQ=255)
  ccAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGa                                                                                                                                       <  3:385832/149‑1 (MQ=255)
    aaGCTACGACCAGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATa                                                                                                                                     <  1:357604/149‑1 (MQ=255)
      gCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATata                                                                                                                                   <  1:703317/149‑3 (MQ=255)
      gCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATata                                                                                                                                   <  1:353559/149‑3 (MQ=255)
      gCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTACGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATata                                                                                                                                   <  2:1279120/149‑3 (MQ=255)
       cTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCACGGACGATatac                                                                                                                                  <  2:1405662/149‑4 (MQ=255)
               tgctgTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGAt                                                                                                                          <  2:106197/149‑1 (MQ=255)
                gctgTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATa                                                                                                                         >  3:616204/1‑149 (MQ=255)
                 ctgTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGACGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATAc                                                                                                                        <  4:1199722/149‑1 (MQ=255)
                     aTCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTc                                                                                                                    <  4:1418900/149‑1 (MQ=255)
                      tCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCt                                                                                                                   >  1:65800/1‑149 (MQ=255)
                       cAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtg                                                                                                                  <  1:992941/149‑1 (MQ=255)
                         ggAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCt                                                                                                                   <  1:1009786/146‑1 (MQ=255)
                          gAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATAc                                                                                                                        >  1:439979/1‑140 (MQ=255)
                          gAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggt                                                                                                                <  1:1462663/148‑1 (MQ=255)
                           aaGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggtgg                                                                                                              <  3:648515/149‑1 (MQ=255)
                           aaGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggtgg                                                                                                              <  4:1048788/149‑1 (MQ=255)
                           aaGAGCTCGATACGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggtgg                                                                                                              >  3:836257/1‑149 (MQ=255)
                            agagCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTggg                                                                                                             >  1:975501/1‑149 (MQ=255)
                             gagCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGc                                                                                                            <  4:448459/149‑1 (MQ=255)
                                   gATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTc                                                                                                                    <  2:1366551/135‑1 (MQ=255)
                                   gATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTc                                                                                                                    >  1:1366551/1‑135 (MQ=255)
                                     tCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    <  1:932971/149‑1 (MQ=255)
                                     tCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    <  4:617295/149‑1 (MQ=255)
                                       cGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCt                                                                                                                   >  2:677347/1‑132 (MQ=255)
                                       cGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCt                                                                                                                   <  1:677347/132‑1 (MQ=255)
                                           ccTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaag                                                                                              <  2:85406/149‑1 (MQ=255)
                                               aCAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCAGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTaa                                                                                          >  1:114826/1‑149 (MQ=255)
                                                    ttATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGaa                                                                                       >  2:1328091/1‑147 (MQ=255)
                                                     tATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGa                                                                                    <  3:287226/149‑1 (MQ=255)
                                                                  tgtTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaag                                                                                              <  3:1464557/126‑1 (MQ=255)
                                                                  tgtTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaag                                                                                              >  4:1464557/1‑126 (MQ=255)
                                                                  tgtTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    <  2:439619/120‑1 (MQ=255)
                                                                    ttAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaa                                                                                                     <  2:510152/117‑1 (MQ=255)
                                                                    ttAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACctctctc                                                                     >  1:164793/1‑149 (MQ=255)
                                                                     tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTc                                                                                                                    <  4:131775/101‑1 (MQ=255)
                                                                     tAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTc                                                                                                                    >  3:131775/1‑101 (MQ=255)
                                                                      aaCTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaag                                                                                              <  2:481950/122‑1 (MQ=255)
                                                                      aaCTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCtctctctc                                                                   >  4:370535/1‑149 (MQ=255)
                                                                       aCTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATAc                                                                                                                        <  2:1445177/95‑1 (MQ=255)
                                                                       aCTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATAc                                                                                                                        >  1:1445177/1‑95 (MQ=255)
                                                                       aCTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACAAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTcc                                                                  >  4:655640/1‑149 (MQ=255)
                                                                       aCTAATCTGGGTGCCGTTGCCGTCGATAACGAGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACACGCGATACGTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTcc                                                                  >  4:900181/1‑149 (MQ=255)
                                                                         tAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    <  2:262077/113‑1 (MQ=255)
                                                                         tAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    >  1:262077/1‑113 (MQ=255)
                                                                          aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGaa                                                                                       >  1:1324602/1‑125 (MQ=255)
                                                                          aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGacaa                                                                                 <  4:1223477/131‑1 (MQ=255)
                                                                          aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGacaa                                                                                 >  3:1223477/1‑131 (MQ=255)
                                                                          aaTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACaaa                                                                             <  1:513549/135‑1 (MQ=255)
                                                                           aTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaa                                                              >  3:856962/1‑149 (MQ=255)
                                                                           aTCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGaa                                                              >  2:1195633/1‑149 (MQ=255)
                                                                                ggTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTg                                                         <  3:370535/149‑1 (MQ=255)
                                                                                 gTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                        <  2:975501/149‑1 (MQ=255)
                                                                                 gTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                        <  3:655640/149‑1 (MQ=255)
                                                                                      gttgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGAGAGGTATATCGTCCATGACGATACCACTCGCGATACTTTCTGGCGGGCAGACAATGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCTGAAACATGgcaaca                                                   >  4:902895/1‑149 (MQ=255)
                                                                                       ttgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCtggt                                                                                                                <  4:1023513/87‑1 (MQ=255)
                                                                                       ttgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAAc                                                                            >  3:175482/1‑123 (MQ=255)
                                                                                       ttgccgtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAAc                                                                            <  4:175482/123‑1 (MQ=255)
                                                                                            gtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                        >  2:646824/1‑138 (MQ=255)
                                                                                            gtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTgg                                                        <  1:646824/138‑1 (MQ=255)
                                                                                            gtCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGaa                                             >  4:813858/1‑149 (MQ=255)
                                                                                               gATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAgg                                          <  4:856962/149‑1 (MQ=255)
                                                                                                 tACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGa                                                                                    <  2:870685/105‑1 (MQ=255)
                                                                                                 tACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGa                                                                                    >  1:870685/1‑105 (MQ=255)
                                                                                                 tACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAAc                                                                            >  4:466670/1‑113 (MQ=255)
                                                                                                 tACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAAc                                                                            <  3:466670/113‑1 (MQ=255)
                                                                                                  aCCAGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCt                                       <  2:114826/149‑1 (MQ=255)
                                                                                                      ggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTggg                                                                                                             >  3:283521/1‑75 (MQ=255)
                                                                                                      ggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTggg                                                                                                             <  4:283521/75‑1 (MQ=255)
                                                                                                      ggATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGa                                                                                    <  1:1405598/100‑1 (MQ=255)
                                                                                                           ttCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAAc                                                            >  2:483839/1‑119 (MQ=255)
                                                                                                           ttCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAAc                                                            <  1:483839/119‑1 (MQ=255)
                                                                                                              aCCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCa                           >  3:307919/1‑149 (MQ=255)
                                                                                                                      ttCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCg                   >  3:139/1‑149 (MQ=255)
                                                                                                                      ttCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCg                   >  2:212164/1‑149 (MQ=255)
                                                                                                                       tCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAGCCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCgg                  <  4:307919/149‑1 (MQ=255)
                                                                                                                       tCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCgg                  <  2:164793/149‑1 (MQ=255)
                                                                                                                           cAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGt                                    >  3:621635/1‑127 (MQ=255)
                                                                                                                           cAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGt                                    <  4:621635/127‑1 (MQ=255)
                                                                                                                           cAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGAcc                                >  2:571927/1‑131 (MQ=255)
                                                                                                                           cAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTcc                    <  2:1083365/143‑1 (MQ=255)
                                                                                                                           cAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCaaa              <  3:900181/149‑1 (MQ=255)
                                                                                                                             aaaGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTaa                                                                                          <  2:1446544/71‑1 (MQ=255)
                                                                                                                                gATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCt         >  2:322491/1‑149 (MQ=255)
                                                                                                                                    aGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTc     <  1:212164/149‑1 (MQ=255)
                                                                                                                                     gTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaaggcaaag                                                                                              <  1:491362/59‑1 (MQ=255)
                                                                                                                                     gTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCa                           >  4:966861/1‑126 (MQ=255)
                                                                                                                                     gTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCa                           <  3:966861/126‑1 (MQ=255)
                                                                                                                                     gTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCg    >  3:580478/1‑149 (MQ=255)
                                                                                                                                      tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTc     <  3:932005/147‑1 (MQ=255)
                                                                                                                                      tataTCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGt   <  1:322491/149‑1 (MQ=255)
                                                                                                                                       gtatCGTCACGGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCg                   >  3:684483/2‑132 (MQ=255)
                                                                                                                                       gtatCGTCACGGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCg                   <  4:684483/131‑1 (MQ=255)
                                                                                                                                       atatCGTCCGTGACGATACCACTCGCGATACTTTGTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCCGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGtt  <  1:1195633/149‑1 (MQ=255)
                                                                                                                                       atatCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGtt  <  4:580478/149‑1 (MQ=255)
                                                                                                                                            gTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    >  2:1010458/1‑46 (MQ=255)
                                                                                                                                            gTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGAcaaag                                                                                                    <  1:1010458/46‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
ACCCAAGCTACGACCTGCTGTATCAGGAAGAGCTCGATCCGAGCCTGACAGGTTATGAGCGCGGGGTGTTAACTAATCTGGGTGCCGTTGCCGTCGATACCGGGATCTTCACCGGTCGTTCACCAAAAGATAAGTATATCGTCCGTGACGATACCACTCGCGATACTTTCTGGTGGGCAGACAAAGGCAAAGGTAAGAACGACAACAAACCTCTCTCTCCGGAAACCTGGCAGCATCTGAAAGGCCTGGTGACCAGGCAGCTTTCCGGCAAACGTCTGTTCGTT  >  NC_000913/3532894‑3533177

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