Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,758,920 C→T 5.4% intergenic (+64/‑65) alpA → / → yfjI CP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,758,9200CT 5.4% 370.3 / 7.6 111intergenic (+64/‑65)alpA/yfjICP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (48/57);  new base T (3/3);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 115 total aligned reads displayed.

GAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCA  >  NC_000913/2758777‑2759067
                                                                                                                                               |                                                                                                                                                    
gAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                                                                 >  3:930926/1‑149 (MQ=255)
  aTGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAAt                                                                                                                                                <  1:370777/148‑1 (MQ=255)
  aTGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATa                                                                                                                                               <  4:678889/149‑1 (MQ=255)
      tCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCgg                                                                                                                                           >  3:818670/1‑149 (MQ=255)
       cGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGa                                                                                                                                         <  4:559573/149‑1 (MQ=255)
          cGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAAt                                                                                                                                       <  2:514726/148‑1 (MQ=255)
          cGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAAt                                                                                                                                       <  2:334976/148‑1 (MQ=255)
              aaaTCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaatagc                                                                                                                                   <  1:324843/148‑6 (MQ=255)
              aaaTCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaatagc                                                                                                                                   <  1:471219/148‑6 (MQ=255)
                 tCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCata                                                                                                                                >  2:323116/1‑148 (MQ=255)
                       ggatggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                                                                 <  4:400747/126‑1 (MQ=255)
                       ggatggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                                                                 >  3:400747/1‑126 (MQ=255)
                       ggatggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTa                                                                                                                         >  3:508726/1‑149 (MQ=255)
                        gatggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTa                                                                                                                         >  2:290704/1‑148 (MQ=255)
                        gatggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTa                                                                                                                         >  1:120005/1‑148 (MQ=255)
                          tggatTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTaaa                                                                                                                       <  2:929547/148‑1 (MQ=255)
                               tGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAttt                                                                                                                 <  4:644201/149‑1 (MQ=255)
                                gAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAttt                                                                                                                 <  1:915883/148‑1 (MQ=255)
                                   gCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCtt                                                                                                              <  1:323116/148‑1 (MQ=255)
                                     cGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTaaa                                                                                                                       >  3:535220/1‑137 (MQ=255)
                                     cGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaa                                                                                                            <  1:570072/148‑1 (MQ=255)
                                     cGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaa                                                                                                            <  1:712401/148‑1 (MQ=255)
                                      gagaTTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatt                                                                                                          >  3:18647/1‑149 (MQ=255)
                                      gagaTTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaat                                                                                                           <  1:75367/148‑1 (MQ=255)
                                         aTTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                                                                       <  3:895685/149‑1 (MQ=255)
                                          ttGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaat                                                                                                      <  3:47992/149‑1 (MQ=255)
                                          ttGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                                                                       <  1:195170/148‑1 (MQ=255)
                                           tgatgaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaat                                                                                                      <  1:622264/148‑1 (MQ=255)
                                                  tGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAAt                                                                                                                                                >  2:625111/1‑100 (MQ=255)
                                                  tGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAAt                                                                                                                                                <  1:625111/100‑1 (MQ=255)
                                                  tGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGc                                                                                                                           >  3:519920/1‑120 (MQ=255)
                                                  tGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGc                                                                                                                           <  4:519920/120‑1 (MQ=255)
                                                  tGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                                                               >  1:739017/1‑148 (MQ=255)
                                                   ggTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAaga                                                                                             <  3:70574/149‑1 (MQ=255)
                                                   ggTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAag                                                                                              <  1:697843/148‑1 (MQ=255)
                                                    gTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                                                            >  4:676941/1‑149 (MQ=255)
                                                    gTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                                                            <  3:924665/149‑1 (MQ=255)
                                                         cacaACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCAt                                                                                        <  2:823567/148‑1 (MQ=255)
                                                           caACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCAt                                                                                        <  2:685606/146‑1 (MQ=255)
                                                            aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGc                                                                                     >  1:242255/1‑148 (MQ=255)
                                                            aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                                                                    >  4:718796/1‑149 (MQ=255)
                                                               gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGt                                                                                 >  4:41485/1‑149 (MQ=255)
                                                                cTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGtt                                                                                <  4:500849/149‑1 (MQ=255)
                                                                 tGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTc                                                                               >  3:605611/1‑149 (MQ=255)
                                                                  gTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                                                                              >  3:208079/1‑149 (MQ=255)
                                                                   tAAACTTATGTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCTCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                                                                              <  2:196054/148‑1 (MQ=255)
                                                                    aaaCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                                                            <  3:460817/149‑1 (MQ=255)
                                                                      aCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                                                                          >  3:545404/1‑149 (MQ=255)
                                                                       cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAg                                                                                                                                              >  4:296913/1‑81 (MQ=255)
                                                                            ttGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTGAATGgtcgtc                                                                     <  1:506038/148‑1 (MQ=255)
                                                                            ttGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTcc                                                                    <  4:545404/149‑1 (MQ=255)
                                                                             tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCt                                                                   <  3:360220/149‑1 (MQ=255)
                                                                                        cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTcc                                                               >  4:208180/1‑142 (MQ=255)
                                                                                        cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTcc                                                               <  3:208180/142‑1 (MQ=255)
                                                                                         aTGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGGATGTTCAATGGTCGTCCTTTCCCTGTAGAt                                                       >  3:776937/1‑149 (MQ=255)
                                                                                          tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt                                                       >  2:909208/1‑148 (MQ=255)
                                                                                          tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt                                                       >  1:107138/1‑148 (MQ=255)
                                                                                           gAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATg                                                      >  2:842836/1‑148 (MQ=255)
                                                                                                 aTACTCAATTGCTGTCATCGGAGTGTAAACATTCTGCTATTTCCACCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATCTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTcc                                               <  3:103633/149‑1 (MQ=255)
                                                                                                 aTACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTcc                                               <  4:618014/149‑1 (MQ=255)
                                                                                                   aCTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTa                                             >  3:536376/1‑149 (MQ=255)
                                                                                                    cTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaa                                            <  3:914367/149‑1 (MQ=255)
                                                                                                     tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAtt                                                                                                                  <  4:804886/78‑1 (MQ=255)
                                                                                                     tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAtt                                                                                                                  >  3:804886/1‑78 (MQ=255)
                                                                                                     tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaaa                                           >  3:323251/1‑149 (MQ=255)
                                                                                                     tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAAGTTAATTAAGTAAGAGGCCTTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaaa                                           >  4:775011/1‑149 (MQ=255)
                                                                                                      cAATTGTTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaaa                                           >  1:839621/1‑148 (MQ=255)
                                                                                                      cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaaa                                           >  1:836953/1‑148 (MQ=255)
                                                                                                      cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTaaa                                           >  1:553798/1‑148 (MQ=255)
                                                                                                       aaTTGCTGTCATCGTAGTGTAAACATCCTGCTATTTCCGCCAGTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTaaa                                                                                                  <  4:2318/92‑1 (MQ=255)
                                                                                                       aaTTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTaaa                                                                                                  >  3:2318/1‑92 (MQ=255)
                                                                                                        aTTGCTGTCATCGGAGTGTAAACATCCTACTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTGTTAATTTAATTAAATAAGAGGCCTTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCAGTTCCTAAAAtt                                        >  4:910020/1‑149 (MQ=255)
                                                                                                           gCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATc                                     <  3:676941/149‑1 (MQ=255)
                                                                                                            cTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATc                                     <  1:207265/148‑1 (MQ=255)
                                                                                                            cTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCACTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATc                                     >  1:295805/1‑148 (MQ=255)
                                                                                                             tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAg                                   <  3:775011/149‑1 (MQ=255)
                                                                                                                cATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGa                                 >  2:560886/1‑148 (MQ=255)
                                                                                                                  tCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAAt                               >  1:239698/1‑148 (MQ=255)
                                                                                                                  tCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATg                              >  3:613070/1‑149 (MQ=255)
                                                                                                                  tCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATg                              >  3:122850/1‑149 (MQ=255)
                                                                                                                      aGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAAt                          <  3:695807/149‑1 (MQ=255)
                                                                                                                          tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTAt                      <  4:613070/149‑1 (MQ=255)
                                                                                                                           aaaCATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTCCCCGGTAGATGCATTTCCTAAAATTATCAGGAATGAAATTTAt                      >  1:452170/1‑148 (MQ=255)
                                                                                                                             aCATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGa                    <  2:502925/148‑1 (MQ=255)
                                                                                                                             aCATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGa                    <  2:744402/148‑1 (MQ=255)
                                                                                                                                  cTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGa              >  4:78963/1‑149 (MQ=255)
                                                                                                                                   tGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                                                            >  2:61398/1‑86 (MQ=255)
                                                                                                                                   tGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                                                            <  1:61398/86‑1 (MQ=255)
                                                                                                                                   tGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGaa             <  3:853075/149‑1 (MQ=255)
                                                                                                                                     gctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAAc            >  3:359249/5‑149 (MQ=255)
                                                                                                                                      tATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCa                                    >  2:343859/1‑123 (MQ=255)
                                                                                                                                      tATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCa                                    <  1:343859/123‑1 (MQ=255)
                                                                                                                                       tattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                                                             >  3:7194/3‑82 (MQ=37)
                                                                                                                                       tattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                                                             <  4:7194/80‑1 (MQ=37)
                                                                                                                                       aTTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAg          >  2:88351/1‑148 (MQ=255)
                                                                                                                                       aTTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAg          <  2:49193/148‑1 (MQ=255)
                                                                                                                                          tCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCAt       >  1:384202/1‑148 (MQ=255)
                                                                                                                                              ccATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCACTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGc   <  2:295805/148‑1 (MQ=255)
                                                                                                                                              ccATTAATAGCGGAATAGCCATATATGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGc   <  1:842836/148‑1 (MQ=255)
                                                                                                                                               cATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCa  >  2:821980/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GAATGCTCGGCGTGAAATCTGTCGGATGGATTGAGGCCGAGATTGATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAAATTATCAGGAATGCAATTTATGAAGTGGAACAGCATACGCA  >  NC_000913/2758777‑2759067

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