Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,425,696 C→A N59K (AAC→AAA priB → primosomal replication protein N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,425,6960CA94.2% 146.0 / ‑3.2 52N59K (AAC→AAApriBprimosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (17/32);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACT  >  NC_000913/4425451‑4425932
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
cGTCGAGATGATTTCGCACACGAAACCACTGATGATGCTGCAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTATCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCAATCAGTAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCAGTTATTGTTATCGGAAACGCAAAACAGGc                                                                                                                                                                                                                                         >  2:307826/1‑251 (MQ=255)
           tttCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTcac                                                                                                                                                                                                                               >  1:286903/1‑250 (MQ=255)
              cGCAAACGAAACCGCGGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCCCCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGta                                                                                                                                                                                                                           <  2:223291/251‑1 (MQ=255)
                 aaaCGAAACCGTTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCTGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATaa                                                                                                                                                                                                                        <  1:397751/251‑1 (MQ=255)
                 aaaCGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCAGGCACCGTGTGCAGTGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGCAAAACAGGCCATTACTCACAGtata                                                                                                                                                                                                                         >  1:154540/1‑250 (MQ=255)
                   aCGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACg                                                                                                                                                                                                                      <  2:384210/251‑1 (MQ=255)
                   aCGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCACCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCCCAGTATAACg                                                                                                                                                                                                                      >  2:272943/1‑251 (MQ=255)
                               atgatgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGc                                                                                                                                                                                                          >  1:200290/1‑251 (MQ=255)
                                             ggggATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCTCCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAgg                                                                                                                                                                                             <  2:58761/250‑1 (MQ=255)
                                                gatgcgaaaGAGTAATTTCAGATGACCAACCGGCTGGTTTTGTCCGGCACCGTGTGCAGGGCCCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACGGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTTTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGtt                                                                                                                                                                                         <  2:177005/244‑1 (MQ=255)
                                                           gTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATGGTTAGCGGAAACGAAAAACAGGCCAGTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGcc                                                                                                                                                                              <  1:223652/251‑1 (MQ=255)
                                                             aaTTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAAc                                                                                                                                                                                                     >  1:310262/1‑226 (MQ=255)
                                                             aaTTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAAc                                                                                                                                                                                                     <  2:310262/226‑1 (MQ=255)
                                                                 tCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCAt                                                                                                                                                                                                        >  1:431381/1‑219 (MQ=255)
                                                                    gatgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCaa                                                                                                                                                                     <  1:329236/251‑1 (MQ=255)
                                                                       gtccAACCGTCGGGTGTTCTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTGCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGa                                                                                                                                                                  <  2:364044/249‑1 (MQ=255)
                                                                        aCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGaa                                                                                                                                                                 <  1:406814/251‑1 (MQ=255)
                                                                              cGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACCAGGCAAAGAACGGACt                                                                                                                                                           >  1:357429/1‑251 (MQ=255)
                                                                              cGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCACCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACt                                                                                                                                                           <  1:272943/251‑1 (MQ=255)
                                                                                         gtCCGGCCCCGTGTGCAGGGCTTCCCTTCGAAAGTTCAGTCCATCAGGAATCCCTCACGGCCAGTTCGTGCTTGGGCATCGTTCTGTGCAGGAGGAGGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATgg                                                                                                                                                <  2:170918/251‑1 (MQ=255)
                                                                                           ccGGCACCGTTTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCAGCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCTCATAACCGGTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtt                                                                                                                                              <  1:442645/251‑1 (MQ=255)
                                                                                             ggccccGTGTTCAGGGCTCCCCTTCGAAAGGTCAGTCCAGCAGGAATTCCTCACTGCCATTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAAGGCCCGTTATTTTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtttt                                                                                                                                            <  1:39948/247‑1 (MQ=255)
                                                                                                aCCGTGTGCAGGGCTCCCCTTCTAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTg                                                                                                                                           <  2:255952/249‑1 (MQ=255)
                                                                                                aCCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTg                                                                                                                                           >  1:255952/1‑249 (MQ=255)
                                                                                                          gggCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGc                                                                                                                                 >  1:329460/1‑249 (MQ=255)
                                                                                                              cccccTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGattg                                                                                                                           <  2:445014/250‑1 (MQ=255)
                                                                                                                cccTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATAGTTCTGTGCAGGAGGAAGCCGGATTTCACCGGCAGGCGTGGTGTCAAATGCCAGTTATTGATAGCGGAAACGAAAAACAGGCCAATACTCACAGCATAAAGGTCGGCAGTCGCATAACCGTTCAGGGGCACATTTCATGCAACAAGGAAAAGAAAGGAATAAGCAAAATGGtgatgaatgccagggagattga                                                                                                                          >  1:426891/1‑229 (MQ=255)
                                                                                                                       aaatgtccgTCCATCAGGAATTCCTCCCTGCCAGTTCTTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCCCCGGCAGGCGTTGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCTAGCAGATTGAATTGATat                                                                                                                  <  1:333358/243‑2 (MQ=255)
                                                                                                                       aaaGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAg                                                                                                                  <  1:333355/251‑1 (MQ=255)
                                                                                                                             cAGTCCATCAGGAATTCCTCACTGCCAGTTCGTTCTTGAGCAGCGTGCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACATGATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGGTTTGCATGCCGAGCAGATTGAATTGATAGATTCTg                                                                                                            <  2:275028/251‑1 (MQ=255)
                                                                                                                               gTCCAGCAGGAATTCCTCCCTGCCCTTTCGTGCGTGAGCATCGTGCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTTAAATGCCCGTTATTGTTATCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTGGCAGGCCGTGCAGATTGAATCGATAGATTCTGga                                                                                                          <  2:300459/251‑1 (MQ=255)
                                                                                                                                cccaacAGGAATTCCCCCCTGCCAGTGCGTGCTTTAGCATCGTTCGGTGCAGGAGGAATCCGGCTTTCACCTGCAGGCGTGGTTTCAAATCCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGgag                                                                                                         <  2:455330/246‑1 (MQ=255)
                                                                                                                                                cattGCCAGTTCGTGTTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTAAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGc                                                                                           <  2:216225/246‑1 (MQ=255)
                                                                                                                                                                               tgCAGGAGGAAGCCGGCTTTCACCGGCAGGAGTGGTGTCAAATGCCAGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAAACGTTCAGGGGTTCATTTAATGCCACAAGGCAAAGAACGGAATGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTTCGTCGTTGAAAGTTTTGGCGttt                                                          >  2:213681/1‑251 (MQ=255)
                                                                                                                                                                                 cAGGAGGAAGCCTGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGGTCTGCCGTTTCa                                                        <  1:348912/251‑1 (MQ=255)
                                                                                                                                                                                 aaggaggAAGCCGGCTTTCACCGGCAGGCGTGGTGGCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCGTTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                        <  1:263098/250‑1 (MQ=255)
                                                                                                                                                                                  aggaggAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCAc                                                       >  1:95706/1‑251 (MQ=255)
                                                                                                                                                                                            cGGCTTTCCCCGGCGGGCGTGGTGTCAAATGCCCGTTATTGTTTGCGGAAACGAAAAACAGGCCATTACTCACCGTATCCCGGTCGGCAGACTCATAACCGTTCACGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCTTGCCGAGCTGATTGAATTGGTAGATTATGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCACAGGTTCTGCCGTTTCACTGCGGCAGGc                                             >  2:207643/1‑251 (MQ=255)
                                                                                                                                                                                                 ttCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                         <  1:213681/250‑1 (MQ=255)
                                                                                                                                                                                                  tCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCAGTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGTAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCa                                        <  1:328658/250‑1 (MQ=255)
                                                                                                                                                                                                   ccccGGCGGGCGTGGGGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACCCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATGGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTGCaa                                       <  2:159568/248‑1 (MQ=255)
                                                                                                                                                                                                       ggCAGGCGGGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTACAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                  <  2:199143/251‑1 (MQ=255)
                                                                                                                                                                                                        gCAGGCGTGGTGGCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCGCAGTATAACGGTCGGCAGTCGGATATCCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGGCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                 <  1:307826/251‑1 (MQ=255)
                                                                                                                                                                                                           ggCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTCTAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTCCATGCCGAGCAGATTGAATTGATAGCTTATGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGAGTTCAAGAGATCGAc                              >  1:302478/1‑251 (MQ=255)
                                                                                                                                                                                                                 gtgtCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGCTCAGGGGTTCATTTCCTGCCACAAGGCAAAGAACGGACTGAGCACAATGGTTTTGCATGCCGAGAAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTATGCCGCTTCACCGCGGCAGGCGATCAAGAGATCGACTCTaa                         >  1:445737/1‑250 (MQ=255)
                                                                                                                                                                                                                        aTGCCCGTTATTGTGAGCGGAAACGAAAAACAGGCCATTACTCAAAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCg                 <  1:320563/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTAc             <  1:461185/251‑1 (MQ=255)
                                                                                                                                                                                                                                tATTGTTAGCGGGAACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGAAAAGAACGGAATGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATCGGATCTGGCGACTAGCCATATGGCAAGTTATTTCCGTCGTCGCACGTTCGGCCGTTACCCCGGGGAAGGCGTTCAAGAGATCGCCTATAAAGatatagcagagcgg         >  2:204919/1‑241 (MQ=255)
                                                                                                                                                                                                                                tATTGTTAGCGGAAACGAAAAACAGGCCATTTCTCACAGTATAACGGCCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGc                <  1:204919/244‑1 (MQ=255)
                                                                                                                                                                                                                                tATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCt          <  1:113610/250‑1 (MQ=255)
                                                                                                                                                                                                                                tATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTg         >  1:4767/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACAGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACt  >  1:414488/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
CGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACT  >  NC_000913/4425451‑4425932

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