Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,920,305 C→A 77.8% A33A (GCG→GCT atpH ← F1 sector of membrane‑bound ATP synthase, delta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,920,3050CA77.8% 267.2 / 89.9 126A33A (GCG→GCTatpHF1 sector of membrane‑bound ATP synthase, delta subunit
Reads supporting (aligned to +/- strand):  ref base C (11/17);  new base A (36/62);  total (47/79)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 128 total aligned reads displayed.

GCCATAACCCGAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGTTTATCCACGATGTC  >  NC_000913/3920162‑3920450
                                                                                                                                               |                                                                                                                                                 
gCCATAACCCGAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAg                                                                                                                                              <  4:526103/149‑1 (MQ=255)
  cATAACCCGAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCa                                                                                                                                            >  2:636701/1‑149 (MQ=255)
  cATAACCCGAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCa                                                                                                                                            >  3:57214/1‑149 (MQ=255)
        ccGAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCt                                                                                                                                      <  3:397057/149‑1 (MQ=255)
          gAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGc                                                                                                                                    <  3:886036/149‑1 (MQ=255)
             tCAGGTTCTGACCGTTTTCGTCCAGTTGGTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAg                                                                                                                                 >  3:564862/1‑149 (MQ=255)
             tCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGACCTGCCAg                                                                                                                                 >  1:525618/1‑149 (MQ=255)
               aGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCg                                                                                                                               >  2:586743/1‑149 (MQ=255)
                 gTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGtt                                                                                                                             <  4:554397/149‑1 (MQ=255)
                  ttCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTc                                                                                                                            >  3:649280/1‑149 (MQ=255)
                   tCTGACCGTTTTCGTCCAGTTGCGCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCt                                                                                                                           <  1:795677/149‑1 (MQ=255)
                       aCCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTTGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTacac                                                                                                                       >  2:569175/1‑149 (MQ=255)
                       aCCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTacac                                                                                                                       >  1:76393/1‑149 (MQ=255)
                         cGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGc                                                                                                                                             >  3:609063/1‑125 (MQ=255)
                         cGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACtt                                                                                                                     <  1:343852/149‑1 (MQ=255)
                         cGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACtt                                                                                                                     <  4:314030/149‑1 (MQ=255)
                         cGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACtt                                                                                                                     <  3:338408/149‑1 (MQ=255)
                         cGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACtt                                                                                                                     <  1:742472/149‑1 (MQ=255)
                          gTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACttt                                                                                                                    >  3:249492/1‑149 (MQ=255)
                          gTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACttt                                                                                                                    <  1:416053/149‑1 (MQ=255)
                           ttttCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTg                                                                                                                   <  4:619335/149‑1 (MQ=255)
                           ttttCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACcttt                                                                                                                   <  4:254322/149‑2 (MQ=255)
                           ttttCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTg                                                                                                                   <  4:663915/149‑1 (MQ=255)
                           ttttCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTg                                                                                                                   <  3:107050/149‑1 (MQ=255)
                             ttCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGgt                                                                                                                 <  4:289395/149‑1 (MQ=255)
                             ttCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGgt                                                                                                                 <  4:649280/149‑1 (MQ=255)
                               cGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGgtgt                                                                                                               <  1:707796/149‑1 (MQ=255)
                                gTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGtt                                                                                                              <  1:249575/149‑1 (MQ=255)
                                 tCCAGTTGCTCACCACAAACTGCGATAAATGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTc                                                                                                             >  1:642558/1‑149 (MQ=255)
                                 tCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTc                                                                                                             >  2:31979/1‑149 (MQ=255)
                                 tCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTc                                                                                                             <  3:45972/149‑1 (MQ=255)
                                 tCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTc                                                                                                             <  1:569175/149‑1 (MQ=255)
                                   cAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGa                                                                                                           <  2:15131/149‑1 (MQ=255)
                                   cAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGa                                                                                                           <  4:165005/149‑1 (MQ=255)
                                    aGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGGGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGAc                                                                                                          >  2:192040/1‑149 (MQ=255)
                                     gTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGAc                                                                                                          <  4:328199/148‑1 (MQ=255)
                                     gTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACg                                                                                                         >  1:112759/1‑149 (MQ=255)
                                              cacaAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCa                                                                                                >  4:719240/1‑149 (MQ=255)
                                              cacaAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCa                                                                                                >  2:124666/1‑149 (MQ=255)
                                              cacaAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCa                                                                                                >  3:427254/1‑149 (MQ=255)
                                              cacaAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCa                                                                                                >  2:221128/1‑149 (MQ=255)
                                                  aaCTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAg                                                                                            <  3:357096/149‑1 (MQ=255)
                                                  aaCTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAg                                                                                            <  2:751435/149‑1 (MQ=255)
                                                    cTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGgt                                                                                                                 >  3:597559/1‑126 (MQ=255)
                                                    cTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGgt                                                                                                                 <  4:597559/126‑1 (MQ=255)
                                                        gATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCtt                                                                                      <  2:76393/149‑1 (MQ=255)
                                                        gATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCtt                                                                                      <  4:511331/149‑1 (MQ=255)
                                                         aTAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCttt                                                                                     <  4:249492/149‑1 (MQ=255)
                                                         aTAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCttt                                                                                     <  3:828051/149‑1 (MQ=255)
                                                            aaCGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGc                                                                                  <  2:723002/149‑1 (MQ=255)
                                                            aaCGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGc                                                                                  <  4:68890/149‑1 (MQ=255)
                                                             atgACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCg                                                                                 <  2:642558/147‑1 (MQ=255)
                                                               gACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTa                                                                               <  4:205771/149‑1 (MQ=255)
                                                                 cTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAgg                                                                             <  1:194411/149‑1 (MQ=255)
                                                                  tCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAg                                                                                                   <  3:673948/126‑1 (MQ=255)
                                                                  tCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAgctg                                                                                         <  4:724296/136‑1 (MQ=255)
                                                                    ggCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGc                                                                          <  3:719240/149‑1 (MQ=255)
                                                                     gCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCg                                                                         >  3:251838/1‑149 (MQ=255)
                                                                        aGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGc                                                                      >  2:605771/1‑149 (MQ=255)
                                                                           gTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTAc                                                                   <  1:221128/149‑1 (MQ=255)
                                                                            tttCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTAcc                                                                  <  4:325810/149‑1 (MQ=255)
                                                                             ttCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCg                                                                 <  3:53421/149‑1 (MQ=255)
                                                                                  gcgcAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGtaata                                                            >  2:367660/1‑149 (MQ=255)
                                                                                       aGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTc                                                       <  1:636701/149‑1 (MQ=255)
                                                                                          gcgcCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGa                                                    >  3:176413/1‑149 (MQ=255)
                                                                                             ccAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACAt                                                 >  1:180337/1‑149 (MQ=255)
                                                                                              cAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATc                                                >  4:171544/1‑149 (MQ=255)
                                                                                              cAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATc                                                <  4:251838/149‑1 (MQ=255)
                                                                                               agagagAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCa                                               >  3:795937/1‑149 (MQ=255)
                                                                                               agagagAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCa                                               >  4:193407/1‑149 (MQ=255)
                                                                                                  gagaAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGa                                                    <  2:641335/141‑1 (MQ=255)
                                                                                                  gagaAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGa                                                    >  1:641335/1‑141 (MQ=255)
                                                                                                   agaAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGccc                                           <  4:57214/149‑1 (MQ=255)
                                                                                                   agaAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTCTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGccc                                           >  4:556304/1‑149 (MQ=255)
                                                                                                        ctctGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAg                                                                                            <  4:797692/95‑1 (MQ=255)
                                                                                                        ctctGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAg                                                                                            >  3:797692/1‑95 (MQ=255)
                                                                                                         tctGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCccctccct                                     <  2:447134/149‑1 (MQ=255)
                                                                                                          ctGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCcctccctc                                    <  4:879996/149‑1 (MQ=255)
                                                                                                          ctGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCcctccctc                                    <  1:124666/149‑1 (MQ=255)
                                                                                                             ccATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTAcc                                                                  <  2:238690/116‑1 (MQ=255)
                                                                                                             ccATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTAcc                                                                  >  1:238690/1‑116 (MQ=255)
                                                                                                             ccATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCtt                                 <  3:151455/149‑1 (MQ=255)
                                                                                                              cATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAggg                                                                            <  4:863333/105‑1 (MQ=255)
                                                                                                              cATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAggg                                                                            >  3:863333/1‑105 (MQ=255)
                                                                                                              cATTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTa                                <  3:256667/149‑1 (MQ=255)
                                                                                                               aTTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTc                                                       <  1:563512/125‑1 (MQ=255)
                                                                                                               aTTTGTTCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTAc                               >  4:712623/1‑149 (MQ=255)
                                                                                                                     tCGTTTTTGGTTACCTCGGCGGCAAACGTCAGCATGTCCTGCCAGCGTTCTACACATTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGCCATCAGCCCCTCCCGCCTTACAGTTCa                         <  2:175762/149‑1 (MQ=255)
                                                                                                                     tCGTTTTTGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCa                         <  3:712623/149‑1 (MQ=255)
                                                                                                                          tttGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGAc                    >  1:583035/1‑149 (MQ=255)
                                                                                                                          tttGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGAc                    <  1:367660/149‑1 (MQ=255)
                                                                                                                            tGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACaa                  <  3:725571/149‑1 (MQ=255)
                                                                                                                            tGGTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACaa                  <  4:339497/149‑1 (MQ=255)
                                                                                                                              gTTACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGt                >  1:249883/1‑149 (MQ=255)
                                                                                                                                tACCTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGttt              <  2:481542/149‑1 (MQ=255)
                                                                                                                                  ccTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGTTTAt            <  2:97267/149‑1 (MQ=255)
                                                                                                                                   cTCGGCGGCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGTTTATc           <  1:691704/149‑1 (MQ=255)
                                                                                                                                     cggcggCAAAAGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGTTTATCCa         >  3:492552/1‑149 (MQ=255)
                                                                                                                                            aaaaGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGTTTATCCACGATGTc  <  4:496142/149‑1 (MQ=255)
                                                                                                                                              aCGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAg                                                     >  2:650467/1‑96 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCCATAACCCGAATCAGGTTCTGACCGTTTTCGTCCAGTTGCTCACCACAAACTGCGATAAACGACTCGGCGAGCGTTTCTGGCGCAAGCGCGCCAGAGAGAAGCTCTGCCATTTGTTCGTTTTTGGTTACCTCGGCGGCAAACGCCAGCATGTCCTGCCAGCGTTCTACACTTTGGTGTTCGACGGCAAAGTCAAAAGCTGCTTTGGCGTAGGGGCGAGCTACCGTAATAAATTCAGACATCAGCCCCTCCCTCCTTACAGTTCAGCGACAAGTTTATCCACGATGTC  >  NC_000913/3920162‑3920450

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