Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 292,302 C→G 7.4% intergenic (+71/+20) yagJ → / ← yagK CP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,3020CG 7.4% 331.8 / 15.5 94intergenic (+71/+20)yagJ/yagKCP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base C (48/39);  new base G (4/3);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

GATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGT  >  NC_000913/292154‑292450
                                                                                                                                                    |                                                                                                                                                    
gATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                                                                      <  1:204114/149‑1 (MQ=255)
 aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACt                                                                                                                                                     <  1:119001/149‑1 (MQ=255)
      cAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGaa                                                                                                                                                <  2:142911/149‑1 (MQ=255)
      cAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGaa                                                                                                                                                <  1:59777/149‑1 (MQ=255)
       aGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                                                                               <  2:76069/149‑1 (MQ=255)
        ggATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                                                                              >  4:298339/1‑149 (MQ=255)
         gATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  4:52006/149‑1 (MQ=255)
                  gCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGa                                                                                                                                    <  2:587062/149‑1 (MQ=255)
                   cTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGAc                                                                                                                                   <  3:468719/149‑1 (MQ=255)
                     aCCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                                                                                 >  3:776166/1‑149 (MQ=255)
                                tGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGAtt                                                                                                                      >  3:41522/1‑149 (MQ=255)
                                  cAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTg                                                                                                                    <  3:194360/149‑1 (MQ=255)
                                   aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  1:225356/123‑1 (MQ=255)
                                   aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             >  2:225356/1‑123 (MQ=255)
                                   aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   >  4:408688/1‑149 (MQ=255)
                                   aGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   <  3:298339/149‑1 (MQ=255)
                                        cGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACa                                                                                                              >  2:794195/1‑149 (MQ=255)
                                         gTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAg                                                                                                             >  3:70604/1‑149 (MQ=255)
                                              ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGGCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   >  1:139353/1‑138 (MQ=255)
                                              ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   <  2:139353/138‑1 (MQ=255)
                                              ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGaa                                                                                                        >  3:465668/1‑149 (MQ=255)
                                              ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGaa                                                                                                        >  2:494105/1‑149 (MQ=255)
                                              ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGaa                                                                                                        <  1:794195/149‑1 (MQ=255)
                                                 gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                                                                               >  4:729376/1‑107 (MQ=255)
                                                 gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                                                                               <  3:729376/107‑1 (MQ=255)
                                                 gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGtt                                                                                                     >  1:9760/1‑149 (MQ=255)
                                                              tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                                                                        >  1:330335/1‑149 (MQ=255)
                                                              tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                                                                        >  3:485664/1‑149 (MQ=255)
                                                              tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                                                                        >  2:495652/1‑149 (MQ=255)
                                                              tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                                                                        >  4:732825/1‑149 (MQ=255)
                                                               gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAc                                                                                       <  3:732825/149‑1 (MQ=255)
                                                                   aCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGAGTTACAGTTCTATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCaaa                                                                                   >  4:59833/1‑149 (MQ=255)
                                                                   aCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCA‑TTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAac                                                                                  >  1:570809/1‑149 (MQ=255)
                                                                       tCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACt                                                                               >  4:136689/1‑149 (MQ=255)
                                                                           tAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGc                                                                           >  2:771179/1‑149 (MQ=255)
                                                                            aaTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCt                                                                          <  1:495652/149‑1 (MQ=255)
                                                                             aTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCtt                                                                         >  4:156042/1‑149 (MQ=255)
                                                                                 ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCaa                                                                     <  1:494105/149‑1 (MQ=255)
                                                                                  taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  3:567004/1‑149 (MQ=255)
                                                                                  taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  4:659222/1‑149 (MQ=255)
                                                                                  taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  4:386806/1‑149 (MQ=255)
                                                                                  taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  3:760494/1‑149 (MQ=255)
                                                                                  taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  1:485074/1‑149 (MQ=255)
                                                                                  taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATCTGGCTACAGCTGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  3:525801/1‑149 (MQ=255)
                                                                                   aGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                                                                            >  2:737256/1‑76 (MQ=255)
                                                                                   aGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                                                                            <  1:737256/76‑1 (MQ=255)
                                                                                       aaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                               <  4:567004/149‑1 (MQ=255)
                                                                                       aaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATCTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                               <  4:525801/149‑1 (MQ=255)
                                                                                           gCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAg                                                           >  2:100663/1‑149 (MQ=255)
                                                                                              tAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTa                                                        <  3:591580/149‑1 (MQ=255)
                                                                                               aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                <  3:35447/140‑1 (MQ=255)
                                                                                               aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                >  4:35447/1‑140 (MQ=255)
                                                                                                       tGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCg                                               >  1:87812/1‑149 (MQ=255)
                                                                                                       tGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCg                                               >  2:350297/1‑149 (MQ=255)
                                                                                                        gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCgg                                              <  3:659222/149‑1 (MQ=255)
                                                                                                         cATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCggg                                             <  2:330335/149‑1 (MQ=255)
                                                                                                           ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGtt                                           <  3:136689/149‑1 (MQ=255)
                                                                                                           ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGtt                                           >  4:440951/1‑149 (MQ=255)
                                                                                                             ttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTaa                                         <  3:386806/149‑1 (MQ=255)
                                                                                                                 tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTaa                                         <  4:344739/145‑1 (MQ=255)
                                                                                                                 tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTaa                                         >  3:344739/1‑145 (MQ=255)
                                                                                                                 tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                     >  4:231962/1‑149 (MQ=255)
                                                                                                                 tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                     <  4:41522/149‑1 (MQ=255)
                                                                                                                 tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                     >  4:566581/1‑149 (MQ=255)
                                                                                                                    ttttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGcc                                  <  4:70604/149‑1 (MQ=255)
                                                                                                                     tttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCt                                 <  3:764963/149‑1 (MQ=255)
                                                                                                                      ttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAGGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCt                                 <  1:697130/148‑1 (MQ=255)
                                                                                                                      ttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAGGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCt                                 >  2:697130/1‑148 (MQ=255)
                                                                                                                      ttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTg                                <  1:380435/149‑1 (MQ=255)
                                                                                                                        ttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGt                              >  2:572317/1‑149 (MQ=255)
                                                                                                                           aaaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGt                           >  3:64750/1‑149 (MQ=255)
                                                                                                                            aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGtt                          <  4:760494/149‑1 (MQ=255)
                                                                                                                                aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      <  4:485664/149‑1 (MQ=255)
                                                                                                                                aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      >  3:504620/1‑149 (MQ=255)
                                                                                                                                 tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAg                                                           >  1:528959/1‑111 (MQ=255)
                                                                                                                                 tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAg                                                           <  2:528959/111‑1 (MQ=255)
                                                                                                                                 tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  1:386972/1‑149 (MQ=255)
                                                                                                                                 tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAATTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  4:554478/1‑149 (MQ=255)
                                                                                                                                     aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  1:100663/149‑1 (MQ=255)
                                                                                                                                     aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  1:350297/149‑1 (MQ=255)
                                                                                                                                     aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                 <  4:465668/149‑1 (MQ=255)
                                                                                                                                      aTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcgc                                   >  3:693067/1‑130 (MQ=255)
                                                                                                                                      aTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcgc                                   <  4:693067/130‑1 (MQ=255)
                                                                                                                                      aTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAg                >  1:259788/1‑149 (MQ=255)
                                                                                                                                      aTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAg                >  3:763472/1‑149 (MQ=255)
                                                                                                                                       tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGt               <  3:156042/149‑1 (MQ=255)
                                                                                                                                         aGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGCAGTTCTGCTGGAAGTCa             <  4:776166/149‑1 (MQ=255)
                                                                                                                                           gCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGt           >  1:259440/1‑149 (MQ=255)
                                                                                                                                           gCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGt           >  1:175975/1‑149 (MQ=255)
                                                                                                                                             cTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCggg                                             <  4:189963/113‑1 (MQ=255)
                                                                                                                                             cTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCggg                                             >  3:189963/1‑113 (MQ=255)
                                                                                                                                              ttCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACt        <  2:570809/149‑1 (MQ=255)
                                                                                                                                                  cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTg    >  3:241272/1‑149 (MQ=255)
                                                                                                                                                   aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg   >  1:55365/1‑149 (MQ=255)
                                                                                                                                                   aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg   >  1:118617/1‑149 (MQ=255)
                                                                                                                                                    cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTg    <  2:405123/147‑1 (MQ=255)
                                                                                                                                                    cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTg    >  1:405123/1‑147 (MQ=255)
                                                                                                                                                    cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGgt  <  3:440951/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGT  >  NC_000913/292154‑292450

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