Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,063,045 C→G 5.4% A20G (GCC→GGC)  argK → membrane ATPase/protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,063,0450CG 5.4% 318.4 / 6.3 90A20G (GCC→GGC) argKmembrane ATPase/protein kinase
Reads supporting (aligned to +/- strand):  ref base C (50/35);  new base G (2/3);  total (52/38)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGCACCCCCGGCGCGGGGAAAAGTAC  >  NC_000913/3062898‑3063192
                                                                                                                                                   |                                                                                                                                                   
gCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGcc                                                                                                                                                    >  1:501015/1‑149 (MQ=255)
  aaGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCcac                                                                                                                                                  >  2:709205/1‑149 (MQ=255)
      gcgcGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTc                                                                                                                                              <  2:619616/149‑1 (MQ=255)
        gcggcgTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGc                                                                                                                                            >  4:884021/1‑149 (MQ=255)
            cgTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAg                                                                                                                                        <  2:540437/149‑1 (MQ=255)
               ggcggcGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGcc                                                                                                                                     <  3:384034/149‑1 (MQ=255)
                 cggcgATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCAt                                                                                                                                   <  2:172950/149‑1 (MQ=255)
                      aTTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGc                                                                                                                              <  1:1138754/149‑1 (MQ=255)
                       tttATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGGCACACTCGCCCAGGCCATGACGCt                                                                                                                             <  1:975511/149‑1 (MQ=255)
                        ttATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAAACTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCtg                                                                                                                            <  3:946411/149‑1 (MQ=255)
                         tATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCtgg                                                                                                                           <  3:884021/149‑1 (MQ=255)
                           tGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCtggtg                                                                                                                         <  1:820434/149‑1 (MQ=255)
                               ccAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGaaa                                                                                                                     <  4:1118896/149‑1 (MQ=255)
                                 aGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGCGAGTGTGGCACGCTCAccct                                                                                                                                         <  1:1086889/127‑2 (MQ=255)
                                 aGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGCGAGTGTGGCACGCTCAccct                                                                                                                                         >  2:1086889/1‑126 (MQ=255)
                                  ggTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATTAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCATGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  <  2:856302/149‑1 (MQ=255)
                                  ggTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  <  2:451937/149‑1 (MQ=255)
                                  ggTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  <  1:543509/149‑1 (MQ=255)
                                          tATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  >  1:857340/1‑141 (MQ=255)
                                          tATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  <  2:857340/141‑1 (MQ=255)
                                           aTGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCcgc                                                                                                         <  2:351135/149‑1 (MQ=255)
                                                 gACAGTGAGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCATCGCTGGCATAAAGGATTGGCCGCTTACGTCAGAGTGAGCGGGCGACACTCGCGCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGtcatc                                                                                                   >  2:82465/1‑149 (MQ=255)
                                                     gtgtGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGc                                                                                               <  2:320452/149‑1 (MQ=255)
                                                        tgCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACt                                                                                            >  4:163600/1‑149 (MQ=255)
                                                          cgcgACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTaa                                                                                          <  2:501015/149‑1 (MQ=255)
                                                             gACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTa                                                                                       >  3:856280/1‑149 (MQ=255)
                                                                 tACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCtggtg                                                                                                                         <  1:511733/111‑1 (MQ=255)
                                                                 tACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCtggtg                                                                                                                         >  2:511733/1‑111 (MQ=255)
                                                                   cTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgc                                                                                 >  2:128113/1‑149 (MQ=255)
                                                                   cTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgc                                                                                 <  1:932407/149‑1 (MQ=255)
                                                                   cTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgc                                                                                 >  2:844021/1‑149 (MQ=255)
                                                                    tGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgct                                                                                >  1:597586/1‑149 (MQ=255)
                                                                     gAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGc                                                                                                                                      <  4:703435/94‑1 (MQ=255)
                                                                     gAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGc                                                                                                                                      >  3:703435/1‑94 (MQ=255)
                                                                     gAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctg                                                                               >  3:644196/1‑149 (MQ=255)
                                                                      aaTCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctgc                                                                              >  4:229700/1‑149 (MQ=255)
                                                                       aTCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctgct                                                                             >  1:887269/1‑149 (MQ=255)
                                                                        tCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGaaa                                                                                                                     <  4:159784/108‑1 (MQ=255)
                                                                        tCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGaaa                                                                                                                     >  3:159784/1‑108 (MQ=255)
                                                                        tCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCtt                                                                            >  2:378613/1‑149 (MQ=255)
                                                                        tCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCtt                                                                            >  4:926354/1‑149 (MQ=255)
                                                                          tGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTTCTTGa                                                                          <  2:524410/149‑1 (MQ=255)
                                                                           gATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGAt                                                                         <  2:1144842/149‑1 (MQ=255)
                                                                           gATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGAt                                                                         >  2:437565/1‑149 (MQ=255)
                                                                           gATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGAt                                                                         >  2:421563/1‑149 (MQ=255)
                                                                            aTAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAg                                                                                                 >  4:434017/1‑124 (MQ=255)
                                                                            aTAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAg                                                                                                 <  3:434017/124‑1 (MQ=255)
                                                                             tAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGc                                                                       >  4:358796/1‑149 (MQ=255)
                                                                             tAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGc                                                                       >  4:914560/1‑149 (MQ=255)
                                                                              aaGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgc                                                                                 >  4:1052295/1‑138 (MQ=255)
                                                                              aaGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgc                                                                                 <  3:1052295/138‑1 (MQ=255)
                                                                               aGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCa                                                                     >  3:202514/1‑149 (MQ=255)
                                                                               aGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCa                                                                     <  3:358796/149‑1 (MQ=255)
                                                                                gCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCAt                                                                    >  1:685549/1‑149 (MQ=255)
                                                                                       tcatGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  >  1:521874/1‑96 (MQ=255)
                                                                                       tcatGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGcc                                                                                                                  <  2:521874/96‑1 (MQ=255)
                                                                                       tcatGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGCCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCg                                                             >  3:921186/1‑149 (MQ=255)
                                                                                          tGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGcc                                                                      >  2:31379/1‑137 (MQ=255)
                                                                                          tGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGcc                                                                      <  1:31379/137‑1 (MQ=255)
                                                                                            aTTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTg                                                        >  3:851463/1‑149 (MQ=255)
                                                                                                 tGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTa                                                   >  3:570924/1‑149 (MQ=255)
                                                                                                 tGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTa                                                   >  3:935505/1‑149 (MQ=255)
                                                                                                  gAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTaa                                                  >  1:1116986/1‑149 (MQ=255)
                                                                                                     gCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAacac                                               >  1:1127691/1‑149 (MQ=255)
                                                                                                      ccACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACAcc                                              >  2:745167/1‑149 (MQ=255)
                                                                                                              gCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACt                                      <  2:1116986/149‑1 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTTAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTgg                                    >  4:273984/1‑149 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTgg                                    >  2:252055/1‑149 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTgg                                    <  1:844021/149‑1 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTgg                                    >  1:840708/1‑149 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTgg                                    >  1:704548/1‑149 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGGAACACCCTGCGACTgg                                    >  1:969920/1‑149 (MQ=255)
                                                                                                                aGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGCCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTgg                                    <  4:921186/149‑1 (MQ=255)
                                                                                                                 gAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACgggg                                   >  3:133383/1‑149 (MQ=255)
                                                                                                                  aaaGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGc                                  >  4:171140/1‑149 (MQ=255)
                                                                                                                     gTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGTTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGtt                               <  4:935505/149‑1 (MQ=255)
                                                                                                                      tATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTa                              >  3:559145/1‑149 (MQ=255)
                                                                                                                        ttCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTAcc                            >  1:229726/1‑149 (MQ=255)
                                                                                                                         tCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCg                           >  4:314669/1‑149 (MQ=255)
                                                                                                                           gccgcTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGc                         >  2:762274/1‑149 (MQ=255)
                                                                                                                                    gTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTc                                                                                                               <  3:646318/54‑1 (MQ=255)
                                                                                                                                    gTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTc                                                                                                               >  4:646318/1‑54 (MQ=255)
                                                                                                                                    gTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctgc                                                                              <  3:62534/87‑1 (MQ=255)
                                                                                                                                    gTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAgctgc                                                                              >  4:62534/1‑87 (MQ=255)
                                                                                                                                    gTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTCGGCGTTACCGGCACCCCCGgc                >  4:830587/1‑149 (MQ=255)
                                                                                                                                     tCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAaca                                                <  2:354186/116‑1 (MQ=255)
                                                                                                                                     tCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAaca                                                >  1:354186/1‑116 (MQ=255)
                                                                                                                                       tgggcgAGTGTGGCACGCTCACCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGCACCCCCGGcgcg             >  4:320977/6‑149 (MQ=255)
                                                                                                                                        gggCGAGTGTGGCACGCTCACCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGCACCCCCGGCGCgg            <  3:320977/149‑1 (MQ=255)
                                                                                                                                         ggTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGCACCCCCGGCGCggg           <  2:1096638/149‑1 (MQ=255)
                                                                                                                                                  gCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGCACCCCCGGCGCGGGGAAAAGTAc  >  3:649967/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTCGACAGTGTGCGCGACGTACTGAATCTGATAAGCCAGCATCATGATTAATGAAGCCACGCTGGCAGAAAGTATTCGCCGCTTACGTCAGGGTGAGCGTGCCACACTCGCCCAGGCCATGACGCTGGTGGAAAGCCGTCACCCGCGTCATCAGGCACTAAGTACGCAGCTGCTTGATGCCATTATGCCGTACTGCGGTAACACCCTGCGACTGGGCGTTACCGGCACCCCCGGCGCGGGGAAAAGTAC  >  NC_000913/3062898‑3063192

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