Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,473,882 T→G 13.1% L94R (CTG→CGG)  yfdP → CPS‑53 (KpLE1) prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,473,8820TG13.1% 152.0 / 3.8 69L94R (CTG→CGG) yfdPCPS‑53 (KpLE1) prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (21/36);  new base G (9/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 8.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCACTCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACTTA  >  NC_000913/2473736‑2474027
                                                                                                                                                  |                                                                                                                                                 
ctcactcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGa                                                                                                                                                 <  2:621506/149‑1 (MQ=255)
 tcactcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGAc                                                                                                                                                >  2:266546/1‑149 (MQ=255)
 tcactcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGAGGTATGTCAATCGCGGTGGTCGCCCGAc                                                                                                                                                >  2:288285/1‑149 (MQ=255)
     tcaACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAGATTTCCGAAGAAGACACCACCGAAACGCCGATGTGTGTCAATCGCGGTGGTCGCCTGACTAtt                                                                                                                                            <  1:949390/149‑1 (MQ=255)
      caACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTc                                                                                                                                           <  1:564807/149‑1 (MQ=255)
          cTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCa                                                                                                                                       >  1:1047046/1‑149 (MQ=255)
          cTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCa                                                                                                                                       <  1:659577/149‑1 (MQ=255)
          cTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCa                                                                                                                                       >  2:232845/1‑149 (MQ=255)
          cTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCa                                                                                                                                       >  2:279821/1‑149 (MQ=255)
           tGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGCGTACACCTAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAg                                                                                                                                      >  1:1363230/1‑149 (MQ=255)
              aaaTCAAACAATTAGATGGGAACCGTCTGAAAATCCAACACAACCCCTCATATGTTCGCCCGACTCATCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGaa                                                                                                                                   <  2:844553/149‑1 (MQ=255)
                   gAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCa                                                                                                                              <  1:1318371/149‑1 (MQ=255)
                    aaCCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACTATTCTGCAGGAAGaaccg                                                                                                                             >  2:15716/1‑145 (MQ=255)
                       cATTTGTTGACAACCGTCGGAAAATCAAACACGATCCTTCATATGGTCGCCCGACTAACCGTGGGAATATTTCCGAAGAAGACACCACCGAATCGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGa                                                                                                                          <  2:1199095/149‑1 (MQ=255)
                       cATTGGATGACAACCTGCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGa                                                                                                                          <  2:1199080/149‑1 (MQ=255)
                        aTTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTAGTCGCCTTACTATTCTGCTGGAAGACCAGGGAc                                                                                                                         >  1:886221/1‑149 (MQ=255)
                        aTTTGACGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTCCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGAc                                                                                                                         <  1:847335/149‑1 (MQ=255)
                          ttGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACaa                                                                                                                       <  1:321136/149‑1 (MQ=255)
                             atgaGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTa                                                                                                                    <  1:179323/149‑1 (MQ=255)
                                agaACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTg                                                                                                                 >  1:1139928/1‑149 (MQ=255)
                                agaACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGTACCATTCTGCAGGAAAACCAGGGCCAATTACTg                                                                                                                 >  2:266624/1‑149 (MQ=255)
                                 gaACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAACCGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGa                                                                                                                <  1:665992/149‑1 (MQ=255)
                                          gAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGGACCATTCAGAAGAAAGACCAGGGACAATTACTGACTCTGGCCg                                                                                                       >  2:624449/1‑149 (MQ=255)
                                          gAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCg                                                                                                       >  1:232143/1‑149 (MQ=255)
                                                 aaaCACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTCCCGGTGAGcc                                                                                                >  1:589433/1‑149 (MQ=255)
                                                  aaCACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCTGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCt                                                                                               >  1:663798/1‑149 (MQ=255)
                                                    cacaAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGa                                                                                             <  1:266546/149‑1 (MQ=255)
                                                     acaAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACTATTACTGACTCTTGCCGGTGAGCCTGAc                                                                                            <  1:874568/149‑1 (MQ=255)
                                                         acccccATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGATACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGaa                                                                                        <  2:886221/144‑1 (MQ=255)
                                                               atatGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAGGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTcc                                                                                  >  1:1300412/1‑149 (MQ=255)
                                                               atatGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTcc                                                                                  <  1:1284635/149‑1 (MQ=255)
                                                               atatGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTcc                                                                                  <  1:881256/149‑1 (MQ=255)
                                                                  tGTTCGCCCGACGAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCcgcg                                                                               <  2:395125/149‑1 (MQ=255)
                                                                    ttCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCACTCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGCCAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCgca                                                                             >  1:122584/1‑149 (MQ=255)
                                                                     tCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGGAAACTCCGCgcag                                                                            >  1:1127645/1‑149 (MQ=255)
                                                                        ctcGACTAACCGTGGGAAATTTTCCGAAGAAGCCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGCCAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGcagcag                                                                         <  2:1096788/147‑1 (MQ=255)
                                                                        cccGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGcagcag                                                                         >  1:957631/1‑149 (MQ=255)
                                                                         ccGACCAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAgg                                                                        <  1:1401583/149‑1 (MQ=255)
                                                                           gACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTc                                                                      <  1:262386/149‑1 (MQ=255)
                                                                            aCTAACCTTGGCAAAATTTCCGAAGAAGACCCCACCGAACCGCCGATGTATGCCAAGCGCGGTGGTCGCCTGACTATTCTCCAGCAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCa                                                                     <  1:406963/149‑1 (MQ=255)
                                                                                aCCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTTCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTaa                                                                 >  1:1124706/1‑149 (MQ=255)
                                                                                 ccGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGTATCATTATAAAAAAAGACCACGGACAATTACTGCATCTTTCCGGGGAGCCTGACGGAAAACTCCGCCCAGCAGGTCATAAat                                                                >  2:624306/1‑149 (MQ=255)
                                                                                      ggAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCAt                                                           <  1:745561/149‑1 (MQ=255)
                                                                                       gAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCAtt                                                          >  1:167196/1‑149 (MQ=255)
                                                                                       gAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGATTATTCAGCAGGAAGACAAGGGACAATTACTGACTCTTGCCGGAGAGCCTGACGGAAAAATCCGCGCAGCAGGTCCTTAATCTCAtt                                                          >  2:353218/1‑149 (MQ=255)
                                                                                         aaaTTTCCGAAGAAGCCACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCt                                                        <  2:1132144/149‑1 (MQ=255)
                                                                                             ttCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCGGACCATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaat                                                    >  2:1106647/1‑149 (MQ=255)
                                                                                              tCCGAAGAGTACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaatt                                                   <  2:22411/149‑1 (MQ=255)
                                                                                              tCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaatt                                                   >  1:805766/1‑149 (MQ=255)
                                                                                                cGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaattaa                                                 <  2:1139928/149‑1 (MQ=255)
                                                                                                cGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACCATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCttaattaa                                                 >  2:1350770/1‑149 (MQ=255)
                                                                                                       gaCACCACCGAAACGCCGATGTATGTCAATCTCGGTGGTCGCCGGATCATTAAAAAAAAAAACCAGGGACAATTACTGACTCTTGCCGGTGAGCGAGACGGAAACCTCCCCGCAGCAGGCCATTAATACCGTTCGTAATTAACTAatta                                          >  2:39833/1‑149 (MQ=255)
                                                                                                         caccacCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATAGTTCTTAATTCACTAattatt                                        >  1:188991/1‑149 (MQ=255)
                                                                                                            caccGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAAttatttat                                     <  1:1235480/149‑1 (MQ=255)
                                                                                                            caccGAAACGCCGATGTATGTCAATCGCGGTGGTCCCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAAttatttat                                     <  2:110051/149‑1 (MQ=255)
                                                                                                              ccGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTAtct                                   <  1:15716/149‑1 (MQ=255)
                                                                                                               cGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAACCTTACAGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTAtctc                                  <  1:279821/149‑1 (MQ=255)
                                                                                                                   tcGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCtcatca                              <  1:104474/148‑1 (MQ=255)
                                                                                                                        gATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGaa                         <  1:1255396/149‑1 (MQ=255)
                                                                                                                          tgtatgtCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAata                       >  1:316489/1‑149 (MQ=255)
                                                                                                                              tgtCAATCTCGGTGGTCGCCGTACCATTAAAAAAAAAAACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAACAACtt                   >  2:621459/1‑149 (MQ=255)
                                                                                                                              tgtCAATCTCGGTGGTCGCCGGATCATTAAAAAAAAAAACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACCCCGCGCAGCAAGACATTAATATCATTCTTAATTAACTAATAATTTATCTCATCACTGAATATCtt                   >  2:621521/1‑149 (MQ=255)
                                                                                                                               gtCAATCGCGGTGGTCGGCGGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTa                  <  1:302664/149‑1 (MQ=255)
                                                                                                                                tCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTaa                 <  1:171929/149‑1 (MQ=255)
                                                                                                                                        ggtggtCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGa         <  1:1225433/149‑1 (MQ=255)
                                                                                                                                        ggtggtCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCAGCACTGAATATCTTAATATAGTGa         <  1:809889/149‑1 (MQ=255)
                                                                                                                                           ggtCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAGAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGa      <  1:941673/149‑1 (MQ=255)
                                                                                                                                           ggtCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGa      <  1:288285/149‑1 (MQ=255)
                                                                                                                                               gCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACtta  >  2:47135/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTCACTCAACCTGGAAATCGAACCATTTGATGAGAACCGTCTGAAAATCAAACACAAACTCTCATATGTTCGCCCGACTAACCGTGGGAAAATTTCCGAAGAAGACACCACCGAAACGCCGATGTATGTCAATCGCGGTGGTCGCCTGACTATTCTGCAGGAAGACCAGGGACAATTACTGACTCTTGCCGGTGAGCCTGACGGAAAACTCCGCGCAGCAGGTCATTAATATCATTCTTAATTAACTAATTATTTATCTCATCACTGAATATCTTAATATAGTGAGGACTTA  >  NC_000913/2473736‑2474027

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