Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 289,256 G→C 5.5% intergenic (‑94/+45) yagI ← / ← argF CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913289,2560GC 5.5% 246.0 / 5.2 73intergenic (‑94/+45)yagI/argFCP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base C (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

TTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCG  >  NC_000913/289110‑289401
                                                                                                                                                  |                                                                                                                                                 
ttGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGCAAAATCAATGAAAGGGAGCCGTTGTTCTCCCTCTCCAtgtg                                                                                                                                                 >  6:124054/1‑149 (MQ=255)
ttGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                 >  8:248215/1‑149 (MQ=255)
  gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTgg                                                                                                                                                <  8:353658/148‑1 (MQ=255)
  gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTgg                                                                                                                                                <  3:42326/148‑1 (MQ=255)
  gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGATGTCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTgg                                                                                                                                                <  2:300701/148‑1 (MQ=255)
     cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  2:259432/149‑1 (MQ=255)
       ggTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCACAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                          <  2:427708/149‑1 (MQ=255)
        gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                         >  7:133179/1‑149 (MQ=255)
        gTCGACGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATACTGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCGCCCTCTCCCACAGGGAGAggg                                                                                                                                         <  6:14872/149‑1 (MQ=255)
         tCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                        >  5:403790/1‑149 (MQ=255)
          cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                       >  4:188136/1‑149 (MQ=255)
             ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                    <  6:403790/149‑1 (MQ=255)
             ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCAAGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                    <  1:66105/149‑1 (MQ=255)
              gATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTTGGAGAGGGCCgggg                                                                                                                                   <  3:71974/149‑1 (MQ=255)
                    caccGCACGGTCACCAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTACGGCGAAAACCTCCCTTTTTCGGGATGAAAAATCAATGACATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAg                                                                                                                               <  1:311915/144‑1 (MQ=255)
                    cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                <  5:109547/146‑1 (MQ=255)
                    cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                             >  4:119457/1‑149 (MQ=255)
                    cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                             <  3:293542/149‑1 (MQ=255)
                      aCGCACGTTCAACAGACTGAATAATCGGCATAATGCCGGCATGTCCACGCTGGAAGTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTGCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                           <  1:112813/149‑1 (MQ=255)
                             ttCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgc                                                                                                                    >  3:131822/1‑149 (MQ=255)
                                aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTATTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAc                                                                                                                          >  7:265520/1‑140 (MQ=255)
                                  cAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGCGAGGGCACCAACGCGCAGc                                                                                                               <  2:396251/149‑1 (MQ=255)
                                     aCTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGa                                                                                                            >  4:193050/1‑149 (MQ=255)
                                       tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc                                                                                                          >  3:10946/1‑149 (MQ=255)
                                       tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc                                                                                                          >  1:308620/1‑149 (MQ=255)
                                         aataatCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAccc                                                                                                         >  6:145784/1‑148 (MQ=255)
                                         aataatCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAccc                                                                                                         >  1:95783/1‑148 (MQ=255)
                                              tCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCAc                                                                                                   >  2:344515/1‑149 (MQ=255)
                                              tCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCAc                                                                                                   >  5:321983/1‑149 (MQ=255)
                                                 gCATATTGGTGTCATGCCCACGCTTTAATTAAGGCTAAAACGTACCTTTTTGGGTATTAAAAATCAATGAATTGTAGCGGGTGTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGTGGGCACCAACGCGCAGCGGACCCAATCACTcc                                                                                                <  6:358632/149‑1 (MQ=255)
                                                 gCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCATTGAACTGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGCCACCAACGCGCAGCGGACCCAATCACTcc                                                                                                <  6:297184/149‑1 (MQ=255)
                                                 gCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                                                                                 >  7:1846/1‑148 (MQ=255)
                                                    tAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcccc                                                                                              >  7:353663/1‑148 (MQ=255)
                                                               tGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                     <  6:185246/98‑1 (MQ=255)
                                                               tGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                     >  5:185246/1‑98 (MQ=255)
                                                               tGTCCACGCTGCAATTAAGGCCAAAACCTACCTTTTTCGGTATTAAAAATCAATCAAATGCACCCGGTGTGCACCCTCTCCCTGTGGTATAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcat                                                                                  <  1:47018/149‑1 (MQ=255)
                                                                  ccACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                                                                                 >  1:401028/1‑131 (MQ=255)
                                                                  caccgcgGGAATTACGGCGAAAACGAACCTTTTTGGGGATGAAAACTCCTTCAATGGTCGCCGGTGTGCCCCCACTCCGTGTGGGTGTGGGGGGGGGGGTGGGCACCAACGCGCAGCGGACCCACTCACTCCCCAAGCGTTGCcatcat                                                                               <  5:124054/146‑1 (MQ=255)
                                                                    ccGCTGTAATTATGGCTAAAACGTACCTTTTTGGGGATTAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                                                                                 <  2:29627/128‑1 (MQ=38)
                                                                    aCGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                                                                                 >  1:29627/1‑129 (MQ=255)
                                                                     cGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAc                                                                            <  3:119457/149‑1 (MQ=255)
                                                                          gAATTAAGGCGAAAACGTACCTTTTTCGGGATTAAAAATCAATCAAATGGAGCCGGTGTTCTCCCTCACCCTGGGGGAGAGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc                                                                        <  3:193050/148‑1 (MQ=255)
                                                                          gAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCt                                                                       >  6:149122/1‑149 (MQ=255)
                                                                              ccgatataaaaCTAACTTTTTTCGGTATTAAAAATCAATCAATTGGCGCCGGTGTGCTCCCTCTCCCTGGGGGTGTGGGCCGGGGTGTGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTaa                                                                    <  2:5265/141‑1 (MQ=255)
                                                                                  gccaaaaGGTACCTTTTTCGGGATCAAAAATCAATGAAAGGGCGCCGGTTTGCCCCCACCCCCTGGGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                               <  5:145784/146‑1 (MQ=255)
                                                                                  gCGAAATCGTACCTTTTTCGGGATTAAAAATCAATTTAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTTAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                               <  2:308620/149‑1 (MQ=39)
                                                                                  gCGAAAACGTACCTTTTTCGGGATCAAAAATCAATCAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                               <  3:90950/149‑1 (MQ=255)
                                                                                   cGAAAACGTACCTTTTTCGGTATTAAAAATCAATGAAATGTAGCCGGTGTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                              <  4:180716/149‑1 (MQ=255)
                                                                                     aaaaCGTACCTTTTTCGGGATGAAAAATCGATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATg                                                             >  4:287333/1‑148 (MQ=255)
                                                                                      aaaCGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTaa                                                                    >  1:5265/1‑140 (MQ=255)
                                                                                          gTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCCCCGCCTTAATCGTATGCATGCg                                                       >  6:1314/1‑149 (MQ=255)
                                                                                          gTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCg                                                       >  7:74034/1‑149 (MQ=255)
                                                                                              cTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGttt                                                   >  3:380904/1‑149 (MQ=255)
                                                                                                        aTGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGt                                         <  1:84313/149‑1 (MQ=255)
                                                                                                         tGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                        <  1:344515/149‑1 (MQ=255)
                                                                                                          gAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGTGCCGGGGTGAGGGCACCTACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                        >  3:260329/1‑148 (MQ=255)
                                                                                                            aaaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGaa                                     >  4:238219/1‑149 (MQ=255)
                                                                                                              aaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCAAGCGTTGCCATCATGACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAaca                                   >  5:48420/1‑149 (MQ=255)
                                                                                                                tCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGt                                         >  2:48007/1‑141 (MQ=255)
                                                                                                                tCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGt                                         <  1:48007/141‑1 (MQ=255)
                                                                                                                  aattaattGGCGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                               <  8:353663/149‑1 (MQ=255)
                                                                                                                  aaTGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCCACGCGCAGCGGACCCCATCACTCCCCCCGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                               >  1:165712/1‑149 (MQ=255)
                                                                                                                  aaTGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGAt                               >  1:171798/1‑149 (MQ=255)
                                                                                                                              ccGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGa                   >  1:96300/1‑149 (MQ=255)
                                                                                                                                    tgctccctctcccgGGGGTCGTGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCaaa             <  2:95783/147‑1 (MQ=255)
                                                                                                                                     tgtccctctccctGTGGTAGTGGGCCGGGGGGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCaaa             <  2:96300/146‑1 (MQ=255)
                                                                                                                                      ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAaca           >  8:174294/1‑149 (MQ=255)
                                                                                                                                       tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAacac          >  3:124479/1‑149 (MQ=255)
                                                                                                                                        ccctctccctGTGGTAGTGGGCGGGGGTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACAcc         <  4:205296/149‑1 (MQ=255)
                                                                                                                                        ccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCAACAACGCGCAGCGGACCCAATCACTACCCAATCGTTGCAATCATAACCGACTTAATCGTATGCATCAGGTTTTCCGCCTGGTCGAACACGATGATCGCAGCAGACTAAAaaac          >  6:223152/1‑148 (MQ=255)
                                                                                                                                         cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCTCTCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt        >  1:217856/1‑149 (MQ=255)
                                                                                                                                         cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACAcc         >  3:205296/1‑148 (MQ=255)
                                                                                                                                         cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt        >  6:310966/1‑149 (MQ=255)
                                                                                                                                           tctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcg      >  4:241855/1‑149 (MQ=255)
                                                                                                                                               cctGTGGGAGAGGGCCGGGGTGATGGCACAAACGCGCAGCGGACCCAAACACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCCCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtccg  >  1:219087/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCG  >  NC_000913/289110‑289401

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: