Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,505 G→T 12.5% E29* (GAG→TAG) ‡ ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,5050GT12.5% 332.6 / 30.9 104E29* (GAG→TAG) ‡ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (42/49);  new base T (5/8);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

Only 100 of 105 total aligned reads displayed.

AACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTC  >  NC_000913/3180357‑3180648
                                                                                                                                                    |                                                                                                                                                
aacaaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGGTACCGATGGAAAGAGTCAGTATTACCGATCGCCCGGACTCGCGTGAGAACGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCg                                                                                                                                                  >  3:562261/1‑149 (MQ=255)
   aaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAg                                                                                                                                                >  1:301025/1‑148 (MQ=255)
      gCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGt                                                                                                                                            >  4:449434/1‑149 (MQ=255)
        aCTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                          >  3:731054/1‑149 (MQ=255)
        aCTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                          <  4:359238/149‑1 (MQ=255)
        aCTATGCAGCGTAGTGCAACCGGTACCTCTTCTAGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCACACGAATACGGTTTCAATTTTCACACCATTTACGGCGAGCCGTa                                                                                                                                           >  1:543671/1‑148 (MQ=255)
          tATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatg                                                                                                                                       <  3:382694/149‑4 (MQ=255)
          tATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatg                                                                                                                                       <  4:297475/149‑4 (MQ=255)
            tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTgg                                                                                                                                      <  2:575638/148‑1 (MQ=255)
             gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCTCGCCGTAcatgg                                                                                                                                      <  2:193971/148‑5 (MQ=255)
                   tAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAg                                                                                                                               <  2:317551/148‑1 (MQ=255)
                    aGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGAt                                                                                                                             <  3:573697/149‑1 (MQ=255)
                      tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATg                                                                                                                            <  2:561630/148‑1 (MQ=255)
                      tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATg                                                                                                                            <  2:729862/148‑1 (MQ=255)
                      cgCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                           >  4:316250/2‑149 (MQ=255)
                         aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCtt                                                                                                                         <  2:235330/148‑1 (MQ=255)
                         aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTta                                                                                                                        >  3:566151/1‑149 (MQ=255)
                         aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTta                                                                                                                        <  4:808882/149‑1 (MQ=255)
                         aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGCGTGAAGATGCtt                                                                                                                         >  2:199780/1‑148 (MQ=255)
                          aCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtac                                                                                                                       <  3:695628/149‑1 (MQ=255)
                           ccGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtac                                                                                                                       >  1:220289/1‑148 (MQ=255)
                              gTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtactac                                                                                                                    >  1:372237/1‑148 (MQ=255)
                               tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACaa                                                                                                                  >  4:103848/1‑149 (MQ=255)
                                aCCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAg                                                                                                                 >  4:524754/1‑149 (MQ=255)
                                 ccTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGt                                                                                                                >  4:165857/1‑149 (MQ=255)
                                    cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTg                                                                                                              >  1:249625/1‑148 (MQ=255)
                                    cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTg                                                                                                              >  2:350708/1‑148 (MQ=255)
                                    cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGa                                                                                                             >  3:811784/1‑149 (MQ=255)
                                     ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGa                                                                                                             >  2:18772/1‑148 (MQ=255)
                                     ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAc                                                                                                            >  4:439399/1‑149 (MQ=255)
                                     ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAc                                                                                                            >  4:732138/1‑149 (MQ=255)
                                     ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAc                                                                                                            >  3:204755/1‑149 (MQ=255)
                                               aTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGGTGTGaaaa                                                                                                                              >  4:565221/1‑119 (MQ=38)
                                               aTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGGTGTGaaaa                                                                                                                              <  3:565221/121‑3 (MQ=38)
                                                   gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTg                                                                                              >  3:468392/1‑149 (MQ=255)
                                                    tGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGa                                                                                             >  4:448499/1‑149 (MQ=255)
                                                         gATACTGATGGAAAGAGTCAGTATTACCGAGAGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatg                                                                                                                                       <  2:728288/102‑4 (MQ=37)
                                                         gATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTg                                                                                                              <  2:915250/127‑1 (MQ=255)
                                                         gATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTg                                                                                                              >  1:915250/1‑127 (MQ=255)
                                                         gATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAg                                                                                         >  1:735228/1‑148 (MQ=255)
                                                          aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                       <  3:103848/149‑1 (MQ=255)
                                                           tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                       <  2:785613/148‑1 (MQ=255)
                                                                 tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                <  3:85032/149‑1 (MQ=255)
                                                                 tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                <  3:449434/149‑1 (MQ=255)
                                                                 tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                <  3:316250/149‑1 (MQ=255)
                                                                 tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaaga                                                                                 <  1:886536/148‑1 (MQ=255)
                                                                  ggAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGaagaag                                                                               >  4:533137/1‑149 (MQ=255)
                                                                             gTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGa                                                                     <  2:372237/148‑1 (MQ=255)
                                                                              tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGc                                                                                        <  3:611724/128‑1 (MQ=255)
                                                                              tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGc                                                                                        >  4:611724/1‑128 (MQ=255)
                                                                              tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGaa                                                                    <  2:743943/148‑1 (MQ=255)
                                                                                 tACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                 <  1:18772/148‑1 (MQ=255)
                                                                                     gAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGCCACCGCCGAACTGCAcc                                                             <  2:301025/148‑1 (MQ=255)
                                                                                         gcCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                         <  1:199780/148‑1 (MQ=255)
                                                                                          cccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                        <  2:882134/148‑1 (MQ=255)
                                                                                            cGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgtg                                                      >  1:655964/1‑148 (MQ=255)
                                                                                              gACTGGCATGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGcc                                                    <  2:87117/148‑1 (MQ=255)
                                                                                                cTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGa                                                 >  4:352079/1‑149 (MQ=255)
                                                                                                 tGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaa                                                >  3:648460/1‑149 (MQ=255)
                                                                                                  ggCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaaa                                               <  3:352079/149‑1 (MQ=255)
                                                                                                     gTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                            <  4:45528/149‑1 (MQ=255)
                                                                                                      tGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGc                                                                >  2:564106/1‑128 (MQ=255)
                                                                                                      tGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGc                                                                <  1:564106/128‑1 (MQ=255)
                                                                                                       gagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                          <  3:732138/149‑1 (MQ=255)
                                                                                                       gagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                           >  2:658987/1‑148 (MQ=255)
                                                                                                          aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGtggtgg                                        >  1:882791/1‑148 (MQ=255)
                                                                                                              cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaa                                    <  1:546408/148‑1 (MQ=255)
                                                                                                              cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaa                                    <  1:675613/148‑1 (MQ=255)
                                                                                                               ccACGAATACGGTTTCAATTTTCAAACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAATAg                                  <  4:566151/149‑1 (MQ=255)
                                                                                                                cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGt                                                                                                >  2:306805/1‑86 (MQ=255)
                                                                                                                cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGt                                                                                                <  1:306805/86‑1 (MQ=255)
                                                                                                                cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                  >  2:123044/1‑148 (MQ=255)
                                                                                                                cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                 <  4:811784/149‑1 (MQ=255)
                                                                                                                cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                 >  4:134846/1‑149 (MQ=255)
                                                                                                                  cGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAGCTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                >  1:296473/1‑148 (MQ=255)
                                                                                                                   gAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGa                               <  2:220289/148‑1 (MQ=255)
                                                                                                                      tACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                              <  1:621747/146‑1 (MQ=255)
                                                                                                                      tACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                              >  2:621747/1‑146 (MQ=255)
                                                                                                                      tACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                             <  4:638163/147‑1 (MQ=255)
                                                                                                                      tACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTAACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGCGGAAAAAGTGATc                             >  3:638163/1‑147 (MQ=255)
                                                                                                                       aCGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                   <  3:532066/140‑1 (MQ=255)
                                                                                                                       aCGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                   >  4:532066/1‑140 (MQ=255)
                                                                                                                           tttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCg                      <  3:448499/149‑1 (MQ=255)
                                                                                                                            ttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                           <  1:601303/127‑1 (MQ=255)
                                                                                                                            ttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                           >  2:601303/1‑127 (MQ=255)
                                                                                                                              cAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                    <  2:249625/148‑1 (MQ=255)
                                                                                                                                    tCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga             <  4:731054/149‑1 (MQ=255)
                                                                                                                                    tCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga             <  3:524754/149‑1 (MQ=255)
                                                                                                                                    tCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg              <  1:658987/148‑1 (MQ=255)
                                                                                                                                    tCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg              <  1:758384/148‑1 (MQ=255)
                                                                                                                                    tCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg              <  1:619878/148‑1 (MQ=255)
                                                                                                                                    tCACACCAGTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGATCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg              >  1:587533/1‑148 (MQ=255)
                                                                                                                                       cacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAGCTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatg           <  2:296473/148‑1 (MQ=255)
                                                                                                                                       cacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCACCGATGAGctgctg           >  1:230093/1‑148 (MQ=255)
                                                                                                                                        acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc         >  4:226578/1‑149 (MQ=255)
                                                                                                                                            agTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                           >  1:858652/2‑112 (MQ=255)
                                                                                                                                            agTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                           <  2:858652/111‑1 (MQ=255)
                                                                                                                                             gTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaa      >  1:556878/1‑148 (MQ=255)
                                                                                                                                             gTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCaa      <  2:556878/148‑1 (MQ=255)
                                                                                                                                                cGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTc  >  2:184830/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTC  >  NC_000913/3180357‑3180648

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: