Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 292,302 C→G 11.1% 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,3020CG11.1% 420.2 / 35.1 127intergenic (+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 (53/59);  new base G (6/8);  total (59/68)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 130 total aligned reads displayed.

ATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGT  >  NC_000913/292155‑292450
                                                                                                                                                   |                                                                                                                                                    
aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                                                                      <  2:37177/148‑1 (MQ=255)
aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGGAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACt                                                                                                                                                     >  3:179717/1‑149 (MQ=255)
  gtgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTga                                                                                                                                                   <  3:295091/149‑1 (MQ=255)
      aGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                                                                               <  3:236555/149‑1 (MQ=255)
        gacgCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                                                                             >  4:384205/4‑149 (MQ=255)
        gATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  3:926851/149‑1 (MQ=255)
        gATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  4:443745/149‑1 (MQ=255)
        gATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                                                                             >  3:747149/1‑149 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  2:560261/148‑1 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGc                                                                                                                                             <  1:816406/148‑1 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                                                                             <  1:484357/148‑1 (MQ=255)
                  cTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGa                                                                                                                                    <  2:875634/148‑1 (MQ=255)
                  cTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGAc                                                                                                                                   <  4:613487/149‑1 (MQ=255)
                     ccGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                                                                                 <  1:887846/148‑1 (MQ=255)
                      cGCTAAAAATGAAAAAACTTTATTAAGTACCGAACTTGCTAAGGCTAATAAGAAACTAAAAAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACa                                                                                                                              <  3:234358/149‑1 (MQ=255)
                          gAAAATGCACAAACGTTATTTAGTAGCGGACTTGCTGATGACTATAAGTAATCTATAGAAAATTGCGTAAAGAAACTGCATTTTAGTAATTTTTTGAAACCATACAAATAAGGCCTTCCCACTGAGAAGGCCATTACCGGTTTACAGtc                                                                                                                           <  4:243727/149‑2 (MQ=255)
                           aaaaTGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTc                                                                                                                          <  3:942761/149‑1 (MQ=255)
                           aaaaTGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTc                                                                                                                          <  4:726949/149‑1 (MQ=255)
                                   gAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                                                                           <  4:552557/124‑2 (MQ=255)
                                   gAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   >  1:689745/1‑148 (MQ=255)
                                    aaaCGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                   >  4:358317/1‑147 (MQ=255)
                                                   gCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTAc                                                                                                   <  2:946374/148‑1 (MQ=255)
                                                   gCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACg                                                                                                  >  3:462323/1‑149 (MQ=255)
                                                    cGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTg                                                                                                                    >  3:382745/1‑130 (MQ=255)
                                                    cGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACg                                                                                                  >  2:823808/1‑148 (MQ=255)
                                                     ggACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACgg                                                                                                 <  1:811523/148‑1 (MQ=255)
                                                      gACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGt                                                                                                >  1:426421/1‑148 (MQ=255)
                                                       aCTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtct                                                                                              >  3:176032/1‑149 (MQ=255)
                                                       aCTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtc                                                                                               >  1:246491/1‑148 (MQ=255)
                                                          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTa                                                                                                    <  1:473116/140‑1 (MQ=255)
                                                          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTa                                                                                                    >  2:473116/1‑140 (MQ=255)
                                                          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTcc                                                                                            >  2:404422/1‑148 (MQ=255)
                                                             tgatgaCTATCTGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                                                                        >  3:762544/1‑149 (MQ=255)
                                                              gatgaCTATCAGTAGTCTATAGAAAATTGCGTAACGATACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAc                                                                                       >  3:69650/1‑149 (MQ=255)
                                                              gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTcg                                                                                        >  1:791470/1‑147 (MQ=255)
                                                               atgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAcc                                                                                      <  3:353380/149‑1 (MQ=255)
                                                                 gaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCaa                                                                                    >  3:377886/1‑149 (MQ=255)
                                                                        aGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACttt                                                                             <  4:201324/149‑1 (MQ=255)
                                                                        aGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACtt                                                                              <  2:174255/148‑1 (MQ=255)
                                                                         gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGa                                                                                                         <  1:441657/120‑1 (MQ=255)
                                                                         gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGa                                                                                                         >  2:441657/1‑120 (MQ=255)
                                                                          tAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACa                                                                                                              >  1:249516/1‑114 (MQ=255)
                                                                          tAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACa                                                                                                              <  2:249516/114‑1 (MQ=255)
                                                                           aaTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGc                                                                           >  1:198010/1‑148 (MQ=255)
                                                                                ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCa                                                                      <  1:404422/148‑1 (MQ=255)
                                                                                ataGAAAATTGCGTAACGAAACTACATTTTAGTAATTTTTTTAAACAATACAAATAAGGCCTTCTCATTGAGGAGGCCATAACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCac                                                                     <  4:188937/149‑2 (MQ=255)
                                                                                 taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  3:151005/1‑149 (MQ=255)
                                                                                  aGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                    >  1:27580/1‑148 (MQ=255)
                                                                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACt                                                                   <  2:807497/148‑1 (MQ=255)
                                                                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACt                                                                   >  2:35585/1‑148 (MQ=255)
                                                                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACttt                                                                 >  4:672428/1‑149 (MQ=255)
                                                                                     aaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                <  3:5565/149‑1 (MQ=255)
                                                                                            gTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAgg                                                          >  1:637202/1‑148 (MQ=255)
                                                                                               aCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTa                                                                                                    >  2:371620/1‑103 (MQ=255)
                                                                                               aCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTa                                                                                                    <  1:371620/103‑1 (MQ=255)
                                                                                               aCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAg                                                       >  2:637015/1‑148 (MQ=255)
                                                                                                  aaaCTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCa                                                   <  4:460232/149‑1 (MQ=255)
                                                                                                   aaCTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCa                                                   <  2:198010/148‑1 (MQ=255)
                                                                                                    aCTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGt                                                 >  3:763344/1‑149 (MQ=255)
                                                                                                     cTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTc                                                >  4:367545/1‑149 (MQ=255)
                                                                                                       gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCg                                               >  1:42979/1‑148 (MQ=255)
                                                                                                        cATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCggg                                             <  4:798839/149‑1 (MQ=255)
                                                                                                          ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAg                                                           >  3:30778/1‑133 (MQ=255)
                                                                                                          ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGtt                                           <  4:561992/149‑1 (MQ=255)
                                                                                                          ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTCACAGTTCGATTTGGCTACAGCCGAAGTTTCGGTCTCCTTCGCCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGt                                            <  2:791470/148‑1 (MQ=255)
                                                                                                           tttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGtt                                           >  1:326076/1‑148 (MQ=255)
                                                                                                            ttAGTAATTTTTTGAAACAATGCAAATAAGGCCTTCTCACTGAGAATGCCAGTACCGACTTACAGTCCGATTTGGCTAGAGCCGAAGTTACGGTCGCCTTCACCGTACACTTTGCTTTTAACTTTGGTCAGGTAGTCGAGTCGGGTTaa                                         <  4:347337/149‑1 (MQ=255)
                                                                                                                tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                     <  3:740890/149‑1 (MQ=255)
                                                                                                                   ttttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTa                                          >  1:244328/1‑141 (MQ=255)
                                                                                                                   ttttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTa                                          <  2:244328/141‑1 (MQ=255)
                                                                                                                   ttttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcgc                                   >  2:77953/1‑148 (MQ=255)
                                                                                                                    tttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGcc                                  <  2:689745/148‑1 (MQ=255)
                                                                                                                        tGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGt                              <  2:704531/148‑1 (MQ=255)
                                                                                                                          aaaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGt                           <  4:428398/149‑1 (MQ=255)
                                                                                                                           aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                               <  1:295122/112‑1 (MQ=255)
                                                                                                                           aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                               >  2:295122/1‑112 (MQ=255)
                                                                                                                           aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGtt                          >  4:743789/1‑149 (MQ=255)
                                                                                                                            aCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGtt                          >  1:734141/1‑148 (MQ=255)
                                                                                                                               aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      >  3:270625/1‑149 (MQ=255)
                                                                                                                               aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctg                       <  1:757156/148‑1 (MQ=255)
                                                                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                      <  1:33588/148‑1 (MQ=255)
                                                                                                                                 aCAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAg                                                       <  2:308190/114‑1 (MQ=255)
                                                                                                                                 aCAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                     >  2:539293/1‑148 (MQ=255)
                                                                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCa                                                            <  4:720011/106‑1 (MQ=255)
                                                                                                                                       atggCCTACTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACTGTCTCCTTCACCAAACACTTTACTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTc              >  4:479287/3‑149 (MQ=255)
                                                                                                                                         ggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCa             >  4:573659/1‑148 (MQ=255)
                                                                                                                                         ggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCa             <  3:573659/148‑1 (MQ=255)
                                                                                                                                         ggCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAg            >  4:327928/1‑149 (MQ=255)
                                                                                                                                          gCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAg            >  1:99643/1‑148 (MQ=255)
                                                                                                                                            cTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTa          >  2:864227/1‑148 (MQ=255)
                                                                                                                                              tctcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACt        >  2:595252/1‑148 (MQ=255)
                                                                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAtt     >  4:814708/1‑149 (MQ=255)
                                                                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAtt     >  4:152699/1‑149 (MQ=255)
                                                                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAt      <  1:539293/148‑1 (MQ=255)
                                                                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAt      <  2:42979/148‑1 (MQ=255)
                                                                                                                                                 cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAtt     >  1:667463/1‑148 (MQ=255)
                                                                                                                                                 cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTg    <  4:462323/149‑1 (MQ=255)
                                                                                                                                                 cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTg    <  3:296849/149‑1 (MQ=255)
                                                                                                                                                   cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg   <  1:77953/148‑1 (MQ=255)
                                                                                                                                                   cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTCACCAAACACTTTTGTTTTAAATTTGGTTAGGTAGTTCAGTCGGGGTTATAGCGCCTGGTTGTTTTTCTGGAAGTCAGGACTATTGgt  >  1:279462/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGT  >  NC_000913/292155‑292450

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