Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 270,206 T→C 6.0% intergenic (‑24/‑72) perR ← / → insN CP4‑6 prophage; putative DNA‑binding transcriptional regulator/pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913270,2060TC 6.0% 207.8 / 3.2 66intergenic (‑24/‑72)perR/insNCP4‑6 prophage; putative DNA‑binding transcriptional regulator/pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/2);  total (32/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 = 8.82e-01

ATGGCATGGCTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTACACCGT  >  NC_000913/270058‑270351
                                                                                                                                                    |                                                                                                                                                 
atggcatggcTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGAt                                                                                                                                                   >  2:134800/1‑149 (MQ=255)
   gcatggcTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGaa                                                                                                                                                >  2:677857/1‑149 (MQ=255)
     atggcTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGt                                                                                                                                              >  1:668874/1‑149 (MQ=255)
      tggcTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTg                                                                                                                                             <  2:475141/149‑1 (MQ=255)
        gcTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGt                                                                                                                                           <  4:37360/149‑1 (MQ=255)
          tGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGttt                                                                                                                                         <  2:567628/149‑1 (MQ=255)
           gATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTa                                                                                                                                        <  1:226059/149‑1 (MQ=255)
                    gggTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTg                                                                                                                               <  1:757116/149‑1 (MQ=255)
                    gggTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTg                                                                                                                               <  1:413485/149‑1 (MQ=255)
                      gTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGc                                                                                                                             <  4:295607/149‑1 (MQ=255)
                           aGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCt                                                                                                                        <  3:313011/149‑1 (MQ=255)
                           aGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCt                                                                                                                        <  1:357397/149‑1 (MQ=255)
                             ttCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGa                                                                                                                      >  1:384014/1‑149 (MQ=255)
                             ttCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGa                                                                                                                      >  3:327405/1‑149 (MQ=255)
                              tCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCCAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGaa                                                                                                                     <  2:668874/149‑1 (MQ=255)
                               ccAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAAc                                                                                                                    >  3:66372/1‑149 (MQ=255)
                                cAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACa                                                                                                                   <  2:384014/149‑1 (MQ=255)
                                   cTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACagag                                                                                                                <  1:366012/149‑1 (MQ=255)
                                   cTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACagag                                                                                                                <  2:77345/149‑1 (MQ=255)
                                    tCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAgg                                                                                                               >  3:292139/1‑149 (MQ=255)
                                              ccAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGGATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTc                                                                                                     <  2:109109/149‑1 (MQ=255)
                                              ccAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTc                                                                                                     >  2:289179/1‑149 (MQ=255)
                                               cAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCa                                                                                                    >  3:668656/1‑149 (MQ=255)
                                                aaCGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACagag                                                                                                                <  1:569548/136‑1 (MQ=255)
                                                aaCGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACagag                                                                                                                >  2:569548/1‑136 (MQ=255)
                                                        aGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaac                                                                                           <  4:775580/149‑1 (MQ=255)
                                                        aGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaac                                                                                           <  4:534736/149‑1 (MQ=255)
                                                         gCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaaca                                                                                          >  1:732472/1‑149 (MQ=255)
                                                              ggTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAAcaca                                                                                     <  1:677857/149‑1 (MQ=255)
                                                              ggTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAAcaca                                                                                     <  1:289179/149‑1 (MQ=255)
                                                                 aCGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgt                                                                                  <  4:327405/149‑1 (MQ=255)
                                                                        gCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGtt                                                                                                                                          <  2:28595/86‑1 (MQ=255)
                                                                        gCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGtt                                                                                                                                          >  1:28595/1‑86 (MQ=255)
                                                                        gCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACGACACAGGTGTCATaa                                                                           <  1:69528/149‑1 (MQ=255)
                                                                        gCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATaa                                                                           >  3:354282/1‑149 (MQ=255)
                                                                        gCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATaa                                                                           >  1:34099/1‑149 (MQ=255)
                                                                                    gcgcGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGaa                                                               >  2:769606/1‑149 (MQ=255)
                                                                                                       gTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTaa                                            >  2:596277/1‑149 (MQ=255)
                                                                                                           ttttGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTCGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAAcgc                                        >  2:798949/1‑149 (MQ=255)
                                                                                                               gCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAAt                                    <  4:155112/149‑1 (MQ=255)
                                                                                                               gCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAAt                                    >  3:155112/1‑149 (MQ=255)
                                                                                                                     agCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTc                              <  1:134800/149‑1 (MQ=255)
                                                                                                                           attcatGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaa                                                                                            >  2:531997/1‑81 (MQ=255)
                                                                                                                           attcatGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaa                                                                                            <  1:531997/81‑1 (MQ=255)
                                                                                                                           attcatGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTgg                        >  1:472459/1‑149 (MQ=255)
                                                                                                                           attcatGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTgg                        >  2:748771/1‑149 (MQ=255)
                                                                                                                               atGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGt                    >  3:190635/1‑149 (MQ=255)
                                                                                                                                tGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGtt                   >  4:472046/1‑149 (MQ=255)
                                                                                                                                tGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGtt                   <  2:34099/149‑1 (MQ=255)
                                                                                                                                 gAATCCTTTTCACTAGGGACGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGata                                                                                                          <  1:110197/61‑1 (MQ=255)
                                                                                                                                 gAATCCTTTTCACTAGGGACGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGata                                                                                                          >  2:314077/1‑61 (MQ=255)
                                                                                                                                 gAATCCTTTTCACTAGGGACGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGata                                                                                                          >  2:110197/1‑61 (MQ=255)
                                                                                                                                 gAATCCTTTTCACTAGGGACGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGata                                                                                                          <  1:314077/61‑1 (MQ=255)
                                                                                                                                   aTCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGAc                <  3:472046/149‑1 (MQ=255)
                                                                                                                                    tCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGAcc               >  2:325152/1‑149 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAAcgc                                        >  2:285868/1‑121 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAAcgc                                        <  1:285868/121‑1 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGc                                <  2:109403/129‑1 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGc                                >  1:109403/1‑129 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGAc                <  1:773177/145‑1 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGAc                >  2:773177/1‑145 (MQ=255)
                                                                                                                                       ttttCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGa            >  3:614925/1‑149 (MQ=255)
                                                                                                                                           cACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTa        >  3:790537/1‑149 (MQ=255)
                                                                                                                                            aCTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTac       >  1:804954/1‑149 (MQ=255)
                                                                                                                                             cTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTaca      <  4:790537/149‑1 (MQ=255)
                                                                                                                                               aGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTacaca    >  4:514670/1‑148 (MQ=255)
                                                                                                                                                 ggATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTACACCGt  <  4:66372/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGGCATGGCTGATCGCACTGGGTGACAGTTCCAACTCAGAGGCGGCCAACGCAAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTCGTTGACCAGAACTACACCGT  >  NC_000913/270058‑270351

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