Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,552,182 G→C 5.1% G183R (GGC→CGC)  uxuA → mannonate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,552,1820GC 5.1% 566.0 / 17.2 156G183R (GGC→CGC) uxuAmannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (74/74);  new base C (5/3);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

Only 100 of 156 total aligned reads displayed.

CTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATG  >  NC_000913/4552036‑4552325
                                                                                                                                                  |                                                                                                                                               
cTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGgc                                                                                                                                               >  2:1411900/1‑149 (MQ=255)
     ttCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGga                                                                                                                                          >  3:522348/1‑149 (MQ=255)
     ttCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGga                                                                                                                                          >  4:49975/1‑149 (MQ=255)
         aaaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaagaa                                                                                                                                      <  2:995220/149‑1 (MQ=255)
           aTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAgg                                                                                                                                    <  3:919982/149‑1 (MQ=255)
           aTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAgg                                                                                                                                    <  4:522348/149‑1 (MQ=255)
             gCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCt                                                                                                                                  <  3:1098297/149‑1 (MQ=255)
              cATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTa                                                                                                                                 <  1:732792/149‑1 (MQ=255)
                  tCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCt                                                                                                                                  >  3:557237/1‑144 (MQ=255)
                       aaaCGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGaagaag                                                                                                                                     >  4:1081350/1‑136 (MQ=255)
                       aaaCGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGa                                                                                                                        >  3:846086/1‑149 (MQ=255)
                         aCGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGAcc                                                                                                                      <  1:1411900/149‑1 (MQ=255)
                          cGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACGccc                                                                                                                            >  2:448020/1‑142 (MQ=255)
                          cGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACGccc                                                                                                                            <  1:448020/142‑1 (MQ=255)
                           gCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAg                                                                                                                    >  3:277943/1‑149 (MQ=255)
                               aGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTcc                                                                                                                >  2:584909/1‑149 (MQ=255)
                                     ggAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           >  1:1132044/1‑148 (MQ=255)
                                     ggAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           <  2:1132044/148‑1 (MQ=255)
                                              ttACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGc                                                                                                                                   <  2:920551/115‑1 (MQ=255)
                                              ttACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGc                                                                                                                                   >  1:920551/1‑115 (MQ=255)
                                                    cGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACa                                                                                           <  1:375926/149‑1 (MQ=255)
                                                      aagaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           <  3:1175176/131‑1 (MQ=255)
                                                      aagaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           >  4:1175176/1‑131 (MQ=255)
                                                        gaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           >  2:161929/1‑129 (MQ=255)
                                                        gaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           <  1:161929/129‑1 (MQ=255)
                                                        gaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACaaa                                                                                         >  2:240793/1‑147 (MQ=255)
                                                        gaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGa                                                                                       <  3:336587/149‑1 (MQ=255)
                                                        gaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGa                                                                                       <  3:264075/149‑1 (MQ=255)
                                                         aagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaa                                                                                                                                         >  2:692437/1‑98 (MQ=255)
                                                         aagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaa                                                                                                                                         <  1:692437/98‑1 (MQ=255)
                                                          agaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGata                                                                                     >  1:966186/1‑149 (MQ=255)
                                                           gaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTc                                                                                                                          <  2:357299/111‑1 (MQ=255)
                                                           gaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGatat                                                                                    <  1:1305502/149‑1 (MQ=255)
                                                             aaTTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTa                                                                                                                                 <  4:385349/102‑1 (MQ=255)
                                                             aaTTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTa                                                                                                                                 >  3:385349/1‑102 (MQ=255)
                                                              aTTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCt                                                                                                                                  >  3:1442901/1‑100 (MQ=255)
                                                              aTTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCt                                                                                                                                  <  4:1442901/100‑1 (MQ=255)
                                                              aTTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGa                                                                                 >  1:781399/1‑149 (MQ=255)
                                                                   tCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAacac                                                                                                        >  2:1178873/1‑121 (MQ=255)
                                                                   tCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAg                                                                            >  4:548264/1‑149 (MQ=255)
                                                                   tCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagac                                                                                                                                      <  1:255643/91‑2 (MQ=255)
                                                                      ggCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGa                                                                         >  3:790290/1‑149 (MQ=255)
                                                                              aaCGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTgcgc                                                                 >  3:470065/1‑149 (MQ=255)
                                                                                 gCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           >  3:682980/1‑104 (MQ=255)
                                                                                      gCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACtt                                                         <  4:790290/149‑1 (MQ=255)
                                                                                       ccACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                        >  3:1093159/1‑149 (MQ=255)
                                                                                       ccACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                        <  4:470065/149‑1 (MQ=255)
                                                                                         aCTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGc                                                      >  1:221268/1‑149 (MQ=255)
                                                                                          cTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGcc                                                     <  4:846086/149‑1 (MQ=255)
                                                                                           tATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCg                                                    >  4:860872/1‑149 (MQ=255)
                                                                                              gAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCg                                                                          >  4:294370/1‑124 (MQ=255)
                                                                                                 cGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCACGAAAACTTTGCCGTCTTcc                                              <  2:781399/149‑1 (MQ=255)
                                                                                                  gATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGGCCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCt                                             <  4:443887/149‑1 (MQ=255)
                                                                                                   aTGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATc                                                                                   >  3:1057563/1‑110 (MQ=255)
                                                                                                     gAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                        >  4:626047/1‑135 (MQ=255)
                                                                                                      aaGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                        <  1:1419918/134‑1 (MQ=255)
                                                                                                       aGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGc                                                      <  1:565994/135‑1 (MQ=255)
                                                                                                        gACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTaaa                                                                                                           >  1:1145255/1‑81 (MQ=255)
                                                                                                          cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                     <  1:917840/149‑1 (MQ=255)
                                                                                                          cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                     >  2:174878/1‑149 (MQ=255)
                                                                                                          cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                     <  4:1403153/149‑1 (MQ=255)
                                                                                                           aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGa                                                                                 >  4:740871/1‑104 (MQ=255)
                                                                                                           aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGa                                                                                 <  3:740871/104‑1 (MQ=255)
                                                                                                           aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGaa                                                                        >  2:1265693/1‑113 (MQ=255)
                                                                                                           aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGaa                                                                        <  1:1265693/113‑1 (MQ=255)
                                                                                                           aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGat                                    >  2:446330/1‑149 (MQ=255)
                                                                                                                gcgTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTc                               <  3:984713/149‑1 (MQ=255)
                                                                                                                  gTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACaaa                                                                             <  1:1023440/101‑1 (MQ=255)
                                                                                                                  gTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACaaa                                                                             >  2:1023440/1‑101 (MQ=255)
                                                                                                                    cTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                             >  3:910224/1‑147 (MQ=255)
                                                                                                                    cTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                             <  4:910224/147‑1 (MQ=255)
                                                                                                                     tGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGtt                          <  2:966186/149‑1 (MQ=255)
                                                                                                                       aCCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGatta                                  >  3:1279664/1‑139 (MQ=255)
                                                                                                                        cccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCt                       <  1:363267/149‑1 (MQ=255)
                                                                                                                         ccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGt                                                                                                              >  3:515845/1‑61 (MQ=255)
                                                                                                                         ccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTg                      >  4:89080/1‑149 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCtt                                                <  2:786157/116‑1 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTg                                            >  3:1422249/1‑120 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTg                                            <  4:1422249/120‑1 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                     <  4:1233069/127‑1 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt               >  3:1043069/1‑149 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt               >  3:323002/1‑149 (MQ=255)
                                                                                                                                atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTCGCTGAAGAAGt               <  3:860872/149‑1 (MQ=255)
                                                                                                                                 tcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCg                                                                          >  2:304499/1‑89 (MQ=255)
                                                                                                                                 tcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCg                                                                          <  1:304499/89‑1 (MQ=255)
                                                                                                                                 tcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTc              >  4:1041478/1‑149 (MQ=255)
                                                                                                                                   atTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAAc                                                           >  4:1367318/1‑102 (MQ=255)
                                                                                                                                   atTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAAc                                                           <  3:1367318/102‑1 (MQ=255)
                                                                                                                                   atTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGaagaag                >  2:1350327/1‑145 (MQ=255)
                                                                                                                                   atTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCgg            >  4:958550/1‑149 (MQ=255)
                                                                                                                                    ttGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGaaa                                         <  2:1451774/119‑1 (MQ=255)
                                                                                                                                    ttGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGCAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGaaa                                         >  1:1451774/1‑119 (MQ=255)
                                                                                                                                      gCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt         <  4:1093159/149‑1 (MQ=255)
                                                                                                                                        tGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                     <  2:1425170/119‑1 (MQ=255)
                                                                                                                                         cttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                             >  3:469136/3‑126 (MQ=255)
                                                                                                                                         cttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcg      >  4:783625/3‑149 (MQ=255)
                                                                                                                                          ttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAg                            <  4:565192/125‑1 (MQ=255)
                                                                                                                                          ttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAg                            >  3:565192/2‑126 (MQ=255)
                                                                                                                                          ttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcgt     >  3:938680/2‑149 (MQ=255)
                                                                                                                                             ttCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATg  >  1:106243/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATG  >  NC_000913/4552036‑4552325

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