Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,353,758 T→A D254V (GAT→GTT)  sapC ← putrescine ABC exporter membrane protein SapC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,353,7580TA97.7% 143.9 / ‑3.3 43D254V (GAT→GTT) sapCputrescine ABC exporter membrane protein SapC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (22/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

TTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATA  >  NC_000913/1353533‑1354003
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
ttCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAgg                                                                                                                                                                                                                              >  1:303968/1‑251 (MQ=255)
 tCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGGCATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAgg                                                                                                                                                                                                                              <  1:46886/250‑1 (MQ=255)
         ggTCATACTTACCCGGTCGACGGCTTTAACCCACTCCTCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAgggg                                                                                                                                                                                                                       >  2:317044/1‑249 (MQ=255)
         ggTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAg                                                                                                                                                                                                                     <  2:190191/251‑1 (MQ=255)
               aCTTACCCGTTCGACGGCTTTAACCCACTCATCAACGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGGCATCGCCAAGGAGGTGAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGTAGTtgtg                                                                                                                                                                                                                <  2:40515/250‑1 (MQ=255)
                cTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTgcg                                                                                                                                                                                                              >  2:162709/1‑251 (MQ=255)
                   aCCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGcc                                                                                                                                                                                                            <  1:321981/250‑1 (MQ=255)
                     ccGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCg                                                                                                                                                                                                           >  1:149925/1‑249 (MQ=255)
                                  ttAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGc                                                                                                                                                                                              <  1:321207/249‑1 (MQ=255)
                                       ccACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAgcgc                                                                                                                                                                                       <  2:95672/251‑1 (MQ=255)
                                              tcaCCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTACCCAACACAACGCAAACCATAATTCCCGCACCGGCCAGCATGACGTTCCACGGGGCGACATAAATAGGTTCCAGCGCAACACGAAGAATGCCTCCCCATCAAGGGCAGGGGCGTGTTGGGCCGCGACTGAGAAAGCCCAGCGCGGCGaac                                                                                                                                                                                >  2:104246/1‑249 (MQ=255)
                                                 ccGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCTGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCg                                                                                                                                                                             >  2:258123/1‑251 (MQ=255)
                                                                        gTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGcgcg                                                                                                                                                      >  2:283283/1‑251 (MQ=255)
                                                                            cGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTg                                                                                                                                                  <  1:233120/251‑1 (MQ=255)
                                                                             gAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATc                                                                                                                                                                                                                                            >  1:245363/1‑160 (MQ=255)
                                                                             gAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATc                                                                                                                                                                                                                                            <  2:245363/160‑1 (MQ=255)
                                                                                tatCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGAtct                                                                                                                                              >  1:222176/1‑251 (MQ=255)
                                                                                 atCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGAtctc                                                                                                                                             <  2:22553/251‑1 (MQ=255)
                                                                                       tAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGCGATCTCGGTGAc                                                                                                                                       >  1:332249/1‑251 (MQ=255)
                                                                                             cATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGACGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGAGCGGCGATATCGAGTATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGccc                                                                                                                                 >  1:197819/1‑251 (MQ=255)
                                                                                                      aCCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGAt                                                                                                                        <  2:198573/251‑1 (MQ=255)
                                                                                                              aataatCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGc                                                                                                                 >  2:40156/1‑250 (MQ=255)
                                                                                                                          acgaACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAAcc                                                                                                    >  2:103527/1‑251 (MQ=255)
                                                                                                                                 ccATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAATCCAGAGAat                                                                                             <  1:86899/251‑1 (MQ=255)
                                                                                                                                    tCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATAtt                                                                                          >  1:84375/1‑251 (MQ=255)
                                                                                                                                        cAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATCTTCaac                                                                                      >  2:189040/1‑249 (MQ=255)
                                                                                                                                          aggaggTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGt                                                                                    >  1:346256/1‑251 (MQ=255)
                                                                                                                                             aggTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTg                                                                                  >  1:68042/1‑250 (MQ=255)
                                                                                                                                              ggTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGAt                                                                                <  1:8636/251‑1 (MQ=255)
                                                                                                                                                ttAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGCTGTTTTGCATCACCGCAAACCAGAGAATATTCCGCGTTGATgg                                                                              >  2:161679/1‑250 (MQ=255)
                                                                                                                                                   aCCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGcc                                                                           >  1:161539/1‑251 (MQ=255)
                                                                                                                                                                tAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcg                                                              >  2:38642/1‑251 (MQ=255)
                                                                                                                                                                 aaTCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGggcg                                                              <  1:49651/250‑1 (MQ=255)
                                                                                                                                                                     aTAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGAt                                                         <  1:305037/251‑1 (MQ=255)
                                                                                                                                                                               cacgtGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATAGCGAGAAGTGCCATCGACATTGCGCGGGGGATCTCGTTGAACAGCCCGTCGGTGATGTTGGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTCctct                                               <  1:217119/247‑1 (MQ=255)
                                                                                                                                                                                 ccTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGAGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCGGAGAATATTCAGCGTTGATGCGCCATCCAGACCGGCGGCGATAAAGTACTCttt                                             >  2:208991/1‑251 (MQ=255)
                                                                                                                                                                                   tGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCTATCCAGACGGGCGGCGATAACGTACTCttttt                                           >  2:110190/1‑251 (MQ=255)
                                                                                                                                                                                      cagcaTGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCa                                        <  1:229969/251‑1 (MQ=255)
                                                                                                                                                                                         caTGACAGTCCACTGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTACCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGGGAGGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCTATAACGTACTCTTTTTCCAGtt                                     <  1:342700/251‑1 (MQ=255)
                                                                                                                                                                                         caTGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTACCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGtt                                     <  1:342702/251‑1 (MQ=255)
                                                                                                                                                                                                     cGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCa                         <  2:103725/251‑1 (MQ=255)
                                                                                                                                                                                                         gCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTGAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCt                     <  2:327416/251‑1 (MQ=255)
                                                                                                                                                                                                                            gcgcAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTGGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATa  <  2:347538/251‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
TTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCATA  >  NC_000913/1353533‑1354003

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