Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,497,954 A→T 6.8% intergenic (+14/‑318) intB → / → insC‑6 KpLE2 phage‑like element; putative integrase/KpLE2 phage‑like element; IS2 insertion element repressor InsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,497,9540AT 6.8% 574.0 / 20.7 184intergenic (+14/‑318)intB/insC‑6KpLE2 phage‑like element; putative integrase/KpLE2 phage‑like element; IS2 insertion element repressor InsA
Reads supporting (aligned to +/- strand):  ref base A (93/78);  new base T (9/4);  total (102/82)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TATATCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACC  >  NC_000913/4497806‑4498102
                                                                                                                                                    |                                                                                                                                                    
tataTCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAAt                                                                                                                                                     >  2:1543231/1‑150 (MQ=255)
  tatCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGc                                                                                                                                                   >  2:1693545/1‑150 (MQ=255)
   atCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGcc                                                                                                                                                  <  2:603611/150‑1 (MQ=255)
   atCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGcc                                                                                                                                                  <  1:2266389/150‑1 (MQ=255)
    tCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGccc                                                                                                                                                 >  2:1098889/1‑150 (MQ=255)
     ccATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCg                                                                                                                                                >  2:1133147/1‑150 (MQ=255)
     ccATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCg                                                                                                                                                >  1:175281/1‑150 (MQ=255)
       aTAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCggg                                                                                                                                              >  2:591988/1‑150 (MQ=255)
         aaaGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGaa                                                                                                                                            <  1:1804094/150‑1 (MQ=255)
          aaGCTGAGCATCTTGAAGAACTGCGATTGATTCTACAGTGGTGGGCCGATTTTCTGGATTTAAACAGTGAAATTTTTATGAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCATTTGCCCGGGatg                                                                                                                                           <  1:2238861/150‑3 (MQ=255)
          aaGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAat                                                                                                                                           >  1:865091/1‑150 (MQ=255)
          aaGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAat                                                                                                                                           <  1:1893817/150‑1 (MQ=255)
               gagCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTc                                                                                                                                      <  2:1784500/150‑1 (MQ=255)
                  cATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAg                                                                                                                                   <  1:590541/150‑1 (MQ=255)
                      tGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGatt                                                                                                                               >  1:1215225/1‑150 (MQ=255)
                       ggAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGatta                                                                                                                              >  1:2439927/1‑150 (MQ=255)
                       ggAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGatta                                                                                                                              >  2:1293018/1‑150 (MQ=255)
                       ggAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGatta                                                                                                                              >  1:1167671/1‑150 (MQ=255)
                       ggAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGatta                                                                                                                              >  1:1866128/1‑150 (MQ=255)
                          agaaCGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTAttt                                                                                                                           >  2:449716/1‑150 (MQ=255)
                          agaaCGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTAttt                                                                                                                           >  2:2058661/1‑150 (MQ=255)
                              cGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCtt                                                                                                                       >  2:1322373/1‑150 (MQ=255)
                               ggCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCttt                                                                                                                      >  2:138319/1‑150 (MQ=255)
                               ggCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCttt                                                                                                                      <  2:2076909/150‑1 (MQ=255)
                               ggCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCttt                                                                                                                      <  2:1377560/150‑1 (MQ=255)
                                gCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTg                                                                                                                     <  2:2289357/150‑1 (MQ=255)
                                gCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTg                                                                                                                     <  2:2190453/150‑1 (MQ=255)
                                      gATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTaaaa                                                                                                               >  1:2105015/1‑150 (MQ=255)
                                       aTGCTACAGTGGTTGGCCGATTTTCTGTATGTAAACAGAGAAATGTTTATCAGTCCATTTTATTATTGAAATATTAATAATCCATTAAAACATTAATCATCCCGTGCAAATGCCCGGGAATTATTGTAGGATTATTTTCTTTGTTaaaaa                                                                                                              <  2:338199/150‑1 (MQ=255)
                                       aTGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTaaaaa                                                                                                              <  1:980054/150‑1 (MQ=255)
                                        tGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTaaaaaa                                                                                                             <  2:1344078/150‑1 (MQ=255)
                                        tGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTaaaaaa                                                                                                             <  1:684635/150‑1 (MQ=255)
                                              aGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACaaa                                                                                                       <  2:2277097/150‑1 (MQ=255)
                                                    gggCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTAt                                                                                                 <  1:1705613/150‑1 (MQ=255)
                                                     gggcGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTAtt                                                                                                >  2:1035667/4‑150 (MQ=255)
                                                     ggCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTAtt                                                                                                >  1:739213/1‑150 (MQ=255)
                                                      gacGATTTTCTGGATGTAAACAGAGAAAGGTATATCAGTCCATTTGAATATGAAAAGATTAATAATCCATTACACCACTAATCATCCCGGGCAAATGCACGGGAATTATTATAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTa                                                                                               >  2:1930179/3‑150 (MQ=255)
                                                       ccGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTaa                                                                                              >  2:875841/1‑150 (MQ=255)
                                                        cGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAAc                                                                                             <  2:1262822/150‑1 (MQ=255)
                                                          aTTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTg                                                                                           <  2:404295/150‑1 (MQ=255)
                                                            tttCTTGATGTAAAAAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGTGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGAt                                                                                         <  1:2374524/150‑1 (MQ=255)
                                                            tttCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGAt                                                                                         <  2:1176709/150‑1 (MQ=255)
                                                             ttCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATg                                                                                        <  1:169849/150‑1 (MQ=255)
                                                              tCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGt                                                                                       >  2:1010993/1‑150 (MQ=255)
                                                               cTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGAc                                                                                                          <  2:543173/130‑1 (MQ=255)
                                                                  gATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtattt                                                                                   <  1:2479045/150‑1 (MQ=255)
                                                                    tGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttac                                                                                 <  2:1517537/150‑1 (MQ=255)
                                                                      tAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttacta                                                                               >  1:363155/1‑150 (MQ=255)
                                                                      tAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttacta                                                                               <  2:175281/150‑1 (MQ=255)
                                                                       aaaCAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttactat                                                                              >  1:284328/1‑150 (MQ=255)
                                                                        aaCAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttactatt                                                                             <  2:2440245/150‑1 (MQ=255)
                                                                         aCAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttactattt                                                                            <  1:1133147/150‑1 (MQ=255)
                                                                            gagaAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTAcc                                                                         <  2:2021638/150‑1 (MQ=255)
                                                                              gaAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCTGGCATATGCCCGGGAATTATTATAGGATTATTTTCTTTGTTAAAAAAGATAAACTGTATTAACTGATGTATTTACTTTTTACCGc                                                                       <  2:822135/150‑1 (MQ=255)
                                                                              gaAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGc                                                                       <  2:311338/150‑1 (MQ=255)
                                                                                     ttATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGc                                                                >  1:1326072/1‑150 (MQ=255)
                                                                                        tCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAgg                                                             <  1:416794/150‑1 (MQ=255)
                                                                                           gTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttact                                                                                >  1:1633599/1‑128 (MQ=255)
                                                                                           gTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGtatttact                                                                                <  2:1633599/128‑1 (MQ=255)
                                                                                           gTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAggggg                                                          >  2:1201470/1‑150 (MQ=255)
                                                                                           gTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAggggg                                                          >  1:267002/1‑150 (MQ=255)
                                                                                             ccATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGTAGGGGGGa                                                        <  1:1293018/150‑1 (MQ=255)
                                                                                             ccATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGa                                                        >  1:935593/1‑150 (MQ=255)
                                                                                              cATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGAt                                                       >  2:2097316/1‑150 (MQ=255)
                                                                                              cATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGAt                                                       >  1:947151/1‑150 (MQ=255)
                                                                                                   gAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATc                                                   <  1:385649/149‑1 (MQ=255)
                                                                                                   gAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATc                                                   >  2:385649/1‑149 (MQ=255)
                                                                                                   gAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCa                                                  >  1:1505343/1‑150 (MQ=255)
                                                                                                    aaTATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAg                                                 >  1:2062855/1‑150 (MQ=255)
                                                                                                      tatGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGaa                                               <  1:220532/150‑1 (MQ=255)
                                                                                                        tGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATc                                             >  1:154274/1‑150 (MQ=255)
                                                                                                           aaaGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCcgc                                          >  1:1116387/1‑150 (MQ=255)
                                                                                                             aGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCgg                                        <  1:1015353/150‑1 (MQ=255)
                                                                                                                   ataatCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTAtta                                  >  1:1409633/1‑149 (MQ=255)
                                                                                                                    taatCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCa                                 <  2:1595263/150‑1 (MQ=255)
                                                                                                                     aatCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCat                                <  1:1333385/150‑1 (MQ=255)
                                                                                                                       tCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCatta                              <  2:2439927/150‑1 (MQ=255)
                                                                                                                        ccATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCattat                             <  2:739213/150‑1 (MQ=255)
                                                                                                                        ccATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCattat                             <  2:2363949/150‑1 (MQ=255)
                                                                                                                         cATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCattatt                            >  2:1857878/1‑150 (MQ=255)
                                                                                                                          aTTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCAttattt                           >  2:78508/1‑150 (MQ=255)
                                                                                                                           ttAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCAttattta                          <  2:71369/150‑1 (MQ=255)
                                                                                                                           ttAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCAttattta                          <  2:934979/150‑1 (MQ=255)
                                                                                                                              aaaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTtatat                       >  1:369475/1‑150 (MQ=255)
                                                                                                                                aCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTtatataa                     >  2:1491414/1‑149 (MQ=255)
                                                                                                                                    tAATCATCCCTGTCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGGGGTTATTCATTATTTATATACATTc                 <  2:1485745/150‑1 (MQ=255)
                                                                                                                                     aaTCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCa                >  1:493537/1‑150 (MQ=255)
                                                                                                                                     aaTCATCCCGGGCAAATGCCCGGGAATTAATCTAGGACTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTAACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATAAATTATATACATTCa                <  2:355848/150‑1 (MQ=255)
                                                                                                                                      atcatcCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCAGCAGGGGGGATTATCAGAATCCGAGGTTATTCATTATTTATATACATTCaa               >  2:1291867/1‑150 (MQ=255)
                                                                                                                                       tcatcCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAAt              <  2:1881851/150‑1 (MQ=255)
                                                                                                                                        catcCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATa             >  2:2483188/1‑150 (MQ=255)
                                                                                                                                        attcCCGGGCATTTGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATa             >  2:1149897/3‑150 (MQ=255)
                                                                                                                                         ttcCCGGGCATTTGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATaa            >  2:2152263/2‑150 (MQ=255)
                                                                                                                                          tcCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATaaa           >  1:455597/1‑150 (MQ=255)
                                                                                                                                           cccGGGCATTTGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATaaaa          >  1:765325/1‑150 (MQ=255)
                                                                                                                                           cccGGGCATTTGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATaaaa          >  2:1894995/1‑150 (MQ=255)
                                                                                                                                           cccGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATaaaa          >  1:798895/1‑150 (MQ=255)
                                                                                                                                            ccGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATaaaaa         >  2:488720/1‑150 (MQ=255)
                                                                                                                                              gtgCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAg       <  1:2297950/148‑1 (MQ=255)
                                                                                                                                              gggCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAg       <  1:1035667/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATATCCATAAAGCAGAGCATCTGGAAGAACGGCGATTGATGCTACAGTGGTGGGCCGATTTTCTGGATGTAAACAGAGAAAGGTTTATCAGTCCATTTGAATATGCAAAGATTAATAATCCATTAAAACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACC  >  NC_000913/4497806‑4498102

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